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1.  INTRODUCTION 


A waveform  generator  has  been  developed  for  the  National  Highway  Traffic 
Safety  Administration  (NHTSA)  that  will  provide  standard  signals  for  testing 
the  performance  characteristics  of  data  acquisition  channels  at  the  facilities 
of  NHTSA  crash  test  contractors.  Initial  development  of  the  waveform  genera- 
tor was  carried  out  under  the  Test-Site  Instrumentation  Study  (Contract  No. 
DOT-HS-8-01936 , Task  Order  No.  3).  Modifications  were  subsequently  made 
under  Phase  II  (Contract  No.  DTNH22-82-C-07041)  and  Phase  III  (Contract  No. 
DTRS-57-84-C-00003 , Task  Order  Nos.  3,  3A,  8,  8A,  and  Purchase  Order  No. 
DTRS-57-86-P-81655)  of  the  study.  This  report  presents  a description  and 
operating  manual  for  the  processing  software  package  that  was  initially 
developed  under  Phase  II  and  modified  under  Phase  III  of  this  study.  Under 
Phase  III , the  software  package  was  modified  to  account  for  the  changes  that 
have  been  made  to  the  format  of  the  waveforms  generated  by  the  waveform 
generator  instrument.  The  software  package  was  also  integrated  into  a "user 
friendly"  package.  The  software  is  implemented  on  the  VAX  11/780  minicomputer 
at  NHTSA. 

This  report  consists  of  five  sections.  The  first  section  discusses  the 
characteristics  of  the  waveforms  and  the  results  that  are  generated  by  the 
software  analysis.  The  second  section  is  a user's  manual  for  the  software 
package.  The  third  section  contains  a narrative  description  and  flow  diagrams 
of  the  software.  The  fourth  section  consists  of  a listing  of  the  software. 
The  fifth  section  describes  the  procedure  for  programming  the  test  waveforms 
into  EPROMs  (Erasable  Programmable  Read  Only  Memory) . Also  included  in  this 
section  is  a narrative  description,  a flow  diagram,  and  a source  listing  of 
the  software  which  was  developed  to  compare  the  actual  content  of  the  EPROMs 
to  idealized  test  waveforms  is  presented.  For  further  information  on  the 
waveform  generator  instrument,  the  reader  is  directed  to  the  Operator's  Manual 
for  Waveform  Generator  Model  RPG-6236-A. 
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2 . WAVEFORM  CHARACTERISTICS 


This  section  describes  the  characteristic  of  each  waveform  that  is 
generated  by  the  waveform  generator  instrument.  A discussion  of  the  results 
produced  by  the  processing  software  for  each  waveform  follows  the  discussion 
of  each  waveform  characteristic. 

2.1  RECTANGLE  WAVEFORM  CHARACTERISTICS 

The  characteristics  of  the  rectangle  waveform  include  positive  and 
negative  full-scale  amplitudes  of  200  g's  in  magnitude.  The  period  of  one 
rectangle  cycle  is  10  msec.  The  are  10  cycles  present  in  the  rectangle 

waveform  with  each  cycle  starting  at  a level  of  +200  g's  and  continuing  for  6 
msec  then  dropping  to  a level  of  =200  g's  and  continuing  for  4 msec  and 
finally  returning  to  the  starting  position  level  of  +200  g's.  (It  should  be 
noted  that  the  10th  cycle  returns  only  to  a level  of  0 g's.)  The  duration  of 
the  rectangle  waveform  is  100  msec. 

Figure  2-1  shows  a typical  rectangle  waveform  sequence. 

2.2  RECTANGLE  WAVEFORM  PROCESSING  RESULTS 

The  rectangle  waveform  is  used  to  check  the  following  characteristics  of 
a data  channel: 

a.  Time  deviation  from  theoretical  time 

b . Time  linearity 

c.  Time  offset 

d.  Steady- state  amplitude  response 

e.  Amplitude  overshoot  in  response  to  a step  input 

f.  Channel- to -channel  time  deviation 
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AMPLITUDE  (G's) 


T (MSEC) 


Figure  2-1  Rectangle  Waveform 


2.2.1  Time  Deviation  from  Theoretical  Time 


Zero  level  amplitude 
crossings 


0 


(Theoretical  times  are 
0,  6,  10,  16,  20,  etc.) 


The  actual  time  of  the  zero  crossing  of  the  recorded  rectangle  waveform 
is  measured  to  determine  the  magnitude  of  variations  in  the  data  time  base. 
The  theoretical  time  of  zero  crossing  is  established  by  the  time  reference 
waveform.  That  is,  the  first  sample  of  the  time  reference  waveform  that 
exceeds  50%  of  its  maximum  value  determines  the  value  of  time  for  which  time 
* 0.  Each  subsequent  zero  crossing  of  the  rectangle  waveform  is  supposed  to 
occur  at  6,  10,  16,  20,  26,  30,  etc.  milliseconds  after  TO  (time  -0).  The 
time  deviation  is  the  difference  between  the  measured  and  the  theoretical  zero 
crossing  times.  In  the  present  rectangle  processing  software,  the  time  of  the 
initial  rise  (at  TO)  is  not  measured.  It  is  extrapolated  to  occur  10  milli- 
seconds before  the  first  measured  zero  crossing  of  the  first  positive  going 
edge  of  the  waveform. 

2.2.2  “* 


The  rectangle  waveform  is  examined  using  the  method  of  least  squares. 
By  fitting  a straight  line  to  a number  of  X-Y  coordinates,  the  slope  and  the 
Y- intercept  of  that  line  give  an  indication  of  the  quality  of  the  data. 

The  X coordinate  used  in  this  analysis  is  the  theoretical  time  the 
waveform  begins  a complete  period  as  defined  by  a 0 g's  amplitude  level.  The 
theoretical  time  value  begins  at  zero  milliseconds  and  increments  by  10  msec 
for  each  cycle.  The  Y coordinate  used  in  this  analysis,  called  the  time 
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deviation  from  theoretical,  is  the  difference  between  the  actual  time  at 
which  the  waveform  crosses  the  0 g's  amplitude  level  and  theoretical  time  of 
the  0 g's  amplitude  level  crossing.  This  value  can  range  from  zero  in  the 
best  case  to  infinity  for  the  worst  case. 

The  slope  of  the  least  squares  straight  line  fit  is  an  indication  of 
the  time  linearity  of  the  data.  A deviation  from  zero  of  + 0.01  is  allowed 
while  acceptability  is  maintained.  This  tolerance  translates  to  + 1%  varia- 
tion in  the  data  acquisition  sample  rate  when  the  sample  rate  conforms  to  the 
minimum  requirements  for  SAE  class  1000  data  (8000  samples  per  second) . The 
ideal  case  would  be  a slope  of  zero  which  would  indicate  perfect  linearity. 

2.2.3  Time  Offset 

The  Y-axis  intercept  of  the  least  squares  straight  line  fit  is  an 
indication  of  the  time  offset.  A deviation  of  + 1.0  msec  is  allowed  while 
still  maintaining  acceptability.  The  ideal  value  for  the  Y-axis  intercept 
would  be  zero  which  would  indicate  perfect  timing  synchronization. 

2.2.4  Steady- State  Amplitude  Deviation  from  Theoretical  Amplitude 

The  steady- state  amplitude  deviation  from  theoretical  amplitude  is 
determined  by  comparing  the  steady  state  average  to  the  theoretical  value  of 
+200  g's  (depending  upon  the  sense  of  the  data)  for  each  of  the  20  occurrences 
of  the  extremes  throughout  the  entire  rectangle  waveform.  The  steady- state 
amplitude  is  calculated  by  averaging  the  data  values  over  the  closed  interval 
from  A to  B shown  below. 


5 


Point  A is  the  center  of  the  constant  amplitude  level  duration  and  point  B is 
one  millisecond  before  the  next  zero  crossing.  The  duration  of  the  averaging 
is  no  less  than  2 msec  for  the  positive  data  levels  and  1 msec  for  the 
negative  data  levels.  Deviation  from  the  theoretical  steady-state  amplitude 
is  expressed  as  a percentage  of  the  theoretical  level  and  is  determined  by 
dividing  the  difference  between  the  actual  and  theoretical  amplitudes  by  the 
value  of  200  (the  absolute  value  of  the  theoretical  amplitude) . 


2.2.5 


)litude  Overshoot  Relative  to  the  Calculated  Steady- State  Amplitude 


The  relative  overshoot  is  determined  by  obtaining  the  data  value  extreme 
which  occurs  over  the  closed  interval  from  C to  D for  the  positive  intervals 
and  from  E to  F for  the  negative  intervals. 

) 


The  value  of  relative  overshoot  is  expressed  as  a percentage  of  the  change  in 
steady- state  amplitude  by  dividing  the  calculated  overshoot  by  the  change  in 
steady- state  amplitude.  (Theoretically,  all  amplitude  changes  are  400  g's 
with  the  exception  of  the  start  of  the  first  cycle  which  rises  from  0 to  200 
g's  and  the  end  of  the  last  cycle  which  only  returns  to  0 g's  from  the  -200 
g's  amplitude  level.) 
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2.2.6  Channel- to -Channel  Time  Difference  for  the  Rectangle  Waveform 


The  channel- to -channel  time  deviation  is  a determination  of  how  each 
individual  channel  is  performing  relative  to  a common  standard.  The  value 
which  is  used  to  make  this  determination  is  the  actual  time  of  occurrence  of 
each  rectangle  cycle  beginning  at  zero  amplitude.  The  common  standard  which 
is  used  is  the  actual  times  of  occurrence  of  each  rectangle  cycle  for  the 
first  channel  processed.  The  channel- to -channel  performance  is  measured  by 
the  difference  between  the  time  of  occurrence  of  corresponding  rectangle 
cycles  for  the  first  channel  and  the  channel  being  examined. 

2.3  STAIR  WAVEFORM  CHARACTERISTICS 

The  characteristics  of  the  stair  waveform  include  amplitude  levels 
ranging  from  0 to  +200  g's  and  from  -200  to  0 g's  with  a 0 g's  level  between 
the  +200  and  -200  g's  extremes.  The  time  duration  of  each  level  is  10  msec 
and  there  are  a total  of  13  amplitude  levels  (40  g increments)  in  the  wave- 
form, of  which  three  are  the  zero  level.  The  duration  of  the  waveform  is  120 
msec;  a 5 msec  zero  level  is  added  to  the  end  of  the  waveform  to  extend  it  to 
125  msec. 

Figure  2-2  shows  a typical  stair  waveform. 

2.4  STAIR  WAVEFORM  PROCESSING  RESULTS 

The  stair  waveform  is  used  to  check  the  following  characteristics  of  the 
data  channel: 

a.  Time  deviation  from  theoretical 

b.  Time  linearity 

c.  Time  offset 

d.  Amplitude  linearity 

e.  Amplitude  offset 

f.  Steady- state  amplitude  deviation  from  theoretical  levels 

g.  Amplitude  overshoot  for  different  amplitude  levels 

h.  Channel- to -channel  time  deviation 
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AMPLITUDE  (G's) 


Figure  2-2  Stair  Waveform 
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2.4.1.  Time  Deviation  from  Theoretical 

Shown  below  is  an  example  of  the  calculation  of  time  deviation  from 
theoretical . 


The  time  deviation  from  theoretical  is  determined  by  subtracting  the  theoreti- 
cal time  the  waveform  should  reach  an  amplitude  which  is  50%  of  the  change  in 
amplitude  for  the  interval  being  considered  from  the  actual  time  the  waveform 
reaches  that  amplitude.  In  this  example,  the  theoretical  time  of  reaching  the 
amplitude  of  20  g's  (50%  of  the  amplitude  change  from  0 g's  to  40  g’s)  is  10 
msec.  In  the  STAIR  summary  report,  the  "TIME  DEVIATION  FROM  THEORETICAL  TIME" 
reported  for  each  interval  corresponds  to  the  measured  time  of  occurrence  of 
each  numbered  rise  shown  in  Figure  2-2. 

2.4.2  Time  Linearity 

The  stair  waveform  is  examined  for  time  linearity  and  offset  using  the 
method  of  least  squares.  When  a straight  line  is  fitted  through  a number  of 
X-Y  coordinate  pairs,  the  slope  and  Y- intercept  of  the  line  give  an  indication 
of  the  quality  of  the  data.  The  X-coordinates  used  are  the  theoretical  times 
at  which  the  waveform  changes  amplitude  levels.  The  values  for  the  X-coor- 
dinate  range  from  10  to  120  msec  in  increments  of  10  msec.  The  Y- coordinates 
used  in  the  analysis  are  the  time  deviation  from  theoretical  values  determined 
for  the  waveform  as  described  above.  The  values  can  range  from  zero  for  the 
best  case  to  infinity  for  the  worst  case. 

The  slope  of  the  least  squares  fit  straight  line  is  an  indication  of  the 
time  linearity  of  the  data.  The  value  for  the  slope  is  allowed  a tolerance  of 
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10  msec 
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+0.01.  This  corresponds  to  an  allowable  deviation  of  +1%.  The  ideal  case 
would  have  a slope  of  zero,  indicating  perfect  linearity. 

2.4.3  Time  Offset 

The  Y=  intercept  of  the  least  squares  fit  straight  line  is  an  indication 
of  the  time  offset.  The  tolerance  for  this  value  is  +1.0  msec.  The  ideal 
time  offset  would  be  zero  msec  indicating  perfect  timing. 


2.4.4 


The  amplitude  linearity  is  determined  using  the  same  approach  that  is 
used  for  the  time  linearity  calculation  described  in  Section  2.4.2.  The 
tolerance  for  amplitude  linearity  is  +2.5%. 


2.4.5 


The  amplitude  offset  is  determined  using  a similar  approach  to  that  used 
for  the  time  offset  calculation  described  in  Section  2.4.3.  The  tolerance  for 


amplitude  offset  is  +2.5%. 


DEVIATION 

FROM 

THEORETICAL 


IDEAL  SLOPE  - 0 
IDEAL  INTERCEPT 


0 


2.4.6 


State  Amplitude  Deviation  from  Theoretical  Amplitude 
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The  data  used  to  determine  the  steady-state  amplitude  occurs  in  the 
closed  interval  from  L to  M shown  above.  L corresponds  to  the  center  of  the 
constant  amplitude  levels  and  M is  a point  which  occurs  one  millisecond  before 
the  next  rise  in  amplitude.  The  steady- state  amplitude  is  calculated  by 
taking  the  average  of  the  amplitude  values  in  this  interval.  Once  the  actual 
steady- state  value  is  determined,  it  is  then  compared  to  the  theoretical 
steady-state  value.  The  difference  between  the  theoretical  and  actual  values 
is  divided  by  the  full-scale  amplitude  (200  g's)  in  order  to  express  the 
deviation  from  theoretical  amplitude  as  a percentage  of  the  full-scale 
amplitude.  The  tolerance  for  this  portion  of  the  analysis  is  +2,5%.  In  the 
STAIR  summary  report,  the  "STEADY  STATE  AMPLITUDE  DEVIATION  FROM  THEORETICAL 
AMPLITUDE"  reported  for  each  interval  corresponds  to  the  percentage  dif- 
ference, described  above,  for  each  numbered  level  shown  in  Figure  2-2. 

2.4.7  Amplitude  Overshoot  Relative  to  Calculated  Steady- State  Value 

The  maximum  data  point  in  the  closed  interval  from  P to  Q shown  below  is 
used  to  determine  the  relative  amplitude  overshoot  by  comparing  the  value  of 
the  point  to  the  calculated  (average)  steady- state  amplitude.  The  difference 
between  the  maximum  and  steady- state  values  is  divided  by  the  difference 
between  the  steady- state  levels  occurring  before  and  after  the  step  and 
expressed  as  a percentage  of  the  step  height.  Overshoot  of  +13%  is  considered 
acceptable.  In  the  STAIR  summary  report,  the  "AMPLITUDE  OVERSHOOT  RELATIVE  TO 
CALCULATED  STEADY  STATE  AMPLITUDE"  reported  for  each  interval  corresponds  to 
the  overshoot  calculated  for  each  numbered  rise  shown  in  Figure  2-2. 


2.4.8  Channel- to -Channel  Time  Difference  for  the  Stair  Waveform 


The  channel- to -channel  time  deviation  is  an  indication  of  how  each 
individual  channel  is  performing  relative  to  a common  standard.  The  value 
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which  is  examined  is  the  actual  time  of  occurrence  of  each  amplitude  change 
during  the  stair  waveform.  Each  value  is  compared  to  the  actual  time  of 
occurrence  of  the  corresponding  amplitude  change  of  the  waveform  that  is 
contained  in  the  first  channel  processed  thus  making  the  first  channel  which 
is  processed  the  standard.  The  tolerance  for  this  part  of  the  analysis  is 
+0.25  msec. 

2.5  HALF- SINE  WAVEFORM  CHARACTERISTICS 

The  half- sine  waveform  is  broken  into  two  components  - an  X-axis  and  a 
Z-axis.  The  two  components  are  produced  on  separate  channels  of  the  waveform 
generator  instrument.  (This  forces  the  requirement  of  processing  the  data  in 
multiples  of  two  channels  per  execution  of  the  processing  software.)  Both 
components  consist  of  half-sine  pulses  which  have  a time  duration  of  2.985 
msec.  The  X-axis  component  has  a peak  amplitude  of  -200  g's  and  the  Z-axis 
component  has  a peak  amplitude  of  +200  g's.  Together,  these  two  components 
form  a resultant  waveform  which,  when  analyzed  in  terms  of  HIC  (Head  Injury 
Criterion),  give  a HIC  value  of  1665.6  with  a HIC  calculation  time  duration  of 
2.0  msec.  It  should  be  noted  that  the  components  are  synchronized  in  time. 
The  half- sine  waveform  was  designed  to  simulate  a short  duration  impact  which 
is  often  found  in  actual  crash  test  data.  The  half-sine  waveform  components 
are  used  to  provide  a means  of  checking  the  HIC  computation  accuracy  and  the 
channel  to  channel  time  difference  between  the  two  channels  containing  the 
waveform  components.  Figure  2-3  shows  the  two  components  of  the  half- sine 
waveform. 

2.6  HALF-SINE  WAVEFORM  PROCESSING  RESULTS 
2.6.1  Channel -to -Channel  Time  Difference 

The  actual  times  for  the  peak  amplitude  of  the  X and  Z components 
are  utilized  to  determine  the  channel- to -channel  time  difference.  The  times 
of  the  peak  amplitudes  are  determined  by  calculating  the  time  at  which  each 
waveform  crosses  an  amplitude  level  which  is  20%  of  full  scale.  It  is  then 
assumed  that  the  peak  amplitude  occurs  at  the  time  which  is  exactly  halfway 
between  the  two  20%  level  crossings.  The  difference  between  the  X-axis  and 
Z-axis  time-of-peak  amplitude  values  is  the  channel- to -channel  time  dif- 
ference. The  tolerance  allowed  for  this  difference  is  +0.1  msec. 
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Figure  2-3  Half -Sine  Waveform 
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2.6.2  Waveform  Component  Time  Deviation 


To  determine  the  time  deviation  for  the  two  components  of  the  waveform, 
the  times  of  peak  amplitude  occurrence  calculated  for  the  channel- to -channel 
time  difference  (described  above)  are  compared  against  theoretical  times  at 
which  the  peaks  should  occur.  The  allowable  deviation  for  each  component  is 
+1.0  msec. 

2.6.3  HIC  Deviation 

The  HIC  deviation  is  determined  by  comparing  the  HIC  value  determined 
from  the  resultant  of  the  recorded  X and  Z components  against  the  theoretical 
value.  The  theoretical  value  in  this  case  may  not  be  the  value  1665.6 
mentioned  above  since  the  theoretical  is  determined  from  curves  which  are 
generated  so  that  they  have  the  same  amplitude  and  time  duration  as  the 
recorded  waveform  components.  To  accomplish  this,  true  half-sine  curves 
(a  - A * sin(wt) ) are  fitted  to  the  peak  amplitude  and  time  values  calculated 
during  the  channel- to -channel  time  difference  determination  (Section  2.6.1). 
This  method  is  used  in  order  to  separate  the  errors  in  time  and  amplitude 
deviation  from  the  HIC  calculation.  (If  the  waveform  components  do  not  match 
the  time  duration  and/or  the  amplitude  of  the  theoretical  waveform  components , 
then  the  • chances  that  the  HIC  calculation  will  fall  within  the  +6%  limit  set 
for  the  HIC  calculation  tolerance  are  greatly  reduced.) 

2.7  CRASH  WAVEFORM  CHARACTERISTICS 

As  with  the  half-sine  waveform,  the  crash  waveform  is  composed  of  two 
components:  X-axis  and  Z-axis.  The  waveform  was  obtained  from  the  head 
accelerations  which  occurred  in  a crash  test  which  was  assigned  test  number 
206  by  NHTSA.  The  resultant  of  the  X and  Z components  yield  a computed  HIC 
value  of  930.  The  crash  waveform  is  used  as  a means  of  checking  the  accuracy 
of  the  data  channels  with  respect  to  actual  crash  test  data.  Figure  2-4  shows 
the  X and  Z components  of  the  crash  waveform. 
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Figure  2-4  Crash  Waveform 


2.8  CRASH  WAVEFORM  PROCESSING  WAVEFORM  RESULTS 


2.8.1  HIC  Deviation 

HIC  value  deviation  is  determined  by  comparing  the  HIC  value  obtained  by 
processing  the  resultant  of  the  recorded  X and  Z components  against  the 
theoretical  HIC  value  of  930.  As  with  the  half -sine  HIC  deviation  analysis, 
the  allowable  difference  is  ±6%. 

2 . 9 SUM-QF-SINES  WAVEFORM  CHARACTERISTICS 

The  sum-of-sines  waveform  is  a test  signal  which  is  used  to  determine  the 
amplitude  frequency  response  of  a data  acquisition  channel. 

The  sum-of-sines  waveform  consists  of  the  sum  of  equal  amplitude  sine 
waves  whose  frequencies  are  successive  integer  harmonics  of  the  fundamental 
(lowest)  frequency  of  the  composite  signal.  For  the  Phase  III  waveform 
generator  instrument,  the  sum-of-sines  composite  waveform  consists  of  the  sum 
of  the  first  14  harmonics  of  273.4375  Hz.  The  peak  level  of  the  composite 
output  is  scaled  to  be  the  full-scale  input  of  the  data  acquisition  channel 

being  tested.  The  sum-of-sines  composite  produced  by  the  Phase  III  waveform 

« 

generator  instrument  is  shown  in  Figure  2-5. 

2.10  SUM-OF-SINES  WAVEFORM  PROCESSING 

The  recorded  data  acquisition  channel  response  to  the  sum-of-sines 
waveform  is  processed  to  determine  the  amplitude  frequency  response  of  the 
channel  at  the  test  signal  frequencies. 

This  data  processing  produces  three  types  of  output  files:  plot  files;  a 
statistical  data  file;  and  a summary  file.  The  plot  file  consists  of 
fourteen  ordered  data  pairs  for  plotting  signal  frequency  vs.  amplitude  ratio 
in  dB.  The  statistical  data  file  consists  of  sequences  of  fourteen  values  of 
amplitude  ratio  (one  value  at  each  signal  frequency) , (one  sequence  for  each 
data  channel  analyzed).  For  each  signal  frequency,  the  summary  file 
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provides  the  signal  frequency,  the  DFT  linet  number,  the  DFT  output  frequency, 
the  input  DFT  value,  the  output  DFT  value,  the  scaled  DFT  ratio,  and  the 
amplitude  ratio  in  dB.  One  table,  of  this  description,  is  provided  for  each 
data  channel  tested. 

It  is  appropriate  to  discuss  the  information  provided  in  the  summary 
file.  The  discrete  Fourier  transform  (DFT)  plays  a central  role  in  the 
determination  of  amplitude  frequency  response.  The  DFT  output  consists  of  a 
sequence  of  numbers.  As  the  DFT  is  used  here,  each  number  in  the  DFT  output 
can  be  considered  the  energy  level  in  the  signal  analyzed,  at  a particular 
frequency.  The  frequency,  to  which  each  DFT  output  value  corresponds,  starts 
at  zero  for  the  first  DFT  output  value  and  increases  monotonically , with  a 
constant  increment,  for  each  subsequent  member  of  the  DFT  output  sequence. 
The  "DFT  line  number"  is  the  number  of  frequency  increments  contained  in  the 
"DFT  output  frequency"  closest  to  the  "signal  frequency"  listed.  The  input 
DFT  value  is  the  magnitude  of  the  output  value  of  the  DFT  (of  the  input 
waveform)  whose  corresponding  frequency  is  closest  to  a signal  frequency.  The 
output  DFT  value  is  the  magnitude  of  the  output  value  of  the  DFT  (of  the 
channel  output  signal)  whose  corresponding  frequency  is  closest  to  a signal 
frequency.  The  DFT  ratio  is  the  ratio  of  the  output  DFT  value  to  the  input 
DFT  value  at  the  same  signal  frequency.  This  ratio  is  multiplied  by  a scale 
factor,  which  is  0.993  multiplied  by  the  ratio  of  the  input  DFT  value  to  the 
output  DFT  value  at  the  lowest  signal  frequency.  The  scaled  DFT  ratio  is  the 
numeric  amplitude  ratio. 
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3.  USER'S  MANUAL 


This  manual  presents  the  procedure  for  processing  the  recorded  waveform 
generator  instrument  output  data.  This  software  package  performs  two  general 
functions.  The  first  function  is  to  prepare  the  channel  output  data  files  so 
that  individual  waveforms  can  be  processed.  The  second  function  is  to 
process  the  data.  This  software  package  is  capable  of  processing  from  2 to  18 
channels  in  a single  execution  pass.  (Note:  If  half- sine  and  crash  waveforms 
are  to  be  processed,  the  sequence  of  data  files  to  be  processed  must  alternate 
between  X (group  1 waveforms)  and  Z (group  2 waveforms)  axis  data,  starting 
with  an  X-axis  data  file.  ) 

Actual  data  recorded  from  the  output  of  a waveform  generator  instrument 
will  be  used  for  examples  throughout  this  manual.  The  data  was  recorded  at 
Calspan  Corporation,  Buffalo,  NY.  The  data  files  which  will  be  utilized  are: 


UDS$WORK: [S000]S0002DA00.015  Time  synchronization 

UDS$WORK: [S000]S0002AA00.012  Channel  12  (X-axis) 

UDS$WORK: [S000]S0002AA00.014  Channel  14  (Z-axis) 

The  examples  to  follow  will  be  for  data  processed  on  the  NHTSA  VAX 
11/780. 

The  remainder  of  this  user's  manual  will  assume  that  the  reader  has  a 
working  knowledge  of  the  use  of  the  NHTSA  VAX  11/780  computer  and  the  Digital 
Command  Language  (DCL) . Program  prompts  will  be  shown  as  bold  characters  in 
this  manual  and  user  responses  shown  will  be  underlined. 

3.1  PROGRAM  EXECUTION 

To  begin  the  processing  procedure  the  user  must  type: 

RUN  WGPROC 
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The  program  will  then  respond  with: 


***  WAVEFORM  GENERATOR  PROCESSING  SOFTWARE  *** 


XO-MAR-88  11:27:51 


DO  YOU  WANT  ANALYSIS  OF  HALF- SINE  AND  CRASH  PULSES  IN  THE  PROCESSING?  (Y  or  N) 

»Y 


OPTION  1 : FINDING  DATA  POINT  INDEX  NUMBERS  CORRESPONDING  TO  EDGES 
OPTION  2 : ENTER  CORRECT  DATA  POINT  INDEX  NUMBERS 
ENTER  DESIRED  OPTION  NUMBER  (1/2) 

»> 


Figure  4-1  (see  Section  4)  presents  a flow  diagram  of  the  processing 
package.  Not  shown  in  Figure  4-1  is  the  first  option  provided  to  the  user, 
the  selection  of  full  or  partial  processing.  The  option  of  eliminating  the 
half- sine  and  crash  pulse  processing  is  provided  so  that  when  the  appropriate 
group  1 and  group  2 signals  are  not  used  as  inputs  to  the  data  acquisition 
system  under  test,  only  the  meaningful  processing  (RECTANGLE,  STAIR,  AND  SUM- 
OF-SINES)  will  be  performed. 

Option  1 takes  the  time  synchronization  file  and  detects  the  rising  or 
falling  edge  data  point  index  numbers.  If  there  are  more  than  four  edges  for 
the  time  zero  and  inverted  time  zero  data  file  types  or  more  than  two  edges 
for  the  delayed  time  zero  and  inverted  delayed  time  zero  data  file  types , the 
program  will  warn  the  user  of  the  occurrence  of  the  extra  edges  and  end 
program  execution.  The  occurrence  of  extra  edges  stems  primarily  from  two 
sources . 

These  sources  are: 

1)  The  waveform  generator  instrument  produces  a pre-trigger  pulse  prior 
to  the  output  of  the  actual  test  waveforms  sequence.  If  the  time 
synchronization  channel  on  the  data  acquisition  system  records  the 
pre- trigger  pulse,  an  extra  edge  will  result. 
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2)  Noise  on  the  time  synchronization  channel  on  the  data  acquisition 
system  can  be  of  sufficient  magnitude  such  that  it  is  mistaken  for 
extra  edges . 

In  case  of  an  occurrence  of  extra  edges,  the  user  should  obtain  a plot  of 
the  time  synchronization  file  in  order  to  ascertain  the  indices  of  the  correct 
edges.  (This  can  be  done  by  using  the  "VTEKPLT"  program  which  resides  in  the 
[ASGPROG]  account.)  Upon  determining  which  of  the  indices  represent  the 
correct  edges,  the  user  can  then  re-execute  the  program  and  make  use  of 
Option  2.  An  example  of  the  occurrence  of  extra  edges  is  given  by  the  data 
chosen  for  use  in  the  following  sample  processing  run. 

Starting  from  the  initial  execution  of  the  processing  program  (assuming 
the  occurrence  of  extra  edges  is  not  known),  the  user  enters: 

RUN  WGPROC 

***  WAVEFORM  GENERATOR  PROCESSING  SOFTWARE  *** 

10-MAR-88  11:27:51 

DO  YOU  WANT  ANALYSIS  OF  HALF- SINE  AND  CRASH  PULSES  IN  THE  PROCESSING?  (Y  or  N) 
»Y 


OPTION  1 : FINDING  DATA  POINT  INDEX  NUMBERS  CORRESPONDING  TO  EDGES 
OPTION  2 : ENTER  CORRECT  DATA  POINT  INDEX  NUMBERS 
ENTER  DESIRED  OPTION  NUMBER  (1/2) 

»>  1 

ENTER  TIME  SYNCHRONIZATION  FILE  NAME 

»UDS$WORK:  rSQOOl  S0002DA00 . 015 
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WARNING! 

ADDITIONAL  EDGES  IN  TIME  ZERO  FILE  HAVE  BEEN  DETECTED 
OBTAIN  A PLOT  OF  TIME  ZERO  FILE 

AND  USE  OPTION  2 TO  ENTER  THE  CORRECT  DATA  POINT  NUMBERS 


END  OF  PROCESSING  RUN 

Upon  running  the  program  "VTEKPLT88  in  the  [ASGPROG]  account  and  obtaining 
a plot  of  the  time  synchronization  file,  it  becomes  clear  that  the  first  edge 
is  not  an  indicator  for  the  leading  edge  of  recorded  data  (see  Figure  2-2) . 
Knowing  this,  the  user  should  then  proceed  as  shown  below: 

RUN  WGPROC 

The  program  will  then  respond  with: 

***  WAVEFORM  GENERATOR  PROCESSING  SOFTWARE  *** 


XO-MAR-88 


11:29:22 


DO  YOU  WANT  ANALYSIS  OF  HALF- SINE  AND  CRASH  PULSES  IN  THE  PROCESSING?  (Y  or  N) 

»Y 
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OPTION  1 : FINDING  DATA  POINT  INDEX  NUMBERS  CORRESPONDING  TO  EDGES 
OPTION  2 : ENTER  CORRECT  DATA  POINT  INDEX  NUMBERS 
ENTER  DESIRED  OPTION  NUMBER  (1/2) 

»>  2 

ENTER  1 IF  TIMING  FILE  IS  A STANDARD  TIME  ZERO  FILE 

2 IF  TIMING  FILE  IS  A DELAYED  TIME  ZERO  FILE 

3 IF  TIMING  FILE  IS  AN  INVERTED  TIME  ZERO  FILE 

4 IF  TIMING  FILE  IS  A DELAYED  & INVERTED  TIME  ZERO  FILE 

1 

ENTER  DATA  POINT  INDEX  NUMBER  1 
»>  3 

ENTER  DATA  POINT  INDEX  NUMBER  2 
»>  1676 

ENTER  DATA  POINT  INDEX  NUMBER  3 
»>  3350 

ENTER  DATA  POINT  INDEX  NUMBER  4 

»>  5024 

**  BEGINNING  OF  INDIVIDUAL  WAVEFORM  SEPARATION  ** 


NOTE:  THE  TOTAL  NUMBER  OF  FILE  NAMES  YOU  ENTER 

MUST  BE  A MULTIPLE  OF  TWO.  THE  FIRST,  THIRD, 
FIFTH,  etc.,  FILE  MUST  BE  AN  X- COMPONENT  FILE 
AND  THE  SECOND,  FOURTH,  SIXTH,  etc.,  FILE 
MUST  BE  A Z-COMPONENT  FILE. 
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ENTER  DATA  FILE  NAME  (XXXXXXXXX.XXX) 

IF  NO  MORE  FILES,  TYPE  "CNTR  Z" 

»UDSSWORK:  [ S000 1 S0002AA00 . 012 

DEL-  7.4686315E-05 

ENTER  DATA  FILE  NAME  (X300KXX30LXXX) 

IF  NO  MORE  FILES,  TYPE  "CNTR  ZB 
»UDS$WORK;  f S00Q 1 S0Q02AA00 . 014 

DEL-  7.4686315E-05 

ENTER  DATA  FILE  NAME  (XXXXXXXXX.XXX) 

IF  NO  MORE  FILES,  TYPE  "CNTR  Z" 

» Exit  (system  response  to  a 


**-*  END  of  PHASE  1 *** 

TEST  ( 1)  - 0002 

TEST  ( 2)  - 0002 

PROCESSING  SEQUENCE: 

• RECTANGLE  WAVEFORM 
STAIR  WAVEFORM 
HALF- SINE  WAVEFORM 
CRASH  WAVEFORM 
SUM- OF -SINE  WAVEFORM 

**  RECTANGLE  WAVEFORM  PROCESSOR  ** 
PROCESSING  COMPLETED  FOR  FILE  S0002AA00 . F12 
PROCESSING  COMPLETED  FOR  FILE  S0002AA00 . F14 

**  STAIR  WAVEFORM  PROCESSOR  ** 

PROCESSING  COMPLETED  FOR  FILE  S0002AA00 . E12 
PROCESSING  COMPLETED  FOR  FILE  S0002AA00 . E14 


"CNTR  Z") 
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STEADY- STATE  AMPLITUDE  DEVIATIONS  FROM  BEST  FIT  STRAIGHT  LINE  (%  FULL  SCALE) 


DEV ( 2)- 

-0.184 

DEV(  3)- 

0.008 

DEV(  4)- 

0.057 

DEV( 

5)  = 

0.144 

DEV(  6)- 

0.014 

DEV ( 7)- 

0.095 

DEV(  8)= 

0.193 

DEV  ( 

9)  = 

0.140 

DEV(IO)- 

-0.155 

DEV(ll)— 

-0.095 

DEV(12)=* 

-0.217 

DEV( 

STEADY- STATE  AMPLITUDE  DEVIATIONS  FROM  BEST  FIT 

STRAIGHT  LINE  (% 

FULL  SCALE) 

DEV(  2)- 

-0.076 

DEV ( 3)- 

-0.004 

DEV ( 4)- 

0.210 

DEV( 

5)- 

0.034 

DEV(  6)- 

-0.063 

DEV ( 7)- 

0.148 

DEV(  8)- 

0.159 

DEV( 

9)- 

0.116 

DEV(IO)- 

-0.098 

DEV(ll)- 

-0.168 

DEV(12)- 

-0.258 

DEV( 

**  HALF- 

SINE  WAVEFORM  PROCESSOR  ** 

RESULT  FILE  FOR  SET  NUMBER  1 

CALCULATING  *HIC*  FOR  SET  NUMBER  1 

PROCESSING  SET  NUMBER  1 

**  CRASH  WAVEFORM  PROCESSOR  ** 

RESULT  FILE  FOR  SET  NUMBER  1 

CALCULATING  *HIC*  FOR  SET  NUMBER  1 

PROCESSING  SET  NUMBER  1 

**  SUM- OF- SINE  WAVEFORM  PROCESSOR  ** 

PROCESSING  COMPLETED  FOR  FILE  S0002AA00 . V12 

PROCESSING  COMPLETED  FOR  FILE  S0002AA00 .V14 

END  OF  PROCESSING  RUN 

As  indicated  above,  once  the  user  specifies  all  of  the  files  to  be 
processed,  there  is  no  further  interaction  with  the  user.  One  point  should  be 
noted.  For  the  processing  of  the  rectangular,  staircase  and  sum-of-sines 
waveforms,  the  program  indicates  that  processing  is  completed  for  files  which 
the  user  did  not  specify.  The  file  names  which  are  indicated  are  hold-overs 
from  a previous  version  of  the  processing  software.  In  the  previous  version, 
each  waveform  of  each  channel  was  separated  from  the  others  of  the  channel 
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and  placed  into  a separate  file.  This  feature  was  removed  from  the  software 
in  order  to  reduce  the  number  of  calls  to  the  system  storage  devices  and  the 
processing  time.  The  file  names  no  longer  have  any  physical  meaning  outside 
of  the  program.  The  use  of  the  file  names  has  been  retained  in  order  to  allow 
the  user  to  track  the  progress  of  the  program.  These  "files"  can  now  be 
thought  of  as  "pseudo = files . n 

3.2  SUMMARY  REPORTS  AND  STATISTICAL  DATA  FILES 

During  the  execution  of  the  program,  a number  of  files  are  generated. 
The  files  which  are  generated  can  be  divided  into  two  groups.  The  first  group 
is  the  summary  files.  The  summary  files  are  files  which  report  on  the 
quality  of  the  data  in  terms  of  the  "properties"  investigated  in  a clear 
formatted  manner.  The  second  group  are  the  statistical  data  files  which 
contain  information  which  can  later  be  utilized  in  other  ways  such  as  in 
statistical  programs.  (The  other  statistical  programs  are  not  part  of  the 
signal  processing  software  package.)  These  statistical  data  files  are 
generally  not  in  a format  which  offers  clearly  meaningful  data  just  by  visual 
inspection.  Unless  the  user  intends  to  conduct  a statistical  analysis  on  a 
number  of  program  runs,  all  questions  relating  to  the  quality  of  the  recorded 
data  can  be  answered  from  the  summary  files. 

In  Table  3 = 1,  the  names  of  the  summary  and  statistical  data  files 
generated  during  program  execution  are  presented.  Also  presented  in  this 
table  is  a brief  description  of  the  contents  of  each  of  the  files.  Examples 
of  these  files  are  presented  following  the  table.  It  should  be  noted  that  in 
order  to  obtain  easily  readable  summary  reports,  the  user  must  have  access  to 
a device  (terminal  or  printer)  which  can  operate  in  132-column  mode. 
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Table  3-1  Files  Generated  During  the  Execution  of  WGPROC 


Waveform 

Rectangle 


Staircase 


Half- sine 
Crash 

Sum- of -sines 

Half -sine 
and  Crash 

Note : 


File  name 


Description 


REC####.SUM 
REC####. 001 

REC####. 002 

REC####. 003 

REC####. 004 


Summary  report 

Data  file  containing  time  deviation  from 
theoretical  time 

Data  file  containing  time  linearity  slope  and 
time  offset  intercept 

Data  file  containing  steady- state  amplitude 
deviation  from  theoretical  amplitude 
Data  file  containing  amplitude  overshoot 
relative  to  calculated  steady- state  amplitude 


STAIR####. SUM 
STAIR####. 001 

STAIR#### . 002 

STAIR#### . 003 

STAIR#### . 004 

STAIR####. 005 


Summary  report 

Data  file  containing  time  deviation  from 
theoretical  time 

Data  file  containing  time  linearity  slope  and 
time  offset  intercept 

Data  file  containing  amplitude  linearity 

slope  and  amplitude  offset  intercept 

Data  file  containing  steady- state  amplitude 

deviation  from  theoretical  amplitude 

Data  file  containing  amplitude  overshoot 

relative  to  calculated  steady- state  amplitude 


HSNE#### . SUM 
HSNE####. 001 
HSNE####. 002 


Summary  report 

Data  file  containing  HIC  deviation 

Data  file  containing  time  deviation  of  X and 

Z components  from  theoretical 


CRASH#### . SUM  Summary  report 

CRASH####. 001  Data  file  containing  HIC  deviation 


SUMSN#### . SUM 
SUMSN#### . 001 

FR####.Vnn 


Summary  report 

Data  file  containing  frequency  of  actual 
signal  and  amplitude  ratio 
Frequency  response  plot  data  files 


RESULT. SUM  File  listing  waveform  component  combinations 

used  in  the  calculation  of  the  X-Z  resultants 


####  is  the  test  number  assigned  to  the  data  by  NHTSA.  This 
corresponds  to  the  four- digit  sequences  displayed  after  the 
"END  OF  PHASE  1"  message  given  during  the  execution  of 
WGPROC.  nn  is  the  channel  number  of  the  data.  This  cor- 
responds to  the  last  two  digits  in  the  extension  of  the 
input  data  file  names. 
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RECTANGLE  WAVEFORM  PROCESSOR  DATA  SUMMARY  •* 
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Figure  3-3  RECTANGLE  WAVEFORM 
Sample  Summary  Report  and  Statistical  Data  Files 
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AMPLITUDE  OVERSHOOT  RELATIVE  TO  CALCULATED  STEADY-STATE  AMPLITUDE  (ALLOWABLE  LIMIT  +/-  13  X) 

INTERVAL  NUMBER 
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Figure  3-3  continued 
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AMPLITUDE  OVERSHOOT  RELATIVE  TO  CALCULATED  STEADY-STATE  AMPLITUDE  (ALLOWABLE  LIMIT  +/-  13  X) 

INTERVAL  NUMBER 
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Figure  3=3  continued 
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INEIS  PROCESSED 
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Figure  3-3  continued 
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A VALUE  OF  +/“  *0®®@  INDICATES  THAT  THE  QUANTITY  COULD . NOT  BE  DETERMINED 
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Figure  3-4  continued 
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Figure  3-5  HALF -SINE  WAVEFORM 
Sample  Summary  Report  and  Statistical  Data  Files 


35 


HSNE0002.002 


I 


I 


> cs;  < 


Figure  3-6  CRASH  WAVEFORM 
Sample  Summary  Report  and  Statistical  Data  Files 


36 


10-MAR— 88  11:32:26 


« 

* 

« 


CM 


Ul 

§ 

h-  — tn 
-!  *-  Q 


I 


8-2 
-IU.I- 
< Q < 

8 “ 


ION*® 

© in  cm  in 


© © in  m- 
T”  (O  CM  00 

an 


N ® ^ O)  *“  © 


cm  co  o»  m r*.  eo 

CM  CM  CM 


I I I I I I 


SMUNNSSooNnnsn 
©mo>iN®®minoiin®m  — oi 
ncocMr>*o»“<o»-inioioco^in 

O>mf^K5Q0'M”'M-tOCS'M-®UJ'M-rO 

aaaaoMoncM^«-ss« 

Q © © © © © © © Q Q © © © Q 


CM 

h-  ui 

3H3 

> 

Q.  U 1 

1-  Q < 

CM 
CM  ® 
»■  ® 
> ® 
• Q£ 
® U. 

s 

3 > 

3 

t~  UI 

CM 

3 1-3 

® 

Q.  U.  _J 

® 

Z O < 

® Ui 

in  ~t 

— > 

«eoo>io»-inin^®in»-cM»-<o 

^oooo®cMiniOh>n®vCM^cM 

®r^r>*in'M-r«»^'M,fO»-r^®cM'M‘ 

O>CMC0aOI^CM©©<O9CM©O)S 
minr^®®rvm»-®roioin®» 
CM»-»-»-OMOCMOOiniOCM*- 
CM  CM  CM  CM  *•»-*- 


«rv®®®^®^®«^incM®m 

~-mCM©t^*'0>*-«f*«.inCM©»- 

■^-r^®-M-®-M-®®CMr»K>®r'Uj 

^■®in^'«-inior«or»cMininv-in 

cointncocococoinincoco©©© 


Lk 


ui 

in  >»  o)  oi  ®®o>©o>®ojoj©F*f^ 


2 

i- 

2 

m 

®®cMMn®«n®m»-r^ 

UI 

o 

tz 

3 

oor^cM®o»fv<oM=c>®<oinM’eM 

Q. 

Q»  Ui 

m 

cm  oo  n ® m 

Mn®w 

&4o 

in 

Ul 

a 

1- 

3 

3 

05 

r>* »-  <o  ® m 

® 

© 

CM 

*-  « 0 

Q£ 

© 

Ul 

r* 

«-  at  to  ■*  »- 

05 

m 

m 

® 

co  in 

m 

K 

(K 

CM 

m 

co  ® in  co  ©s 

® 

cm  in 

© 

2 

> 

U. 

*-*-*“*- 

CM 

CM 

CM 

m 

n in 

m 

h 

U 

t- 

Ui 

# 

CM 

■M- 

inMn-Ki 

in 

<o 

« 

© 

CM  -M- 

to 

-J 

3 

u. 

z 

o 

M- 

ra 

cm  co  ® in  ® 

K> 

to 

® M- 

00 

3 

o 

M 

z 

— — CM  CN  CM 

m 

m 

■M- 

M- 

■n  in 

m 

in 

UJ 

_l 

UI 

£g 

0£ 

U. 

> 

m 

® 

m ® m ® m 

s 

in 

® 

in 

s m 

© 

-1 

r>» 

m 

CMOMnCM 

® 

m 

CM 

® 

m 

2 

✓“V 

M 

m oo  cm  co 

® 

K) 

r* 

in  oo 

CM 

UJ 

N 

•M- 

© 

n t-  © s 

m 

Oi 

m 

00 

cm  to 

o 

3 

X 

• 

• 

• • • • • 

• 

• 

• 

• 

• • 

• 

IO 

o 

® n iv  ® m- 

h» 

® 

■M- 

1^ 

— M> 

© 

in 

Ui 

r» 

■M- 

CM  05  ® M"  -- 

© 

to 

m 

© 

ao  in 

CM 

E 

CM 

in 

© ® in  © o» 

-r 

r». 

s 

cm  in 

© 

CM 

CM 

CM 

m 

m m 

m 

Figure  3-7  SUM-0F- SINES  WAVEFORM 
Sample  Summary  Report  and  Statistical  Data  Files 


37 


RESULT  FOR  FILE  S®@@2AA®@. VI 4 

FREQUENCY  RESPONSE  FILE  FR@@@2„V14 


I 


S 

fco.  e=aj  00 

**=3  fr==  © 


Sn=2 

JU.N" 

3°3 


t=  kJ 

£E3 


tz§£ 

“ t-°  r 

3i 


©P'»P5PS®©0»©©PSO4(^'*>'* 
®04  04©^p«»'*>oj®©*“04©p» 

co*oo«o*e*o*oo 


s e 


OB  06  04  04 

III" 


I 


<4=pj<<>«©®  ® «*• 
99  n ® ^ © 05  as 


©©P^PoOJPS'*. 

as  as as  as  as  r»  © n c4 


«=>  © i 

© as  ® ( 


PJ©©©f***»®®.>~ 

NMOXOSeQ'tW 

'*^s=r 


® <Sf  <o  cy 
sesoioffi^ 


as®®^©^*”©^"*’ 


04  I 
06  • 

NNNN»=' 


9&P»<®©S©<*>®<=-©*0©04eSP} 

®=p'jcy®fe»'#©i'*=eoi^>m<y«©^ 

‘*”F*®<>®<4=®®e4f»*n®6^© 

<*@i©»-<*=©PSP*P»04©©^© 

00000000*0000* 

ieuinwtetewifiisteeeceio 

no  <30  as  SB  CXK>  ro  •=*=  CT=>  03  030  cno  <=Z=>  (lO  UIE 


©v~®©04r*>na»®nat©<?»p* 

©f**o4®asp‘©'*’as©©©«oo« 

©*-04®©a&©«-’'^l*»P50S©<5- 

^®|v^.®®io®®tyt^<j=s©<^ 

08©®®ps©as<.“®p*®04©® 


oo  ee 
O © 
o*  to 


eo  oo 
CM  CV 

I I 

s cy 


T-  © 

(^  © 
04  04 
i I 

o”  n 

<*  ® 
as  r» 
<*  04 

P5  PS 

7*1* 

%% 

p*  *— 

os « 

77 

gg 

f<*  04 
® © 
• * 

© © 

77 

« « 
»-  ® 
04  ■*■ 
m m 

04  04 

77 

■*=  04 
® ® 
♦ •* 
® 04 

TT 


i i 


r<»  ® 
os  rv 
® 04 
IO  © 

TT 


m o» 
r» « 
P»  04 
04  04 


tu  • 
Z © 
O < ■■=  2 


lie 

il5- 

£ 


04-*©tV0»«-P)©©©©04«*>© 
*»®04®©©0»I*S  f*»  «2  ® >4  00 

*-^04  04  04P»P5*0<Om©« 


©®©©©©©®©©©®©® 

6*«.©04®|V©04©P»©O6©P>© 

P>P»*-©®04©®PSP»<J=©ffi04 

“OegPSP<.^ffi®©a&PS®04©'>~ 

PS®®P)f^®<^|4®<^fo,v4>® 

r^^cyos®^^®®pjffi®®«y 

04©®®PS©09»-'>«"f>»©04©® 

»=^^^e40404P5lPJP1lPJ 


s TT 

i 

oS 

® to  <o 

© Pi  04 

1 TT 


© © 
® ® 

ff 


8 

£3 

© © 
04  ■*• 

ij 

• • 

m 

© © 

© 

1 1 

pj 

® 

o 

© 04 

® 

Pi  ® 

® © 

© O- 

1 

© 

® 

04  04 

| 

© 

1 1 

"O 

V» 

© 04 

1 

® 

€4  ® 

© © 

* 


04 

> 

04 


<"“®©®»-»-*-*-040404040404 

<P>  <Q^  i^>  ^p)  <Q  fQ  ^ 

i , f i.|  i,f  i.i  | !|  1 1 j q,j  i.i  I’f)  |(j  qtj  |,j 

«CMKJ®m®'4-'4-CS|CnvO'4-«-^- 

»-®OT©04P>»-©04t^P>0S»“0J 

®o>r'*o4©®in*-o>®o»>4°<s-"© 

*-©P>©»-©04®C4©OSPir'».® 

©P104©^-04««*-»=<s-040404 

Q Q ^ Q ^ Q Q @ ® ® ^ 

n«ll)44'414'4'4'4444'4 

Q Q ^p8  @ <p>  0 ^5) 

ujyuiyyyyuliijyyujuy 
^OO^-OtJC^^^-lOO^-GOi'Of^^ 
VOIOKJPS+NS^N!-*® 
Pi©®OS©<4°«-©©Pi©©©C4 
f*»*<04®P>©a>»-^f>*©04©© 
04«®*-*“*=  — 040404P3P>P}PS 

®®®o®s©sss®®®s 


Figure  3-7  continued 
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Figure  3-8  Resultant  Summary  Report 
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4.  SOFTWARE  MODULE  DESCRIPTIONS 


This  section  presents  a narrative  description  and  flow  diagram  for  each 
program  module  incorporated  in  the  waveform  signal  processing  software 
package. 

4.1  PROGRAM  "WGPROC" 

"WGPROC"  is  the  main  program  module  for  the  signal  processing  software 
package.  Its  main  purpose  is  to  interact  with  the  user  and,  based  upon  the 
response  of  the  user,  control  the  overall  flow  of  the  processing  procedure. 
The  flow  diagram  of  this  module  is  shown  in  Figure  4-1.  Not  shown  in  Figure 
4-1  is  the  option  to  select  full  or  partial  processing. 

4.2  SUBROUTINE  "LED" 

Subroutine  "LED"  contains  the  routine  which  is  used  to  detect  the  type  of 
time  synchronization  file.  If  the  time  synchronization  file  is  either  a 
standard  time  zero  file  or  a delayed  time  zero  file,  subroutine  "VWGEN1"  is 
called.  If  the  time  synchronization  file  is  inverted  (either  a full  or 
delayed  time  zero  file) , subroutine  "VWGEN2"  is  called.  The  flow  diagram  of 
this  module  is  shown  in  Figure  4-2. 

4.3  SUBROUTINES  " WGEN1 " AND  "WGEN2" 

Subroutines  "VWGEN1"  and  "VWGEN2"  are  used  to  detect  the  data  point  index 
number  corresponding  to  the  edges  in  the  time  synchronization  file.  ("VWGEN1" 
is  used  to  detect  the  leading  edges  and  "VWGEN2"  is  used  to  detect  the 
trailing  edges.)  As  indicated  in  Section  4.2,  "VWGEN1"  is  used  for  the 
standard  and  delayed  time  zero  files  and  "VWGEN2"  is  used  for  the  inverted 
standard  and  inverted  delayed  time  zero  files.  The  flow  diagrams  for 
"VWGEN1"  and  "WGEN2"  are  shown  in  Figures  4-3  and  4-4  respectively. 
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Figure  4-1  Flow  Diagram  of  "WGPROC" 
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Figure  4-1  continued 
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Figure  4-1  continued 
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Figure  4-1  continued 
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Figure  4-1  continued 
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Figure  4-2  Flow  Diagram  of  "LED" 
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Figure  4-3  Flow  Diagram  of  "VWGENl" 


48 


START 


XCHK  = 0.5  * (YMAX  - YMIN)  + YM IN 


Figure  4-4  Flow  Diagram  of  "VWGEN2" 
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4.4  SUBROUTINES  "CHAl"  AND  " CHA2 ” 


Subroutines  "CHAl"  and  "CHA2"  are  used  to  obtain  the  file  names  for  the 
data  to  be  processed  and  to  mark  the  starting  points  of  the  individual 
waveforms.  In  previous  versions  of  the  processing  software,  the  individual 
waveforms  were  physically  separated  into  individual  files.  In  this  version 
of  the  software,  the  data  is  stored  in  large  arrays  in  memory  with  the  index 
numbers  of  the  starting  points  marking  the  "separations."  "CHAl"  is  used  for 
the  standard  and  inverted  time  zero  files  and  "CHA2"  is  used  for  the  delayed 
standard  and  inverted  delayed  time  zero  files.  The  flow  diagrams  for  "CHAl" 
and  "CHA2"  are  shown  in  Figures  4-5  and  4-6  respectively.  The  E1(I)  used  in 
CHAl  and  CHA2  are  the  array  indices  of  the  data  points  in  the  time  reference 
waveform  that  first  exceed  one -half  the  distance  between  the  minimum  and 
maximum  values  of  that  waveform  for  a positive  (negative)  going  edge  in  a 
standard  (inverted)  time  reference  waveform. 

4.5  SUBROUTINES  "OPNFXL"  AND  "0PNFIL1" 

Subroutines  "OPNFXL"  and  "0PNFXL1"  are  used  to  create  the  names  for  the 
"pseudo -files"  which  are  displayed  in  order  to  allow  the  user  to  monitor  the 
progress  of  the  processing  procedure.  The  names  for  the  "pseudo -files"  are 
created  by  changing  the  first  character  of  the  extension  of  the  data  channel 
file  name  (entered  by  the  user  in  either  "CHAl"  or  "CHA2")  to  a character 
which  represents  the  waveform  being  processed.  (The  characters  used  are  "F" 
for  the  rectangle  waveform;  "C"  for  the  half-sine  waveform;  "E"  for  the 
staircase  waveform;  "V"  for  the  sum-of-sines  waveform;  and  "D"  for  the  crash 
waveform.)  Flow  diagrams  for  "OPNFXL"  and  "0PNFXL1"  are  shown  in  Figures  4-7 
and  4-8  respectively. 

4.6  SUBROUTINE  "CDNUM" 

Subroutine  "CDNUM"  is  used  to  extract  the  test  number  from  the  data 
channel  file  name  and  to  ensure  that  all  of  the  files  to  be  processed  are 
from  the  same  test.  If  data  files  are  specified  from  more  than  one  test,  a 
message  is  given  and  a flag  is  set  which  causes  the  program  to  cease  execu- 
tion. The  flow  diagram  of  this  module  is  shown  in  Figure  4-9.  In  this 
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Figure  4-5  Flow  Diagram  of  "CHAl" 
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Figure  4-5  continued 
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Figure  4-6  Flow  Diagram  of  "CHA2" 
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Figure  4-6  continued 
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Figure  4-7  Flow  Diagram  of  "OPNFIL" 
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Figure  4-8  Flow  Diagram  of  "0PNFXL1" 
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Figure  4-9  Flow  Diagram  of  "CDNUM" 
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subroutine,  CHANUMB  is  an  integer  that  is  incremented  by  one  each  time  an 
input  data  file  is  read  into  a storage  array  in  CHAl  or  CHA2 . 

4.7  SUBROUTINE  "RESULT" 

Subroutine  "RESULT"  generates  a data  array  which  is  the  resultant  of  X 
and  Z component  input  data  files.  This  module  also  creates  an  output  file 
called  "RESULT. SUM"  which  contains  a listing  of  the  pairings  of  data  files 
which  were  used  to  create  the  resultant  data  arrays.  As  with  the  files 
"created"  in  the  separation  of  the  individual  waveforms,  the  output  resultant 
files  do  not  physically  exist  but,  are  "pseudo-files."  The  flow  diagram  of 
this  module  is  shown  in  Figure  4-10. 

4.8  SUBROUTINE  "HIC3" 

Subroutine  "HIC3"  calculates  the  HIC  value  of  a resultant  data  array. 
This  module  is  a modified  version  of  the  HIC3  routine  written  for  NHTSA  by 
S . Mentzer  of  Automated  Sciences  Group . The  modifications  made  remove  the 
read/write  calls  t©  the  system  storage  device  and  the  creation  of  a report 
file  which  stated  the  value  of  the  HIC  calculated  for  the  data  file.  The 
flow  diagram  of  this  module  is  shown  in  Figure  4-11. 

4.9  SUBROUTINE  "PART" 

Subroutine  "PART"  performs  partitioning  of  possible  optimal  regions  found 
during  the  execution  of  the  HIC3  subroutine.  The  flow  diagram  of  this  module 
is  shown  in  Figure  4-12. 

4.10  SUBROUTINE  "EXTRACT" 

Subroutine  "EXTRACT"  extracts  the  file  name  out  of  the  data  file 
specification.  (The  data  file  specification  can  include  any  or  all  of  the 
following:  communications  node,  logical  unit  device  name,  account  name,  file 

name,  extension,  and  version  number.)  The  flow  diagram  of  this  module  is 
shown  in  Figure  4-13. 
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Figure  4-10  Flow  Diagram  of  "RESULT" 
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Figure  4-11  Flow  Diagram  of  "HIC3 
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Figure  4-13  Flow  Diagram  of  "EXTRACT" 
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4.11  SUBROUTINE  "HSINE" 


Subroutine  "HSINE"  performs  four  calculations/comparisons  on  the  half- 
sine waveform.  These  calculations  are: 

1.  Channel- to -channel  time  deviation.  This  comparison  is  carried  out 
between  the  two  data  channels  which  contain  the  X-axis  and  Z-axis 
half- sine  waveform  components. 

2.  X-axis  component  time  deviation  from  theoretical. 

3.  Z-axis  component  time  deviation  from  theoretical. 

4.  HIC  deviation.  This  comparison  is  carried  out  against  a calculated 
theoretical  HIC  value  which  is  calculated  in  subroutine  "THCALC." 

The  results  of  the  comparisons  made  in  this  module  are  written  to  a 
summary  file  (see  Table  3-1  for  the  file  naming  convention)  . The  flow 
diagram  of  this  module  is  shown  in  Figure  4-14. 

4.12  SUBROUTINE  " ERR2 " 

Subroutine  "ERR2"  is  used  to  interact  with  the  user  if  an  error  occurs 
while  a data  file  is  being  read.  The  flow  diagram  of  this  module  is  shown  in 
Figure  4-15. 

• 4.13  SUBROUTINE  "CRASH" 

Subroutine  "CRASH"  performs  a comparison  between  the  HIC  value  obtained 
from  the  resultant  of  the  recorded  X-axis  and  Z-axis  crash  data  and  a value 
of  930  which  was  calculated  from  the  original  crash  test  data  using  HIC3.  The 
results  of  this  comparison  are  then  written  to  a summary  file  and  a statisti- 
cal data  file.  (See  Table  3-1  for  the  file  naming  convention  and  the 
content  of  the  statistical  data  file.)  The  flow  diagram  of  this  module  is 
shown  in  Figure  4-16. 
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Figure  4-16  Flow  Diagram  of  "CRASH" 
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4.14  SUBROUTINE  "FRSP3" 


Subroutine  "FRSP3"  computes  the  amplitude  frequency  response  of  a data 
acquisition  channel  in  the  following  sequence: 

1.  The  data  file  of  the  recorded  data  acquisition  channel  response  to 
the  sum- of- sines  waveform  is  read  into  an  array. 

2.  The  index  of  the  response  array  is  displaced  by  1 from  "0  to  N+l"  to 

"1  to  N+2"  for  compatibility  with  the  fast  Fourier  transform  (FFT) 

processing. 

3.  The  input  waveform  data  is  recreated,  sampled  at  the  same  frequency 
as  was  the  response  data. 

4.  Subroutine  "SGAUSS"  is  called  to  generate  an  N point  Gaussian  data 
window . 

5.  The  waveform  and  response  arrays  are  multiplied  by  the  Gaussian 
window  weighting  function. 

6.  The  discrete  Fourier  transforms  (DFT)  of  the  weighted  waveform  and 
response  arrays  are  computed  using  the  FAST*  fast  Fourier  transform 
subroutine . 

7.  The  index  of  each  DFT  output  array  element  whose  corresponding 
frequency  is  closest  to  each  signal  frequency  is  determined.  ' 

8.  The  magnitude  of  each  DFT  output  array  element  that  corresponds  to  a 
signal  frequency  is  determined. 

9.  The  scaling  factor  is  computed. 

*The  FAST  FFT  algorithm  and  implementing  software  were  developed  by  G.D. 
Bergland  and  M.T.  Dolan  of  Bell  Laboratories,  Murray  Hill,  N.J.  and  published 
in  Programs  for  Digital  Signal  Processing  by  the  IEEE  Press,  New  York,  © 
1979  IEEE. 
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10. 


The  scaled  ratio  of  the  DFT  of  the  response  to  the  DFT  of  the 
waveform  at  each  signal  frequency  is  computed.  This  is  the 
numeric  amplitude  ratio  at  each  signal  frequency. 

11.  The  amplitude  ratio  (in  dB)  is  computed  at  each  signal 
frequency . 

12.  The  frequency  response  data  files  are  written. 


The  flow  diagram  of  this  module  is  shown  in  Figure  4-17. 
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Figure  4-17  Flow  Diagram  of  "FRSP3" 
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4.15  SUBROUTINE  "SGAUSS" 


Subroutine  "SGAUSS"  generates  an  N point  Gaussian  data  window. 

4.16  SUBROUTINE  "FAST" 

Subroutine  "FAST"  is  the  fast  Fourier  transform  routine  called  by 
"FRSP3 . " 

4.17  SUBROUTINE  "FR2TR" 

Subroutine  "FR2TR"  is  a radix  2 iteration  subroutine. 

4.18  SUBROUTINE  "FR4TR" 

Subroutine  "FR4TR"  is  a radix  4 iteration  subroutine. 

4.19  SUBROUTINE  " F0RD1" 

Subroutine  "F0RD1"  is  the  first  of  two  in-place  reordering  subroutines 
called  by  "FAST." 

4.20  SUBROUTINE  "FORD 2" 

Subroutine  "F0RD2"  is  the  second  of  two  in-place  reordering  subroutines 
called  by  "FAST." 

4.21  SUBROUTINE  "RECTAN" 

Subroutine  "RECTAN"  is  used  to  process  the  rectangular  waveform. 
"RECTAN"  separates  the  rectangular  waveform  into  ten  cycles.  Each  of  the 
cycles  is  then  inspected  with  respect  to  the  steady- state  amplitude  and 
amplitude  overshoot  on  both  the  positive  and  negative  sections  of  the  cycle 
as  well  as  for  time  .shift  within  the  channel.  The  flow  diagram  of  this 
module  is  shown  in  Figure  4-18. 


69 


Figure  4-18  Flow  Diagram  of  "RECTAN” 
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4.22  SUBROUTINE  "ZERO 2" 


Subroutine  "ZER02"  calculates  the  maximum  overshoot  and  the  steady-state 
amplitude  values  for  any  given  cycle  of  the  rectangular  waveform.  The  flow 
diagram  of  this  module  is  shown  in  Figure  4-19. 

4.23  SUBROUTINE  "COMPRE" 

Subroutine  "COMPRE"  performs  comparisons  between  the  results  of  the 
steady- state  amplitude  and  overshoot  and  time  shift  calculations  and  pre- 
defined tolerances.  There  are  four  time  comparisons  and  two  amplitude 

comparisons.  The  time  comparisons  are: 

1.  Deviation  between  actual  times  and  corresponding  theoretical  times 
and  between  the  deviations  and  the  tolerance. 

2.  Time  linearity.  (The  slope  of  a least  squares  fit  straight  line 
through  data  pairs  consisting  of  theoretical  time  as  the  X coordinate 
and  calculated  time  shift  as  the  Y coordinate . ) 

3.  Time  offset.  (The  Y-axis  intercept  of  a least  squares  fit  straight 
line  through  data  pairs  consisting  of  theoretical  time  as  the  X 
coordinate  and  calculated  time  shift  as  the  Y coordinate . ) 

4.  Channel- to -channel  time  deviation. 

The  amplitude  comparisons  are: 

1.  Steady-state  amplitude  deviation  from  theoretical  amplitude. 

2.  Amplitude  overshoot  relative  to  the  calculated  steady- state  ampli- 
tude . 

The  flow  diagram  of  this  module  is  shown  in  Figure  4-20. 


Figure  4-19  Flow  Diagram  of  "ZER02" 


Figure  4-20  Flow  Diagram  of  "COMPRE" 


73 


4.24  SUBROUTINE  "STAIR” 

Subroutine  "STAIR"  separates  the  staircase  waveform  into  12  periods 
(sections)  by  detecting  the  12  leading  edges  of  the  waveform.  Each  of  the  12 
periods  is  then  inspected  with  regards  to  the  steady- state  amplitude  and 
overshoot  and  the  time  shift  within  the  channel.  The  flow  diagram  of  this 
module  is  shown  in  Figure  4-21. 

4.25  SUBROUTINE  "ZER03" 

Subroutine  "ZER03"  performs  the  same  functions  for  the  staircase 
waveform  processing  as  subroutine  "ZER02"  does  for  the  rectangular  waveform 
processing.  For  this  reason,  no  discussion  or  flow  diagram  will  be  presented 
for  subroutine  "ZER03". 

4.26  SUBROUTINE  "COMSTR" 

Subroutine  "COMSTR"  performs  comparisons  between  the  results  of  the 
steady- state  amplitude  and  overshoot  and  time  shift  calculations  and  pre- 
defined tolerances.  There  are  three  time  comparisons  and  four  amplitude 
comparisons.  The  time  comparisons  are: 

1.  Deviation  between  the  recorded  waveform  and  the  theoretical  time. 

2.  Time  linearity.  (The  slope  of  a least  squares  fit  straight  line 
through  data  pairs  consisting  of  theoretical  time  as  the  X coordinate 
and  calculated  time  shift  as  the  Y coordinate.) 

3.  Time  offset.  (The  Y-axis  intercept  of  a least  squares  fit  straight 
line  through  data  pairs  consisting  of  theoretical  time  as  the  X 
coordinate  and  calculated  time  shift  as  the  Y coordinate . ) 
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Figure  4-21  Flow  Diagram  of  "STAIR" 
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The  amplitude  comparisons  are: 


1.  Amplitude  linearity.  (The  slope  of  a least  squares  fit  straight 
line  through  data  pairs  consisting  of  theoretical  amplitude  as  the  X 
coordinate  and  measured  steady- state  amplitude  deviation  as  the  Y 
coordinate . ) 

2.  Amplitude  offset.  (The  Y-axis  intercept  of  a least  squares  fit 
straight  line  through  data  pairs  consisting  of  theoretical  amplitude 
as  the  X coordinate  and  measured  steady- state  amplitude  deviation  as 
the  Y coordinate . ) 

3.  Steady-state  amplitude  deviation  from  theoretical  amplitude. 

4.  Amplitude  overshoot  relative  to  the  calculated  steady- state  ampli- 
tude . 

In  addition  to  the  comparisons , this  subroutine  also  generates  a summary 
file  and  five  statistical  data  files  which  contain  the  results  of  the 
comparisons.  (See  Table  2-1  for  the  file  naming  convention  for  these 
files.)  The  flow  diagram  of  this  module  is  shown  in  Figure  4-22. 

4.27  SUBROUTINE  " C0MPR2 " 

Subroutine  "G0MPR2"  performs  a comparison  of  channel- to -channel  time 
deviation  for  each  period  of  the  staircase  waveform.  The  results  of  this 
comparison  is  added  to  the  summary  file  originally  created  in  subroutine 
"COMSTR".  The  flow  diagram  of  this  module  is  shown  in  Figure  4-23. 

4.28  SUBROUTINE  "UDSIO" 

Subroutine  "UDSIO"  is  a modification  of  the  UDS  (User  Data  Set) 
read/write  utility  routine  developed  for  NHTSA.  The  modifications  made 
relate  to  the  error  handling  routines  found  within  the  subroutine.  The  flow 
diagram  of  this  module  is  shown  in  Figure  4-24. 
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Figure  4-22  Flow  Diagram  of  "COMSTR" 
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Figure  4-23  Flow  Diagram  of  "C0MPR2" 
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Figure  4-24  Flow  Diagram  of  "UDSIO" 
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4.29  SUBROUTINE  "GETWAVE" 


Subroutine  "GETWAVE”  is  a routine  which  allows  the  reduction  of  system 
input/output.  Instead  of  writing  to  and  reading  from  the  system  storage 
device,  the  data  is  placed  into  large  arrays  thus  forming  a crude  "virtual 
disk."  All  necessary  functions  pertinent  to  the  replacement  of  the  read/write 
process  are  handled  within  this  subroutine.  The  flow  diagram  of  this  module 
is  shown  in  Figure  4-25. 

4.30  SUBROUTINE  "GETIND" 

Subroutine  "GETIND"  determines  the  index  of  the  point  in  the  data  array 
which  is  closest  to  the  calculated  time  of  occurrence  of  the  maximum  data 
value  in  the  half-sine  waveform  components.  The  approach  taken  is  to  average 
the  times  the  data  passes  through  a certain  amplitude  level  which  is  ex- 
pressed as  a percentage  of  the  theoretical  maximum  (+200  g's  depending  upon 
the  component  being  considered) . Interpolation  is  used  to  improve  the  time 
calculation.  The  flow  diagram  of  this  module  is  shown  in  Figure  4-26. 

4.31  SUBROUTINE  "THCALC" 

Subroutine  "THCALC"  obtains,  through  multiple  steps,  the  value  of  the 
theoretical  HIC  if  the  curve  were  to  be  made  up  of  two  components  (X  and  Z) 
which  have  the  same  amplitudes  and  duration  as  the  actual  (recorded)  data. 
The  steps  taken  to  obtain  the  theoretical  HIC  value  are: 

1.  Determine  the  equations  of  the  two  half-sine  curves  which  pass 
through  the  points  determined  in  subroutine  "GETIND" . 

2 . Generate  arrays  of  data  which  are  the  equivalent  of  the  equations  by 
sampling  the  equations  at  the  same  sample  rate  as  was  used  to  record 
the  waveform  data. 

3.  Obtain  the  resultant  of  the  two  arrays. 
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NO 


YES 


▼ 

RETRIEVE  VALUES 
OF  INDICIES 
DATA  CORRESPONDING 
TO  INDIVIDUAL 
WAVEFORM 
SPECIFIED 


YES 


▼ 

STORE  VALUES  OF 
INDICIES  MARKING 
START  AND  END  OF 
INDIVIDUAL 
WAVEFORM 
SPECIFIED 


Figure  4-25  Flow  Diagram  of  "GETWAVE" 
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FIND  TIMES  AT  WHICH  THE 
AMPLITUDE  CROSSES  THE  SPECIFIED  LEVEL 


CALCULATE  TIME  OF  MAXIMUM  AMPLITUDE 
(AVERAGE  OF  CROSSING  TIMES) 


CALCULATE  TIME  DURATION  THE  AMPLITUDE  IS 
"GREATER"  OR  EQUAL  TO  THE  SPECIFIED  LEVEL 


DETERMINE  INDEX  OF  DATA  POINT  CLOSEST  TO 
THE  CALCULATED  TIME  OF  MAXIMUM  AMPLITUDE 


Figure  4-26  Flow  Diagram  of  "GETIND" 
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4.  Calculate  HIC  and  the  upper  and  lower  time  boundaries  (T1  and  T2) 
using  a modified  version  of  the  HIC3  subroutine. 

The  flow  diagram  of  this  module  is  shown  in  Figure  4-27. 


Figure  4-27  Flow  Diagram  of  "THCALC" 
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5 . PROGRAM  LISTINGS 


This  section  presents  the  source  code  listings  for  the  waveform 
generator  signal  processing  software  package.  The  package  is  divided  into 
five  segments  with  each  segment  containing  one  or  more  program  modules.  The 
program  listings  will  be  presented  in  the  order  as  they  appear  in  the 
segments . For  those  readers  interested  in  locating  the  source  code  on  the 
NHTSA  VAX  11/780,  the  segments  are  located  in  user  area  [JOHNN.WGPROG]  and 
have  the  following  file  names:  WGPROG.FQR,  FILMGT.FOR,  WGSUBS.FOR,  REG. FOR, 
and  STR.FOR. 

The  first  segment,  "WGPROC",  contains  only  the  main  program  module.  The 
source  code  listing  is: 
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Q •k'k-k-k-k*-k'k-k'k-k-k'k-k-k-k-k&'k-k 

C **  MAIN  PROGRAM  ** 

Q ******■*■****•**■**•*■**** 

PROGRAM  WGPROC 
C 

C VARIABLE  DECLARATION 
C 

DIMENSION  Y( -1000: 10000) 

INCLUDE  'DISK$OVR: [ASGPROG] COMVAR. LIS ' 

CHARACTER  ITIMEZ  *30,DUM7  *1 
DIMENSION  IEDGE(IO) 

INTEGER  OPTION , CHANUMB , CNTFL , EFLAG ( 9 ,18) , DFLAG , DELFLAG 
REAL  T(18 , 23) 

CHARACTER  *4  TSTCDNUM 

CHARACTER  *30  CHAFIL(108) ,RESFIL(72) 

CHARACTER  *9  DAY .HOUR 
Character*l  ANSWER 
C 

C COMMON  BLOCKS 
C 

COMMON  /CBFIL/  CHAFIL 

* /CBCH/  CHANUMB 

* /CBRES/  RESFIL 

* /CBCDNUM/  TSTCDNUM 
C 

C Initialization 
C 

DO  1 1-1,108 
- CHAFIL(I)— ' 

1 CONTINUE 

DO  2 1-1,72 
RESFIL(I)— ' 

2 CONTINUE 
C 

WRITE (6 , 10) 

10  FORMAT (//IX, ’ ***  WAVEFORM  GENERATOR  PROCESSING  SOFTWARE  ***'//) 
CALL  DATE (DAY) 

CALL  TIME (HOUR) 

WRITE (6,12) DAY , HOUR 
12  FORMAT (10X,A9,10X,A9) 

WRITE  (6,101) 

101  FORMAT  (//3x,  'DO  YOU  WANT  ANALYSIS  OF  HALF-SINE  AND  CRASH  PULSES 

* IN  THE  PROCESSING?  (Y  or  N) ' , /lx, '»',$) 

READ  (5,102)  ANSWER 

102  FORMAT  (A) 

1000  WRITE(6 , 100) 

100  FORMAT (//IX, 'OPTION  1 : FINDING  DATA  POINT  INDEX  NUMBERS  ', 

& 'CORRESPONDING  TOEDGES' , /IX, 'OPTION  2 : ENTER  CORRECT  DATA  » , 

6c  'POINT  INDEX  NUMBERS ', /IX, ' ENTER  DESIRED  OPTION  NUMBER  (1/2)', 

& /IX,  '»>  ',$) 

READ (5, 200)  OPTION 
200  FORMAT(Il) 
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IF  ( OPTION. NE. 1. AND. OPTION. NE. 2)  THEN 
WRITE (6 , 3) 

3 FORMAT (/IX, ' ***  ERROR  - INCORRECT  OPTION  NUMBER  ***  ' , 

& /IX,’ ' ***  TRY  AGAIN  ***') 

GOTO  1000 
END  IF 

IF  (OPTION. EQ.l)  THEN 
C 

C OPTION  1 
C 

CALL  LED ( IFLAG2 , IEDGE , XTFL) 

IF  (IFLAG2  .EQ.  0 ) GOTO  2000 
IF  (XFLAG2  „EQ.  10)  GOTO  9999 

C 

C plot  of  time  zero  is  necessary 

C 

WRITE (6, 4) 

4 FORMAT (/IX, ' WARNINGS ' , 

& /IX, ' ADDITIONAL  EDGES  IN  TIME  ZERO  FILE  HAVE  BEEN' , 

& ' DETECTED' , /IX, ! OBTAIN  A PLOT  OF  TIME  ZERO  FILE  ' 

& /IX,’  AND  USE  OPTION  2 TO  ENTER  THE  CORRECT8, 

& ' DATA  POINT  NUMBERS8) 

GOTO  9999 
C 

C correct  number  of  edges  detected 
C 

2000  IF  (ITFL.EQ.1.0R.ITFL.EQ.3)  GOTO  2010 
C 

C delayed  TIMEZ  / inverted  delayed  TIMEZ 
C 

CALL  CLOSE (1) 

GALL  CHA2 (IEDGE, IFLAG2) 

IF  (IFLAG2  .EQ.  10  ) GOTO  9999 
CALL  CDNUM(IFLAG2) 

IF  (I FLAG 2 .EQ.  10)  GOTO  9999 
DELFLAG  -1 
GOTO  4000 
C 

C TIMEZ  / inverted  TIMEZ 
C 

2010  CALL  CLOSE (1) 

CALL  CHA1 ( IEDGE, IFLAG2) 

IF  (I FLAG 2 .EQ.  10)  GOTO  9999 
CALL  GDNUM(IFLAG2) 

IF  (I FLAG 2 .EQ.  10)  GOTO  9999 
DELFLAG  - 0 
GOTO  4000 
END  IF 
C 

C OPTION  2 
C 

IF  (OPTION. EQ. 2)  THEN 
5000  WRITE(5 , 9000) 

9000  FORMAT (/IX, 'ENTER  1 IF  TIMING  FILE  IS  A STANDARD  TIME  ZERO', 
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& ' FILE’, /IX,'  2 IF  TIMING  FILE  IS  A DELAYED  TIME  ZERO', 

& ' FILE', /IX,'  3 IF  TIMING  FILE  IS  AN  INVERTED  TIME  ZERO’, 

& ' FILE', /IX,'  4 IF  TIMING  FILE  IS  A DELAYED  & INVERTED', 

& ' TIME  ZERO  FILE ' ) 

READ(5 , 9001)  ID1 
9001  FORMAT (12) 

IF  (ID1.EQ.2.0R.ID1.EQ.4)  ICOUNT  - 2 
IF  (ID1.EQ.1.0R.ID1.EQ.3)  ICOUNT  - 4 
DO  5010  I - 1, ICOUNT 
WRITE(6 , 130)  I 

130  FORMAT (/IX, 'ENTER  DATA  POINT  INDEX  NUMBER  ',12, 

& /IX , ' »>  ',$) 

READ(5 , 210)  IEDGE(I) 

210  FORMAT (14) 

5010  CONTINUE 

IF  (ID1.EQ.1.0R.ID1.EQ.3)  GOTO  5050 
C 

C delayed  / inverted  delayed  TIMEZ 
C 

CALL  CLOSE(l) 

CALL  CHA2 ( IEDGE , IFLAG2 ) 

IF  (IFLAG2  .EQ.  10)  GOTO  9999 
CALL  CDNUM(IFLAG2) 

IF  (IFLAG2  .EQ.  10)  GOTO  9999 
DELFLAG  - 1 
GOTO  4000 
C 

C TIMEZ-/  inverted  TIMEZ 
C 

5050  CALL  CLOSE(l) 

CALL  CHA1 ( IEDGE, IFLAG2) 

IF  (IFLAG2  .EQ.  10)  GOTO  9999 
CALL  CDNUM(IFLAG2) 

IF  (IFLAG2  .EQ.  10)  GOTO  9999 
DELFLAG  - 0 
END  IF 

4000  WRITE(6 , 135) 

135  FORMAT (/IX, ' PROCESSING  SEQUENCE: • , 

& /10X, 'RECTANGLE  WAVEFORM' ,/10X, 'STAIR  WAVEFORM' , 

& /10X, 'HALF- SINE  WAVEFORM’ ,/10X, 'CRASH  WAVEFORM' , 

& /10X, 'SUM -OF -SINE  WAVEFORM') 

CNTFL  - 1 
C 

C AUTOMATIC  PROCESSING 
C 

4010  DFLAG  - 0 

IF  ( DELFLAG. EQ.l)  GOTO  4500 
C 

C RECTANGLE  PROCESSOR 
C 

CALL  RECTAN(ICHAN, DFLAG, CNTFL) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 
CALL  COMPRE(ICHAN, DFLAG) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 
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oci  ooo  ooo  oo  ooo 


STAIR  PROCESSOR 
DFLAG  - 0 

CALL  STAIR(ICHAN, DFLAG, CNTFL) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 
CALL  COMSTR ( T , I CHAN , EFLAG , DFLAG ) 
IF  (DFLAG  .EQ.  1)  GOTO  9990 
CALL  C0MPR2(T,ICHAN, EFLAG) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 


IF  (ANSWER  .EQ. 
IF  (ANSWER  .EQ. 
GO  TO  4600 
4601  CONTINUE 

HALF- SINE  PROCESSOR 


'Y')  GO  TO  4601 
9y' ) GO  TO  4601 


DFLAG  - 0 
WRITE (6 , 6) 

6 FORMAT (/IX, * **  HALF-SINE  WAVEFORM  PROCESSOR  **  '/) 
IKJNUM  - -3 

IWF-2 

XSLT-1 

CALL  RESULT  ( DFLAG , IKJNUM , CNTFL , IWF , DELFLAG ) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 
CALL  HIC3 (DFLAG , CNTFL , ISLT) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 
CALL  HSINE (DFLAG, CNTFL) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 

CRASH  WAVEFORM  PROCESSOR 

4500  DFLAG  - 0 
WRITE (6 , 7) 

7 FORMAT (/IX, ' **  CRASH  WAVEFORM  PROCESSOR  **  '/) 
IKJNUM  » 0 

IWF“5 

ISLT-2 

CALL  RESULT (DFLAG, IKJNUM, CNTFL, IWF, DELFLAG) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 
CALL  HIC3 (DFLAG, CNTFL, ISLT) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 
CALL  CRASH (DFLAG, CNTFL) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 

SUM- OF -SINE  WAVEFORM  PROCESSOR 


4600  WRITE (6 , 8) 

8 FORMAT (/IX, ' **  SUM-OF-SINE  WAVEFORM  PROCESSOR  **'/) 

DFLAG  - 0 

CALL  FRSP3 (DFLAG, CNTFL) 

IF  (DFLAG  .EQ.  1)  GOTO  9990 
GOTO  9999 


ERROR 
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c 

9990  WRITE (6 , 9) 

9 FORMAT (//IX, '*****  CRITICAL  ERROR  *****'/) 

9999  CONTINUE 

WRITE(6 , 11) 

11  FORMAT (/IX/) 

STOP  'END  OF  PROCESSING  RUN' 

END 
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The  second  segment,  "FXLMGT"  contains  the  following  modules: 


LED 

VWGEN1 

VWGEN2 

CHA1 

OPNFIL 

GHA2 

0PNFIL1 

GDNUM 

UDSIO 

GETWAVE 


The  source  code  listing  for  these  modules  is: 
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Q •kk-k-k'krkk-k-kirk-k'kirk'kir/ck'kick-k-kirkirk-k-JrJrk-k-kk-kk'/rkk'k'k'k-k-k-k-k'k-k-k 

C * SUBROUTINE  LED  * 

C ■kk-k-kk-kirk-krk-k-k-k-k-k-k-k-k-k’k-k-kick-k-k'kk-k-k'k-k-k-k-k'k-k-k'k-k-k-kk-k-k-k-k-k-k-k 

SUBROUTINE  LED ( IFLAG2 , IED1 , ITFL) 

CHARACTER  IFLNM1*60 , DUM  *1 
DIMENSION  IED1(10) 

INTEGER  BLNK , IUNITS ( 2 ) 

INCLUDE  ' DISK$OVR: [ ASGPROG ] COMVAR .LIS' 

REAL  Y( -1000: 10000) , YMA, YMI 


DATA  BLNK/'  '/..NO/ 'NO'/ 

10  IFLAG2  - 0 
WRITE (6 , 20) 

20  FORMAT (/3X, 'ENTER  TIME  SYNCHRONIZATION  FILE  NAME' ,/lX, '»' ,$) 

READ (5,25) IFLNM1 
FORMAT (A60) 

WRITE(6 , *) ' ' 

CALL  UDSIO('R' , 'F' , IFLNM1 , 1 , Y, IERR) 

IF  (IERR  .EQ.  0)  GOTO  30 

ERROR 

WRITE(6 , *)  ' ERROR!!!' 

WRITE (6,*)  • ***  FILE  CANNOT  BE  ACCESSED  ***  ’ 

GOTO  1000 

NOPTS  - (NLP  - NFP)  + 1 

FINDING  MINIMUM  AND  MAXIMUM 

YMA  - -1.E10 
YMI  - 1.E10 
DO  3000  I - NFP, NLP 
YMA  - AMAX1 (YMA , Y ( I ) ) 

YMI  - AMIN1(YMI ,Y(I) ) 

3000  CONTINUE 
ID1  - 1 
ITFL  * ID1 

IF  (ID1  .EQ.  1)  GOTO  100 
IF  (ID1  .EQ.  2)  GOTO  150 
IF  (ID1  .EQ.  3)  GOTO  200 
IF  (ID1  .EQ.  4)  GOTO  250 

ERROR  IN  HEADING  INFORMATION  FOR  TIME  ZERO  CHANNEL  OPTION  USED 
WRITE (6 , *)  ’ ERROR!!!’ 

WRITE (6,*)  * ***  HEADING  INFORMATION  DOES  NOT  HAVE  CORRECT’ , 

* ’ TIME  ZERO  OPTION  IDENTIFIER  ’ 

GOTO  1000 


100  WRITE ( 6 , *)  ’ **NOTE  : TIME  ZERO  HAS  BEEN  RECORDED**’ 

GOTO  300 
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150  WRITE (6,*)  ' **NOTE  : DELAYED  TIME  ZERO  HAS  BEEN  RECORDED**' 
GOTO  300 

200  WRITE (6 , *)  ' **NOTE  : INVERTED  TIME  ZERO  HAS  BEEN  RECORDED**’ 

GOTO  400 

250  WRITE(6 ,*)  ' **NOTE  : INVERTED  DELAYED  TIME  ZERO  HAS  BEEN 

* 'RECORDED**' 

GOTO  400 
C 

C TIME  ZERO  AND  DELAYED  TIME  ZERO 
C 

300  GLOSE(l) 

CALL  UDSIO('R'  , 'F'  , XFLNM1 , 1 , Y,  I ERR) 

IF  (IERR  .NE.  0)  GOTO  27 

CALL  VWGEN1 (NOPTS , Y , IFLAG1 , IED1 , YMA , YMI , NFP , NLP) 

IF  (XFLAG1  .GT.  4)  IFLAG2  - 1 

IF  (I FLAG 1 .GT.  2 .AND.  ID1  .EQ.  2)  I FLAG 2 - 1 
GOTO  1100 
C 

C INVERTED  TIME  ZERO  AND  INVERTED  DELAYED  TIME  ZERO 
C 

400  CLOSE (1) 

CALL  UDSIOCR'  , 'F*  , IFLNM1 1 1 , Y,  IERR) 

IF  (IERR  .NE.  0)  GOTO  27 

CALL  VWGEN2 (NOPTS , Y , I FLAG 1 , IED1 , YMA , YMI , NFP , NLP) 

IF  (IFLAG1  .GT.  4)  I FLAG 2 - 1 

IF  (IFLAG1  .GT.  2 .AND.  ID1  .EQ.  4 ) IFLAG2  * 1 
GOTO  1100 
C 
C 
C 

1000  CLOSE (1) 

1002  WRITE (6 , 1005) 

1005  FORMAT (3X, 'ENTER:  [Y]  TO  TRY  AGAIN' ,/10X, 

* ' [N]  TO  TERMINATE  THE  PROCESSES'  9/lX,  '»'  ,$) 

READ (5  , 1007)  DUM 

1007  FORMAT (Al) 

IF  (DUM  .EQ.  'Y')  GOTO  10 
IF  (DUM  .EQ.  'N')  GOTO  1105 
GOTO  1002 

1105  IFLAG2  - 10 
1100  CONTINUE 
RETURN 
END 
C 

C * SUBROUTINE  VWGEN1  * 

C ****^^*^**ilHt'****-jlr*************************************'**-***** 

SUBROUTINE  VWGEN1 (NOPTS , X , IFLAG1 , IED2 , YMAX , YMIN , NFP , NLP) 

CHARACTER* 3 C3 

DIMENSION  X( -1000: 10000) 

DIMENSION  IED2 (10) 

REAL  YMAX, YMIN 
C 

XCHK  =0.5  *( YMAX- YMIN) +YMIN 
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I FLAG 1 - 0 
J - 1 
ISW-0 

DO  100  I - NFP.NLP 
XCHK2  - X(I) 

IF  (ISW  .EQ.  1 .AND.  XCHK2  .LT.  XCHK)  ISW  - 0 

IF  (ISW  .EQ.  1 .AND.  XCHK2  .GT.  XCHK)  GO  TO  100 

IF  (ISW  .EQ.  0 .AND.  XCHK2  .LT.  XCHK)  GO  TO  100 

ISW  - 1 

K3  - I 

IED2 (J)  - K3 

J - J + 1 

RK3-ABS ( FLOAT (K3 ) /10 . 0 ) 

RK3-RK3 - ABS ( FLOAT ( INT (RK3 ) ) ) 

KKK3-JNINT (RK3*10 . 0 ) 

G3— ' th’ 

IF(KKK3 . EQ. 1)  C3-'  st' 

IF(KKK3 . EQ. 2)  C3-'  nd' 

IF(KKK3 . EQ. 3)  C3-*  rd’ 

WRITE (6 , 1)  K3.C3 

1 . FORMAT (3X, 'LEADING  EDGE  AT* ,16, A3, ' DATA  POINT') 

IFLAG1  - IFLAG1  + 1 
100  CONTINUE 
CLOSE (1) 

RETURN 

END 

C 

C * SUBROUTINE  VWGEN2  * 


C 


SUBROUTINE  VWGEN2 (NOPTS , X , I FLAG 1 , IED2 , YMAX , YMIN , NFP , NLP ) 

CHARACTERS  C3 

DIMENSION  X( *1000: 10000) 

DIMENSION  IED2(10) 

REAL  YMIN, YMAX 

XCHK  - 0.5* (YMAX  - YMIN)  + YMIN 
ISW  - 0 
TFLAG1  - 0 
J - 1 

DO  100  I - NFP.NLP 
XCHK2  - X(I) 


IF 

(ISW 

.EQ. 

1 

.AND. 

XCHK2 

.GT. 

XCHK) 

ISW  - 0 

IF 

(ISW 

.EQ. 

1 

.AND. 

XCHK2 

.LT. 

XCHK) 

GOTO 

100 

IF 

(ISW 

.EQ. 

0 

.AND. 

XCHK2 

.GT. 

XCHK) 

GOTO 

100 

ISW  - 1 
K3  - I 

IED2(J)  - K3 
J - J+  1 

RK3-ABS ( FLOAT (K3 ) /10 . 0 ) 
RK3-RK3 -ABS ( FLOAT ( INT (RK3 )) ) 
KKK3-JNINT (RK3*10 . 0 ) 

C3— ' th' 

IF(KKK3 . EQ. 1)  C3-'  st' 
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IF(KKK3 . EQ. 2)  C3-'  nd' 

IF(KKK3 . EQ. 3)  C3-'  rd’ 

WRITE (6 , 1)  K3,C3 

1 FORMAT (3X, 'TRAILING  EDGE  AT', 16, A3,'  DATA  POINT') 

IFLAG1  - IFLAG1  + 1 
100  CONTINUE 
CLOSE (1) 

RETURN 

END 


C 

c 

c *■ 

G * 

C * 

c * 
c *■ 

SUBROUTINE  CHAl (El , IFLAG2 ) 

DIMENSION  ¥(-1000:10000) 

INTEGER  NPT ( 5 ) , BOUND ( 6 ) , El ( 10 ) , CHANUMB 
INCLUDE  ' DISK$QVR: [ASGPROG]COMVAR„LIS' 
CHARACTER  IFL*60 , DUMMY* 1 
CHARACTER  *30  CHAFIL(108) 

CHARACTER  JNFIL*30 , BLANK , I BUM 


COMMON 

/GBY/ 

Y 

* 

/GBE  / 

E, NPT, BOUND 

* 

/CBNPT/ 

NPTS 

* 

/CBFIL/  CHAFIL 

* 

/CBCH  /CHANUMB 

C 

I BLAG 2 - 0 
CHANUMB  - 0 
KCH  - 0 
INDX-0 


SUBROUTINE  CHANNELl 
TIMEZ , INVERTED  TIMEZ 


■k 

■k 

■k 


BLANK**'  ' 


GLOSE(l) 

WRITE(6 , 700) 

700  FORMAT (/IX, '** 

WRITE(6 , 701) 

701  FORMAT (//IX, 'NOTE: 

& /IX,' 

& /IX,' 

& /IX,' 

& /IX,' 

10  WRITE(6 , 710) 

710  FORMAT (//IX, 'ENTER  DATA  FILE  NAME  (XXXXXXXXX.XXX)  ' , 

* /IX, ’IF  NO  MORE  FILES,  TYPE  "CNTR  Z"  ’ ,/lX,  '»'  , $) 

READ(5 , 520  , END  - 4000)  IFL 
520  FORMAT (A60) 


BEGINNING  OF  INDIVIDUAL  WAVEFORM  SEPARATION  **') 


THE  TOTAL  NUMBER  OF  FILE  NAMES  YOU  ENTER' , 
MUST  BE  A MULTIPLE  OF  TWO.  THE  FIRST,  THIRD,' , 
FIFTH,  etc . FILE  MUST  BE  AN  X- COMPONENT  FILE' , 
AND  THE  SECOND,  FOURTH,  SIXTH,  etc.,  FILE', 
MUST  BE  A Z = COMPONENT  FILE.'//) 


C 

C READING  HEADING  FROM  THE  FILE  CONTAINING  5 WAVEFORMS 
C 


CALL  UDSIO( ' R' , 'F' , IFL, 1 , Y, IERR) 
IF  (IERR  .NE.  0 ) THEN 
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o r>  o o o o o noo 


ERROR  IN  INPUT  FILE 

20  WRITE (6 , *)  • ERROR!!!' 

WRITE (6,*)  ’ ***  FILE  CANNOT  BE  ACCESSED  ***' 

GOTO  8000 
END  IF 

CALL  EXTRACT (IFL) 

CORRECT  UDSIO  FILE 
BOUND  GENERATION 

30  INDX-INDX+1 

CALL  GETWAVE ( INDX ,0,0, NFP , NLP , DEL , ' I ' ,Y) 

CHANUMB  - CHANUMB  + 1 
KCH  - KCH  + 1 
CHAFIL(KCH)  - IFL(1:30) 

BOUND ( 1 ) —El ( 1 ) 

IF(E1(1) .LT.l)  BOUND (1)—1 
B0UND(3)-E1(2) 

B0UND(4)-E1(3) 

BOUND (5)— El (4) 

NOPTS  - (NLP  - NFP)  + 1 
BOUND  ( 6 ) —NLP 
NPTS-E1(2)-E1(1) 

BOUND ( 2 ) -100 . 0/125 . 0*NPTS+E1 ( 1 ) + 1 
DEL-0. 375/(El(4)-El(l)) 

WRITE (6 ,*) ' DEL-', DEL 
DO  40  I -1,5 

NPT ( I ) — BOUND (1+1)  - BOUND ( I ) 

40  CONTINUE 
CLOSE (1) 

DO  1000  1-1,5 

IK  - I 

CALL  OPNFIL(IK, IFL, JNFIL) 

NFP  - 1 

NLP  - (NFP  + NPT ( I ) ) - 1 

STORING  NECESSARY  VALUES  FOR  INDIVIDUAL  WAVEFORM 

CALL  GETWAVE ( INDX , BOUND ( I ) , I , NFP , NLP , DEL , ' A ' ,Y) 
400  KCH  - KCH  + 1 

CHAFIL(KCH)  - JNFIL 
1000  CONTINUE 
GOTO  10 
4000  CONTINUE 

WRITE(6 , *)  ' ' 

WRITE(6 , *)  '***  END  OF  PHASE  1 **** 

WRITE (6 , *)  ' ' 

GOTO  8010 
8000  CALL  CLOSE (1) 

8002  WRITE(6 , 8005) 

8005  FORMAT (3X, 'ENTER:  [Y]  TO  TRY  AGAIN’ ,/10X, 

*'  [N]  TO  TERMINATE  THE  PROCESSES'  ,/lX,  '»'  ,$) 
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READ(5 , 8007)  DUMMY 

8007  FORMAT (Al) 

IF  (DUMMY  .EQ.  'Y')  GOTO  10 
IF  (DUMMY  . EQ.  'N')  GOTO  8008 
GOTO  8002 

8008  I FLAG  2 - 10 
8010  CONTINUE 

RETURN 

END 

C 

^ A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A 

c * * 

C * SUBROUTINE  OPNFIL  * 


C 

C 

c 


c 


c 

c 

G 

C 

C 

C 


C 


* * 
A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A 


SUBROUTINE  OPNFIL(IK, IFL, JNFIL) 
DATA  CODES  /• F! , ' C( , ' Ee , ' Vs , 1 D' / 
CHARACTER  CODES ( 5 ) , IFL*60 
CHARACTER  JNFIL*30 , BLANK* 1 


JNFIL— IFL(1 : 30) 

IDOT-IND EX (JNFIL, ' . ' ) 

JNFIL ( IDOT+1 : IDOT+1 ) -CODES ( IK) 

RETURN 

END 


■k-k-k-k-k'k-k'kif-k'k-k-^-k-k-k-krk-kir-kic-k^-kieic-kie’ikick-k'kickieic-k-kirkirk-k'kicie-Jcirkick'k’k 
* * 

* SUBROUTINE  CHANNEL2  * 

* DELAYED  TIMEZ , INVERTED  TIMEZ  * 

'kJc-kicJe&Jck'kie'k-jrk’k’k’k&’feie’kic’kie-Sc’kielck'kic'k'k’k'&'kirk'k-jc’kieirk'kic'kie'kic'kie'kieicie 


SUBROUTINE  CHA2 ( El , IFLAG2 ) 

DIMENSION  Y( -1000: 10000) 

INTEGER  NPT(3) ,B0UND(4) ,E1(10) .CHANUMB 
INCLUDE  ' DISK$OVR: [ASGPROG] COMVAR.LIS ' 
CHARACTER  IFL*60, DUMMY  *1 
CHARACTER  JNFIL*30 , BLANK , IDUM 
CHARACTER  *30  CHAFIL(108) 

COMMON  /CBY/  Y 

* /CBE  / E,NPT, BOUND 

* /CBNPT/  NPTS 

* /CBFIL/  CHAFIL 

* /CBCH  / CHANUMB 


I FLAG 2 - 0 
CHANUMB  - 0 
KCH  - 0 
INDX-Q 
BLANK-'  ' 

CLOSE(l) 

WRITE(6 , 700) 

700  FORMAT (/ IX, '**  BEGINNING  OF  INDIVIDUAL  WAVEFORM  SEPARATION 
WRITE(6 , 701) 


**'  ) 
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o o o oooo  oooooo 


THE  TOTAL  NUMBER  OF  FILE  NAMES  YOU  ENTER’ , 

MUST  BE  A MULTIPLE  OF  TWO.  THE  FIRST,  THIRD, ' , 
FIFTH,  etc.  FILE  MUST  BE  AN  X- COMPONENT  FILE', 
AND  THE  SECOND,  FOURTH,  SIXTH,  etc.,  FILE', 
MUST  BE  A Z-COMPONENT  FILE.'//) 


701  FORMAT (//IX, 'NOTE: 

& /IX,' 

& /IX,' 

& /IX,' 

& /IX,' 

10  WRITE(6 , 710) 

710  FORMAT (//IX, 'ENTER  DATA  FILE  NAME  (XXXXXXXXX.XXX) ' , 

* /IX, 'IF  NO  MORE  FILES,  TYPE  " CNTR  Z"  ' , /IX,  ’ »' , $) 

READ(5 , 520  , END  - 4000)  IFL 
520  FORMAT (A60) 


READING  HEADING  FROM  THE  FILE  CONTAINING  5 WAVEFORMS 


CALL  UDSIO('R' , 'F' , IFL, 1 , Y, IERR) 

IF  (IERR  .EQ.  0 ) THEN 

ERROR  IN  INPUT  FILE 

20  WRITE (6 ,*)  ’ ERROR!!!' 

WRITE (6,*)  ' ***  FILE  CANNOT  BE  ACCESSED  ***' 

GOTO  8000 
END  IF 

CALL  EXTRACT (IFL) 

CORRECT  UDSIO  FILE 
BOUND  GENERATION 

INDX-INDX+1 

CALL  GETWAVE ( INDX ,0,0, NFP , NLP , DEL , ' I ' ,Y) 

CHANUMB  - CHANUMB  + 1 
KCH  - KCH  + 1 
CHAFIL(KCH)  - IFL(1:30) 

KCH-KCH+3 
BOUND (1) -El (1) 

IF(E1(1) .LT.l)  BOUND (1)—1 
B0UND(2)-E1(2) 

NOPTS  - (NLP  - NFP)  + 1 
BOUND  ( 3 ) —NLP 
DEL-0. 125/(E1(2)-E1(1)) 

WRITE (6 ,*) ' DEL-' , DEL 
DO  40  I -1,2 

NPT(I)-  BOUND (1+1)  - BOUND ( I ) 

40  CONTINUE 
CLOSE(l) 

DO  1000  1-1,2 

IK  - I 

CALL  OPNFIL1 ( IK , IFL , JNFIL) 

NFP  - 1 

NLP  - (NFP  + NPT(I))  - 1 

STORING  NECESSARY  VALUES  FOR  INDIVIDUAL  WAVEFORM 

CALL  GETWAVE( INDX, BOUND(I) ,I,NFP, NLP, DEL, 'A' ,Y) 
KCH  - KCH  + 1 


400 
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CHAFIL(KCH)  - JNFIL 
1000  CONTINUE 
GOTO  10 
4000  CONTINUE 

WRITE (6,*)  '***  END  OF  PHASE  1 ***' 

GOTO  8010 
8000  CLOSE (1) 

8002  WRITE(6 , 8005) 

8005  FORMAT (3X, ’ENTER:  {Y]  TO  TRY  AGAIN’ ,/10X, 

*’[N]  TO  TERMINATE  THE  PROCESSES  ' , /IX , , $) 

READ (5 , 8007)  DUMMY 

8007  FORMAT (Al) 

IF  (DUMMY  .EQ.  ’Y’)  GOTO  10 
IF  (DUMMY  .EQ.  ’N’)  GOTO  8008 
GOTO  8002 

8008  I FLAG 2 - 10 
8010  CONTINUE 

RETURN 

END 

C 

c * * 

G * SUBROUTINE  OPNFIL1  * 

C * * 

C ieikic^rkirkieic^'kirk^ieieirkieieirkieieirieiei^irkieik'kie-kiHeieieieie-kie'kie-irk'kie'k^e^s’kie 

c 

SUBROUTINE  OPNFIL1 ( IK , IFL , JNFIL) 

DATA  CODES  /’V' ,'D'/ 

CHARACTER  CODES (2) , IFL*60 
CHARACTER  JNFIL*30 , BLANK* 1 

C 

JNFIL-IFL(1 : 30) 

IDOT-INDEX( JNFIL , 1 . ’ ) 

JNFIL(IDOT+l: IDOT+1) “CODES (IK) . 

RETURN 


END 

C 

C 

c * * 

C * SUBROUTINE  CDNUM  * 

C * * 

C •kieieiiieifieieieie'kic’kirk'kicirkirkirkie'kickieieie 

C 


* 

* 


* 


SUBROUTINE  CDNUM(IFLAG2) 
CHARACTER  *4  TSTCDNUM 
CHARACTER  *30  CHAFIL(108) ,XFIL 
CHARACTER  *4  TST(18) 

CHARACTER  BLANK 
INTEGER  CHANUMB 


COMMON  BLOCKS 


C 


COMMON  /CBCH/CHANUMB 
/CBFIL/  CHAFIL 
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* /CBCDNUM/  TSTCDNUM 

C 

I FLAG 2 - 0 
KX  - 1 
KXCH  -1 
BLANK  - ' ' 

KNUMX  - CHANUMB 
5 XFIL  - CHAFIL(KXCH) 

DO  10  J -1,30 

IF  (XFIL(J : J)  .NE.  BLANK)  GOTO  20 
10  CONTINUE 

20  TSTCDNUM(1:4)  - XFIL(J+1  :J+4) 

TST(KX)  -TSTCDNUM 

KX  « KX  + 1 

KNUMX  - KNUMX  - 1 

IF  (KNUMX  .EQ.  0)  GOTO  100 

KXCH  - KXCH  + 6 

GOTO  5 

COMPARING  ALL  THE  TEST  CODE  NUMBERS 

100  KNUMX  -CHANUMB 
TSTCDNUM  - TST(l) 

DO  200  I — 1 ,KNUMX 

WRITE(6, 101)1, TST(I) 

101  FORMAT (IX, ’TEST  (',13,')  - ’ ,A4) 

IF  (TSTCDNUM  .NE.  TST(I))  GOTO  220 

00  CONTINUE 

NO  ERROR 
GOTO  300 
ERROR 

220  WRITE (6 ,*)  'ERROR!!!' 

WRITE (6,*)  '***  DATA  FILE  NAMES  ARE  FROM  MORE  THAN  ONE  TEST 
WRITE (6,*)  '**  PROGRAM  TERMINATED 

IFLAG2  - 10 
300  RETURN 
END 


SUBROUTINE  UDSIO (ROW , SOF , FILNAM , LUN, Y, IOS) 


C 

C 

C 

C 

C 

C 

C 

C 

C 


UDSIO  - UDS  File  Read  / Write  Utility  Routine 

Principal  Variables: 

ROW  -Read  or  Write  flag 

SOF  -Specs  only  or  Full  flag 

FILNAM  -UDS  File  Name 
LUN  -Logical  Unit  for  UDS  I/O 

Y -Measurement  data  array 


***' 

***' 
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IOS  -Status  and  message  suppression  flag 


INTEGER  SYS$GETMSG 
DIMENSION  Y( -1000 :*) 

CHARACTER  R0W*1 , S0F*1 , FILNAM*(*) ,M*256 

INCLUDE  ' DISK$0VR: [ASGPROG] COMVAR.LIS ' 

IMSG-IOS 

IOS-0 

IF  (ROW.EQ.’R’)  THEN 
C Read  a UDS  file 

IF  (FILNAM.NE. ' !)  OPEN (LUN , FILE-FILNAM, STATUS-' OLD' , 

* FORM- ’ UNFORMATTED ' , READONLY , IOSTAT-IOS , ERR-104) 

READ (LUN , ERR— 105) TSTSRC , TSTNUM , TSTPRF , TSTREF , TSTCFN , CLSSPD , 

* VEHNO , MAKE , MODEL , YEAR , BODY , ENGINE , VEHTWT , 

* OCCTYP , OCCAGE , OCCSEX , OCCWT , DUMSIZ , RESTR1 , RESTR2 , 

* HIC ,T1 , T2 , CLIP3M, CSI , AIS , 

* CURNUM , DATTYP , UNITS , AXIS , SENLOC , SENATT , STATUS , FORM , DESC , 

* NFP , NLP , DEL , INIVEL , PREF , FCUT , FCOR , FSTP , YO , YD , YCAL , 

* ID1 , ID2 „ ID3 , ID4 , IDS , RD1 , RD2 , RD3 , RD4  s RD5 , GDI , CD2 

IF  (SOF.NE.'S')  THEN 

G Read  Measurement  Data  Also 

IF  ( (NFP.LT. -1000) .OR. (NLP.LT„NFP) ) THEN 
IF  (IMSG.NE. -99)  WRITE(6,993) 

IOS— 3 
RETURN 
END  IF 

IF  ( (NSIZE.GT.O) .AND. (NSIZE.LT.NLP) ) THEN 
NLP— NS I ZE 
END  IF 

IF  (NLP. GT. 10000)  THEN 
NLP-10000 
END  IF 

C Read  Measurement  Data  from  Input  File 

DO  101  IR— 1 , (NLP-NFP) /1000+1 
IS-NFP+1000*(IR-1) 

READ ( LUN , ERR-10  6 ) (Y(I) , I-IS,MIN0(IS+999 ,NLP) ) 

101  CONTINUE 

IF  (NFP.GT.O)  THEN 

DO  102  I-O.NFP-l  • 

Y(I)— 0 . 0 

102  CONTINUE 

NFP— 0 

ELSE  IF  (NLP.LT.O)  THEN 
DO  103  I-NLP+1,0 


100 


Y(I)-0 . 0 

103  CONTINUE 

NLP-0 

END  IF 

IF  ( (FORM. EQ. 'DFT' ) .AND. (NFP. LT.O))  THEN 
NFP-0 
END  IF 
END  IF 
RETURN 

104  IF  (IMSG.NE. -99)  THEN 

WRITE(6 , 991) 

CALL  ERRSNS ( , IRMS , , , ) 

IF  (IRMS.NE.O)  THEN 

1ST— SYS$GETMSG ( %VAL( IRMS ) , LEN , M , , ) 

IF  (LEN.GT.O)  WRITE(6,*)M( :LEN) 

END  IF 
END  IF 
RETURN 

105  IF  (IMSG.NE. -99)  WRITE(6,992) 

IOS— 2 

RETURN 

106  IF  (IMSG.NE. -99)  WRITE(6,995) 

IOS— 4 

RETURN 

ELSE  IF  (ROW.EQ.'W)  THEN 
C Write  a UDS  File 

IF  (FILNAM.NE. ' ’)  THEN  • 

IF  (SOF.NE. 'S' ) THEN 

OPEN(LUNsFILE-FILNAM, STATUS-' NEW' , FORM- ’ UNFORMATTED ' , 

* IOSTAT— IOS , ERR— 110) 

ELSE 

OPEN (LUN,FILE-FILNAM, STATUS- 'UNKNOWN' , FORM- ’ UNFORMATTED ’ , 

* IOSTAT— IOS , ERR— 110) 

END  IF 

END  IF 

IF  (SOF.NE. 'S' ) THEN 

IF  (NFP.LT. -1000)  THEN 
NFP— 1000 
END  IF 

IF  (NLP. LT. NFP)  THEN 
NLP-NFP 
END  IF 

IF  ((NSIZE.GT.O) .AND. (NSIZE.LT. NLP))  THEN 
NLP-NSIZE 
END  IF 

IF  (NLP. GT. 10000)  THEN 
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NLP-10000 
END  IF 

IF  (NFP.GT.Q)  THEN 
DO  107  I-O.NFP-l 
Y ( I ) -0 . 0 

107  CONTINUE 

NFP-0 

ELSE  IF  (NLPoLT.O)  THEN 
DO  108  I-NLP+1 , 0 
Y(I)-0 . 0 

108  CONTINUE 

NLP-0 

END  IF 

IF  ((FORM. EQ. 'DFT' ) .AND. (NFP.LT.O))  THEN 
NFP-0 
END  IF 
END  IF 

WRITE (LUN) TSTSRC , TSTNUM , TSTPRF , TSTREF , TSTCFN , CLSSPD , 

* VEHNO , MAKE , MODEL , YEAR , BODY , ENGINE , VEHTWT , 

* OCGTYP , OCCAGE , OCCSEX , OCCWT , DUMSIZ , RESTR1 , RESTR2 , 

* HIC , T1 , T2 , CLIP3M , CSX , AIS , 

* CURNUM , DATTYP , UNITS , AXIS , SENLOC , SENATT , STATUS , FORM , DESC , 

* NFP , NLP , DEL , INIVEL , PREF , FCUT , FCOR , FSTP , YO , YD , YCAL , 

* ID1 , ID2 , ID3 , ID4 , IDS , RD1 „ RD2 , RD3 , RD4 , RD5 , GDI , CD2 

IF  (SOF.NE. ' S " ) THEN 

DO  109  IR— 1 , (NLP-NFP) /1000+1 
I S-NFP+1000* ( IR- 1 ) 

WRITE (LUN) (Y(I) , I-XS ,MIN0(IS+999 ,NLP) ) 

109  CONTINUE 
END  IF 
RETURN 

110  IF  (IMSG.NE. -99)  THEN 

WRITE(6 , 991) 

CALL  ERRSNS(,IRMS, , ,) 

IF  (IRMS.NE.O)  THEN 

1ST— SYS$GETMSG ( %VAL( IRMS ) ,LEN,M, ,) 

IF  (LEN.GT.O)  WRITE (6 ,*)M( :LEN) 

END  IF 
END  IF 
RETURN 

ELSE 

IF  (IMSG.NE. -99)  WRITE(6S996) 

IOS— 5 
RETURN 
END  IF 

991  FORMAT (/IX, '***  UDSIO  ERROR  - File  open  failure'/) 

992  FORMAT (/ IX, ' ***  UDSIO  ERROR  - Specification  data  read  error'/) 

993  FORMAT ( /IX, ’ ***  UDSIO  ERROR  - Illegal  NFP  and/or  NLP  value'/) 
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995  FORMAT ( /IX, ' ***  UDSIO  ERROR  - Measurement  data  read  error’/) 

996  FORMAT ( /IX , ’ ***  UDSIO  ERROR  - Bad  Read  or  Write  flag  value'/) 

END 

C 

C Subroutine  GETWAVE  is  a routine  which  stores  the  full  (unseparated) 
C set  of  waveforms  and  then  directs  the  portion  of  interest  into 
C the  array  Y when  called  to  do  so 

SUBROUTINE  GETWAVE ( INDX , BOUND , IWF , NFP , NLP , DEL , OT , Y) 

DIMENSION  IPT (18,5) , INFP (18,5), INLP (18,5) ,TDEL(18,5) 

DIMENSION  YFULL( 18, -1000:10000) 

DIMENSION  Y( -1000: 10000) 

INTEGER  BOUND 
CHARACTER*!  OT 

IF  (OT.EQ. 'I' .OR.OT.EQ. * i ' ) THEN 
DO  20  I -NFP , NLP 
YFULL(INDX,I)-Y(I) 

20  CONTINUE 
GOTO  99 
END  IF 

IF  (OT.EQ. 'A' .OR.OT.EQ. 'a' ) THEN 
IPT ( INDX , IWF) -BOUND 
INFP (INDX, IWF) -NFP 
INLP ( INDX , IWF) -NLP 
TDEL( INDX , IWF) -DEL 
GOTO  99 
END  IF 

IF  (OT.EQ. 'O' .OR.OT.EQ. 'o' ) THEN 
J-IPT (INDX, IWF) -1 
NFP-INFP( INDX, IWF) 

NLP- INLP ( INDX , IWF) 

DEL-TDEL( INDX, IWF) 

DO  30  I— NFP , NLP 

Y ( I ) -YFULL ( INDX , J+I ) 

30  CONTINUE 
END  IF 
99  RETURN 
END 
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The  third  segment  , "WGSUBS" , contains  the  following 


modules : 


RESULT 

HIC3 

PART 

EXTRACT 

HSINE 

ERR2 

CRASH 

FRSP3 

SGAUSS 

FAST 

FR2TR 

FR4TR 

FORD1 

FORD  2 

GETIND 

THCALC 


The  source  code  listing  for  these  modules  is: 


ooo  ddd  ddd 


SUBROUTINE  RESULT ( DFLAG , IKJNUM , CNTFL , IWF , DELFLAG ) 

C 

C * * 

C * SUBROUTINE  RESULT  * 

C * OBJECTIVE  :TAKE  THE  SQRT  (SUM  OF  SQURES  OF  X * 

C * AND  Z COMPONENT  OF  HALF- SINE  AND  * 

C * CRASH  DATA  FILE)  TO  PRODUCE  THE  * 

C * RESULTANT  FILE.  * 

C * PRINCIPAL  VARIABLES:  * 

C * NOPTS  - NO.  OF  POINTS  FOR  THE  DATA  SET  * 

C * RES  - RESULTANT  VECTOR  * 

C * FILNAM  - INPUT  FILE  NAMES  (X,Y,Z)  * 

C * OUTFIL  - OUTPUT  FILE  NAME  (RESULTANT)  * 

Q jrk,kicirje-k-jriticicjrk’k-frk-k-&-jrici(icic/eieicjcie-k'kickick*icJrk-k-kickicjcjck-k 

C 

C DECLARATION  OF  VARIABLES 
C 

INCLUDE  ' [ASGPROG]COMVAR.LIS' 

REAL  RES (-1000:4500) ,Y(-1000  : 10000) ,Y1( -1000 : 10000) 

CHARACTER  *30  FILNAM(3) 

CHARACTER  BLANK*1 , OUTFIL  *30,DUMFILE  *30, DUMMY  *1 
CHARACTER  DUMFIL1*30 
CHARACTER  AXIS3(2) 

INTEGER  NOPTS, SETNUM, CHANUMB , CNTFL, DELFLAG 
REAL  XMIN , XMAX , YMIN , YMAX 
CHARACTER  *30  CHAFIL(108) ,RESFIL(72) 

CHARACTER  DUMFILE1  *30,RESF*30 
DATA  AXIS3/'X' , 'Z'/ 

COMMON  BLOCKS 

COMMON  /CBFIL/  CHAFIL 

* /CBRES  /RESFIL 

* /CBCH/ CHANUMB 

COMMON/RSLTNTS/RESPF( 2, 9, -1000:10000) ,NFPRES(2,9) ,NLPRES(2,9) 

START  OF  PROGRAM 

DFLAG  - 0 
KJ  - 0 

KNUM  - CHANUMB 
SETNUM  - 0 
BLANK  ' 

KKJNUM-0 
ISLT-1 

IF  (IWF.EQ.5)  ISLT-2 
SUMMARY  REPORT  FILE 

IF  ((DELFLAG. EQ.l) .OR. ((DELFLAG. EQ.O) . AND . (IWF . EQ . 2) ) ) 

& OPEN (4, FILE-' RESULT. SUM1 , STATUS-' NEW' ,ERR-1000) 

WRITE  (4,400) 

400  FORMAT (/ 34X, '**  RESULTANT  SUMMARY  REPORT  **') 
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ooo  ooo  oooo  ooo  ooo  ooo 


WRITE  PROGRAM  OUTPUT  HEADER 

20  CONTINUE 

IF  (KNUM  .EQ.  0 .AND.  CNTFL  .EQ.  1)  GOTO  300 
SETNUM  - SETNUM  +1 
WRITE (6 , 21)SETNUM 

21  FORMAT (IX, 'RESULT  FILE  FOR  SET  NUMBER  ',14 ) 

BEGIN  PROCESSING 

DO  30  I— 1000, 10000 

IF  (I.LT.4501)  RES (I)  - 0.0 
RESPF( ISLT , SETNUM , I ) -0 . 0 

30  CONTINUE 

INPUT  FILE  NAMES 

DO  50  L-1,2 

IKJNUM  * IKJNUM  + 6 
FILNAM(L)  - CHAFIL (IKJNUM) 

KJ  - &J  + 1 
RESFIL(KJ)  - FILNAM(L) 

DUMFILl  - FILNAM(L) 

DO  523  I - -1000,10000 
Y(I)  - 0.0 
CONTINUE 

READ  DATA  FROM  ARRAY  (X-AXIS  ON  THE  FIRST  PASS, 

Z-AXIS  ON  THE  SECOND  PASS) 

KKJNUM-KKJNUM+1 

CALL  GETWAVE (KKJNUM , K , IWF , NFP , NLP , DEL , ' 0 ' ,Y) 
525  NOPTS  - (NLP  - NFP)  + 1 

KNUM  » KNUM  - I 
IF (NLP  .GT.  4500)  NLP-4500 
DO  70  I -NFP , NLP 

RES (I)  - RES(I)+(Y(I)*Y(I) ) 

70  CONTINUE 

CALL  CLOSE(l) 

50  CONTINUE 

PERFORM  CALCULATION 

DO  100  I— NFP , NLP 

RES (I)  - SQRT(RES (I) ) 

RESPF(ISLT, SETNUM, I)  - RES(I) 

100  CONTINUE 

NLPRES (ISLT , SETNUM) -NLP 
NFPRES (ISLT, S ETNUM) -NFP 

OUTPUT  RESULTANT  FILE 

540  K - 0 
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DO  545  J-1,26 

IF(DUMFIL1(J:J)  .EQ.  ' ')  GOTO  545 
K - K + 1 

RESF(K:K)  - DUMFIL1(J:J) 

545  CONTINUE 

DO  550  J -1,26 

IF  (RESF(J : J)  .EQ.  *.')  GOTO  555 
550  CONTINUE 

555  RESF (J+l  : J+l)  - 'R' 

RESF(J+2  : J+2)  - ’E’ 

RESF(J+3  : J+3)  - ’S’ 

OUTFIL  - RESF 
KJ  - KJ  + 1 
RESFIL(KJ)  - OUTFIL 


WRITE (4, 1100)  SETNUM 
1100  FORMAT (/2X, ’SET  NO  : \2X,I3) 

WRITE(4 , 1200)  FILNAM(l) ,FILNAM(2) 

1200  FORMAT (2X, 'INPUT  FILES (X  AND  Z COMPONENT)-:  ' ,5X, 

1 A30 , 10X, A30) 

WRITE(4, 1300)  OUTFIL 

1300  FORMAT (2X, 'RESULTANT  FILE  : '.5X.A30) 

GO  TO  20 
300  CONTINUE 

IF  (IWF.EQ.5)  CALL  CL0SE(4) 

GOTO  2000 

1000  WRITE (6 ,*)  'ERROR  !!!' 

WRITE (6,*)  ' **  "RESULT. SUM"  CANNOT  BE  CREATED  **' 
WRITE (6,*)  'PLEASE  TRY  AGAIN' 

DFLAG  - 1 
2000  CONTINUE 


C 

C ADD  ENTRIES  TO  REPORT 
C 


RETURN 

END 


C 


SUBROUTINE  HIC3 (DFLAG , CNTFL, ISLT) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


Objective: 


Major  variables: 


Programmer: 


ASGI  - S.  Mentzer 
Latest  update 
UDSIO  calls  removed 


To  find  the  maximum  Head  Injury  Criterion  value 
and  the  corresponding  time  interval  using  an 
efficient  global  branch  and  bound  algorithm. 


R(I)  - Resultant  acceleration  array,  1-0,1 NLP 

S ( I ) - Summed  resultant  array  (integral) 

JL(I)  - Lower  J boundary  index  array 

KL(I)  - Lower  K boundary  index  array 

FMAX  - Incumbent/max  value  of  HIC  subfunction  F 


7/83 

3/84 

11/86 
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C JMAX  - J index  of  FMAX  value 

C KMAX  - K index  of  FMAX  value 

C - ----------- - 

DIMENSION  R( -1000 : 10000) ,3(0:10000) ,JL(15000) ,KL(15000) 

DOUBLE  PRECISION  RC,RM,SC,SM 
CHARACTER  INFXL*30 , CURD  AT* 9 , DUMMY* 1 
INTEGER  CNTFL  , DFLAG , KNUM , CHANUMB , SETNUM 
INCLUDE  ' [ ASGPROG ] C0MVAR .LIS' 

CHARACTER  *30  RESFIL(72) 

COMMON  /CBRES/RESFIL 
* /CBCH/CHANUMB 

COMMON/RSLTNTS/RESFF( 2, 9, -1000:10000) ,NFPRES(2,9) ,NLPRES(2,9) 
COMMON/RHICS/RHIC(2 , 9) ,RT1(2,9) ,RT2(2,9) 

DATA  NAR/X50Q0/ 

C 

DFLAG  - 0 
XKJ  - 0 

KNUM  - CHANUMB 
SETNUM  - 0 
CALL  DATE(CURDAT) 

G OPEN ( 3, FX LE~ ' HI C3. REF  \ STATUS” 'NEW'  , ERR  -*  1000) 

C WRITE (3 , 931)CURDAT 

C 931  FORMAT (//IX, 5 5XS 'HIC3  REPORT' , 30X,A/1X,54X,13(1H-)/// 

C * IX, 2X, ' TSTREF’ ,5X, 'VEHICLE' ,22X, 'OCCUPANT' ,3X, 'HIC' , 

G * 4X/T1  (msec)'  ,1X,  'T2  (msec) ' , 2X, ' T2  - Tl' , 3X, ' FILE’/ 

C * 1X,10(1H-) ,2X,29(XH=) ,1X,8(1H-) ,1X,8(1H-) ,1X,9(1H-) ,1X, 

C * 9(1H-) ,1X,9(1H=) ,2X,40(1H-)) 

10  CONTINUE 

IF  (KNUM  .EQ.  0 .AND.  CNTFL  .EQ.  1)  GOTO  400 
SETNUM  - SETNUM  + 1 

WRITE  (6,*)  'CALCULATING  *HIC*  FOR  SET  NUMBER' , SETNUM 
509  IKJ  — IKJ  + 3 

INFIL  - RESFIL(IKJ) 

DO  516  I * -1000,10000 
R(X)  » RESPF(ISLT, SETNUM, I) 

516  CONTINUE 

512  KNUM  - KNUM  -2 

NFP-NFPRES ( ISLT , SETNUM) 

NLP«NLPRES ( ISLT , SETNUM) 

C Compute  sum  (integral)  array  using  trapezoid  rule 
S(0)-0.0 
SM”0 . DO 
RM”ABS (R(0) ) 

DO  102  I— 1 , NLP 
RG—ABS (R(I) ) 

SC-SM+.5*(RC+RM) 

S(I)-SC 
SM-SC 
RM-RC 
102  CONTINUE 

C Initialize  variables 

FMAX-0.0 
JMAX-0 
KMAX— 0 
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JL(l)-0 

KL(1)  — 1 

NWU-NLP 

IBR-1 

IER-1 

INC— 1 

IBL-NAR 

C Begin  a round  of  region  evaluations 
103  IL-IBL-INC 

DO  104  IR— IBR , IER , INC 
JLI-JL(IR) 

IF  (JLI.GT.0)  THEN 
NJ-NWU 
ELSE 

JLX-AES(JLI) 

NJ-MAX0(NWU-1,0) 

END  IF 
KLI-KL(IR) 

IF  (KLI.GT.O)  THEN 
NK— NWU 
ELSE 

KLI-ABS(KLI) 

NK-MAXO (NWU- 1,0) 

END  IF 
KUI-KLI+NK 

IF  ((NJ.EQ.O) .AND. (NK.EQ.O))  THEN 
C This  is  a single  point  region  - Evaluate 

F-(S(KUI)-S(JLI))*(KUI-JLI)**(- .6) 

IF  (F.GT.FMAX)  THEN 
FMAX-F 
JMAX-JLI 
KMAX-KUI 
END  IF 
ELSE 

C This  is  a finite  region  with  feasible  points  - Bound 

JUI-JLI+NJ 
FS— S (KUI ) - S ( JLI ) 

IF  (FS*MAXO(KLI-JUI, !)**(- „6) . GT.FMAX)  THEN 
C Region  is  possibly  optimal  - Evaluate  and  partition 

F-FS*(KUI-JLI)**(- .6) 

IF  (F.GT.FMAX)  THEN 
FMAX-F 
JMAX-JLI 
KMAX-KUI 
END  IF 

IF  ( (IL+4*INC)*INC . GT . IR*INC)  THEN 
WRITE(6 , 993) 

993  FORMAT ( /IX, '***  ERROR  - HIC3  array  size  exceeded'/) 

CALL  CLOSE(3) 

DFLAG  -1 
GOTO  670 
END  IF 

CALL  PART( JLI , JUI , NJ , NWU-NJ , KLI , KUI , NK , NWU-NK , IL , INC , 

* JL,KL) 
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END  IF 
END  IF 

104  CONTINUE 

IF  (ABS(IBL-IL-INC).GT.O)  THEN 

C Some  regions  remain  - Set  up  boundary  index  arrays 
ITMP-IBL 
IBL-IER 
IER-ITMP 
IBR-IL 
INC«-1*INC 
NWU-NWU/2 
GO  TO  103 
ELSE 

C Done  - No  regions  remain  - Report  results 
HIC-DEL*FMAX**2 . 5 
T1-X0QQ*JMAX*DEL 
T2-1000*KMAX*DEL 
RHIC ( ISLT , SETNUM) -HIC 
RT1 ( I SLT , S ETNUM) -T1 
RT2 ( ISLT , SETNUM) -T2 

C Reporting  of  HIC  (normally  done  here)  has  been  removed 
CALL  EXTRACT (INFIL) 

CLOSE(2) 

GOTO  10 
END  IF 
CONTINUE 
CALL  CLOSE(3) 

GOTO  670 

WRITE (6 ,*)  ’ERROR  ! ! ! ! • 

WRITE (6,*)  '***  "HIC. REP"  CANNOT  BE  OPENED  ***' 
DFIAG  - 1 
CONTINUE 
RETURN 
END 


SUBROUTINE  PART ( JLI , JUI , NJ , NJD , KLI s KUI , NK , NKD , IL , INC , JL , KL) 

Performs  partitioning  of  possibly  optimal  regions 

DIMENSION  JL(*),KL(*) 

J2M— 2*M0D(NJ , 2) 

NJ2-NJD*J2M 
JMU— (JUI-NJ/2)*(1-NJ2) 

K2M— 2*M0D ( NK , 2 ) 

NK2-NKD*K2M 
KMU-(KUI-NK/2)*(1»NK2) 

C First  region 

IL”IL+INC 
JL(IL)-JMU 
KL(IL)— KMU 

C Second  region  if  finite  width  along  J axis 
IF  (NJ.NE.O)  THEN 


400 


1000 


670 


C 

C 


C 

c 

C 
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JLI-JLI*(J2M-NJ2-1) 

IL— IL+INC 
JL(IL)-JLI 
KL(IL)-KMU 
END  IF 

IF  (NK.NE.O)  THEN 

KLI— KLI*  ( K2M  - NK2  - 1 ) 

C Third  region  if  finite  width  along  J and  K axes 
IF  (NJ.NE.O)  THEN 
I L— IL+INC 
JL(IL)-JLI 
KL(IL)— KLI 
END  IF 

C Fourth  region  if  feasible  and  finite  width  along  K axis 
IF  (ABS(KMU) -l.GT.ABS(JMU) ) THEN 
I L— IL+INC 
JL(IL)— JMU 
KL(IL)— KLI 
END  IF 
END  IF 
RETURN 
END 


SUBROUTINE  EXTRACT (FILNAM) 
CHARACTER* (*)  FILNAM 

LOC-INDEX( FILNAM, ' : ' )+l 
FILNAM-FILNAM(LOC : ) 

LOC-IND EX (FILNAM, ' ] ' )+l 
FILNAM-FILNAM(LOC : ) 

LOC-INDEX( FILNAM, • ; ' ) 

IF  (LOC.GT.O)  FILNAM-FILNAM ( : LOC - 1 ) 

RETURN 

END 


*************************************** 

* HALF- SINE  PROCESSING  * 

*************************************** 

SUBROUTINE  HSINE(DFLAG , CNTFL) 

DECLERATION  OF  VARIABLES 

INTEGER  CNTFL , CHANUMB , DFALG , SETNUM 
INTEGER  EFLAG(18 ,4) , IND(4) ,NOPTS 
CHARACTER  *4  TSTCDNUM 
CHARACTER  *30  HSNE1 , HSNE2 , HSNE3 
CHARACTER  *30  IFLNM2(18) 

CHARACTER  * 2 C(18) 

CHARACTER  DAY*9,HOUR*9 
CHARACTER  *30  RESFIL(72) 

CHARACTER  IFLNM1*30  , CHANUM*2  , DUMMY  *1 
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REAL  Y(- 1000  : 10000 ) , NPMS , YMN , YMX 
C 

C COMMON  BLOCKS 

C 

INCLUDE  ' [ASGPROG]COMVAR.LXS’ 

COMMON  /CBRES/  RESFIL 

* /CBCH/CHANUMB 

* /CBCDNUM/TSTCDNUM 

COMMON  /HSTMX/TOMX , YMXX , TOMZ , YMNZ , TCLEV1 , TCLEV2 , PCT , TDACLX , TDACLZ 
GOMMON/RSLTNTS/RESPF(2 , 9 , -1000 : 10000) ,NFPRES(2,9) ,NLPRES(2,9) 
COMMON/RHICS/RHIC (2,9) ,RT1(2,9) ,RT2(2,9) 

COMMON/HSSTAT2/HSDEVX(18) ,HSDEVZ(18) 

COMMON/HSSTAT3/HSHDEV ( 18 ) 

C 

C HALF-SINE  WAVEFORM  PROCESSOR 
C 

DFXAG  - 0 
SETNUM  - 0 
KMIM  - CHANUMB 
N—Q 

IKJ  — 0 
KKJ-0 

DO  30  1-1,18 
DO  30  J— 1,4 
EFLAG ( I , J ) —0 
30  CONTINUE 
C 

C REPORT  FILE 

C 

HSNE1  - ’ HSNE ’ //TSTCDNUM (1:4)//' .SUM' 

HSNE2  - ’ HSNE ’ //TSTCDNUM (1 : 4) // • . 001 1 
HSNE3  HSNE ’//TSTCDNUM ( 1:4)//' .002’ 

OPEN (3 , FILE-HSNE1 , STATUS-’ NEW* , ERR-3000) 

OPEN ( 7, FXLE-HSNE2, STATUS- 'NEW' , ERR-3000) 

OPEN (8 , FILE-HSNE3 , STATUS-' NEW’ , ERR-3000) 

WRITE(3 , 210) 

210  FORMAT ( 43X , ’**  HALF-SINE  WAVEFORM  PROCESSOR  DATA  SUMMARY  **'// 

1 IX, ’THE  * SYMBOL  INDICATES  THAT  THE  ALLOWABLE  LIMIT  HAS  ' 

2 , 'BEEN  EXCEEDED. '//IX, ’TIME  QUANTITIES  ARE  EXPRESSED  IN  ' 

3 , ’MILLISECONDS’ ) 

CALL  DATE (DAY) 

CALL  TIME (HOUR) 

WRITE(3 , 9995)  DAY , HOUR 
9995  FORMAT (//, SOX, A9,10X,A9) 

700  CONTINUE 

IF  (KNUM  .EQ.  0 .AND.  CNTFL  .EQ.  1)  GOTO  800 
SETNUM  - SETNUM  + 1 

WRITE (6 , *)  ’PROCESSING  SET  NUMBER  SETNUM 
DO  35  1-1,4 
35  IND(I)— ' ' 

C 

C PROCESS  HALF- SINE  "X"  COMPONENT 
C 

701  IKJ  - IKJ  + 1 
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KKJ-KKJ+1 

IFLNM1  - RESFIL(IKJ) 

705  DO  706  1—1000,10000 

Y(I)  - 0.0 

706  CONTINUE 
C 

C READING  "X"  COMPONENT  INPUT  FILE (DATA, HEAD) 

C 

CALL  GETWAVE(KKJ,0,2,NFP,NLP,DEL, 'O' ,Y) 

6 N-N+2 

KNUM  - KNUM  - 1 
IFLNM2(N-1)— IFLNM1 
C 

C GETTING  THE  CHANNEL  NUMBER 
C 

DO  900  I - 1,30 

IF  (IFLNM1(I : I)  .EQ.  GOTO  910 

900  CONTINUE 

910  CHANUM(1 : 1)  - IFLNMl(I+2  : 1+2) 

CHANUM(2 : 2)  - IFLNMl(I+3  : 1+3) 

C(N-l)  - CHANUM 
NPMS-1 . 0/DEL/1000 . 

C 

C DETERMINE  ARRAY  NO.  WHICH  CORRESPONDS  TO  THE  MINIMUM  VALUE  IN  DATA 
C 

NOPTS  - (NLP  - NFP)  + 1 
PCT-20 . 0 

CALL  GETIND  (NFP , NLP , PCT , Y , DEL , NX , TOMX , DNX , YMXX , TDACLX) 

C 

C PROCESS  HALF- SINE  "Z"  COMPONENT 
C 

CALL  CLOSE (1) 

801  IKJ  - IKJ  + 1 
KKJ-KKJ+1 

IFLNM1  - RESFIL(IKJ) 

805  DO  807  I —1000,10000 
Y(I)  - 0.0 
807  CONTINUE 
C 

C READING  ”Z”  COMPONENT  INPUT  FILE (DATA, HEAD) 

C 

CALL  GETWAVE (KKJ ,0,2, NFP , NLP , DEL , 'O' ,Y) 

18  IFLNM2 (N)-IFLNMl 
KNUM  - KNUM  - 1 
C 

C GETTING  THE  CHANNEL  NUMBER 
C 

I— INDEX( IFLNM1 , ' . ' ) 

960  CHANUM(1 : 1)  - IFLNMl(I+2  : 1+2) 

CHANUM(2 : 2)  - IFLNMl(I+3  : 1+3) 

C(N)  -CHANUM 
C 

C DETERMINE  ARRAY  NO.  WHICH  CORRESPONDS  TO  THE  MAXIMUM  VALUE  IN  DATA 
C 
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NOPTS  - (NLP  - NFP)  +1 
PCT-20.0 

CALL  GETIND (NFP , NLP , PCT , Y,  DEL , NZ , TOMZ , DNZ , YMNZ , TDACLZ ) 

C 

C COMPUTE  TIME  DEVIATION  FROM  THEORETICAL  AND  CHANNEL -TO -CHANNEL  TIME  DIFFERENCE 
C 

CALL  CLOSE (1) 

DIF— TOMX- TOMZ 
CCTD— DIF/1000 . 

IF(ABS(DIF) .LE.0.1)  GOTO  200 
EFLAG(N-1 t 1)— 1 
EFLAG(N,1)-1 
IND ( 1 ) - 1 *' 

200  DEVX— TOMX- 12 . 5 
DEVZ-TOMZ-12 . 5 

IF(ABS(DEVX) .GTol.O)  EFLAG(N-1 , 2)-l 
IF(ABS(DEVZ) .GT.1.0)  EFLAG(N,3)»1 
IF(ABS(DEVX) .GTol.O)  IND(2)-'  *' 

IF(ABS(DEVZ) .GTol.O)  IND(3)-'  *' 

WRITE (3 , 11)  IFLNM2 (N-l) ,C(N-1) ,IFLNM2(N) ,C(N) 

11  FORMAT ( //IX , "RESULTS  FOR  "X"  COMPONENT  FILE  : ' ,A30,5X, 

1 ' CHANNEL  NO . 6 , A2 , /9X, ' AND  "Z"  COMPONENT  FILE  : • , 

2 A30 » 5X s 9 CHANNEL  NO . • ,A2) 

WRITE (3, 13)  TSTPRF 

13  FORMAT ( /IX , 8 FACILITY  : 9 ,A25) 

WRITE (3  , 1000)  DIF, XND(l)  • 

1000  FORMAT (/IX, 8 1 „ CHANNEL - TO - CHANNEL  TIME  DIFFERENCE  - 8 , F7 . 3 , 1A2 , 

1 IX, ' (ALLOWABLE  LIMIT  - +/-  °-l  MSEC)') 

WRITE (3, 1010)  DEVX , IND ( 2) 

'1010  FORMAT (/IX, 9 2.  "X™  COMPONENT  TIME  DEVIATION  - 8,F7.3,1A2, 

1 IX, ! (ALLOWABLE  - +/-  1.0  MSEC)8) 

WRITE(3 , 1020)  DEVZ, IND(3) 

1020  FORMAT (/IX,  '3.  "Z*8  COMPONENT  TIME  DEVIATION  - ' .F7.3.1A2, 

1 IX, ' (ALLOWABLE  LIMIT  - +/-  1.0  MSEC)') 

HSDEVX ( SETNUM) -DEVX 
HSDEVZ ( SETNUM) -DEVZ 
WRITE (7,1001) DEVX , DEVZ 

1001  FORMAT ( 2X , F7 . 3 , 5X , F7 . 3 ) 

810  IKJ  - IKJ  +1 

IFLNM1  - RESFIL(IKJ) 

820  DO  825  1—1000,10000 

Y(I)  - RESPF(1 , SETNUM, I) 

825  CONTINUE 

HIC-RHIC(1, SETNUM) 

T1-RT1 ( 1 , SETNUM) +100 

T2-RT2 ( 1 , SETNUM) +100 

T2T1-T2-T1 

HICEXP— 0 . 0 

TT1-0 . 0 

TT2-0 . 0 

CALL  THCALC( CCTD, HICEXP, TT1.TT2) 

T21EXP-TT2-TT1 

HICDEV— (HIC/HICEXP-1)*100 . 0 

IF(ABS (HICDEV) .LE.6.0)  GOTO  51 
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EFLAG(N-1,4)-1 
EFLAG(N,4)-1 
IND(4)— ' *' 

51  WRITE(3 , 31)  HICDEV , IND(4) 

31  FORMAT (/IX, '4.  HIC  DEVIATION  = ' ,F6.2,lA2, 

1 * (ALLOWABLE  LIMIT  - +/-  6%)') 

HSHDEV ( SETNUM) -HICDEV 
WRITE(8, 1002) HICDEV 
1002  FORMAT  (2X-.F6. 2) 

WRITE(3 , 110)  HIC,T1,T2,T2T1 

110  FORMAT (/IX, 'THE  CALCULATED  HIC  NO.  - • ,F6.1,'  FOR  T1  - ' , F7 . 3 , 
1 ' AND  T2  - ' , F7 . 3 , ' WITH  T2-T1  - ' ,F7.3) 

WRITE (3 , 20)  HICEXP , T21EXP 

20  FORMAT (/IX, 'THE  THEORETICAL  HIC  NO.  - ',F6.1,'  AND  THE  ', 

1 'THEORETICAL  VALUE  OF  T2-T1  - ',F7.3) 

GOTO  700 

800  WRITE(3 , 1060) 

1060  FORMAT (//5X, 'FILE  NAME' , 19X, ’ CHANNEL  NO .', 20X, ' STATUS ' ) 

NA-0 

NU-0 

DO  160  I— 1 ,N 
KND-0 

DO  170  J-1,4 

IF ( EFLAG ( I , J ) . EQ . 0 ) GOTO  170 
KND-1 

170  CONTINUE 

IF(KND.EQ.O)  NA-NA+1 
IF(KND.EQ.l)  NU-NU+1 

IF(KND.EQ.O)  WRITE(3 , 1070)  IFLNM2(I) ,C(I) 

IF(KND.EQ.l)  WRITE(3 , 1080)  IFLNM2(I) ,C(I) 

1070  FORMAT ( IX, A30,10X,A2,25X, 'PASS' ) 

1080  F0RMAT(1X, A30 , 10X, A2 , 25X, ' FAIL' ) 

160  CONTINUE 

WRITE(3 , 1050)  N.NA.NU 

1050  FORMAT (//IX, 'TOTAL  NO.  OF  CHANNELS  PROCESSED  - ',12 

1 /IX, 'NUMBER  OF  ACCEPTABLE  CHANNELS  - ',12 

2 /IX, 'NUMBER  OF  UNACCEPTABLE  CHANNELS  - ',12) 

CALL  CLOSE(3) 

GOTO  9999 

3000  WRITE (6 ,*)  'ERROR  !!!' 

WRITE (6,*)  ' **  "HSINE . SUM"  CANNOT  BE  CREATED  **' 

9999  CONTINUE 

CALL  CLOSE(7) 

CALL  CLOSE(8) 

RETURN 

END 

C 

SUBROUTINE  ERR2 (DUMMY) 

C 

CHARACTER  DUMMY* 1 
C 

CALL  CLOSE(l) 

8002  WRITE(6 , 8005) 

8005  FORMAT (3X, 'ENTER:  [Y]  TO  TRY  AGAIN' ,/10X, 
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* ' [N]  TO  TERMINATE  THE  PROCESSES  ' , /IX,  , $) 

READ(5 , 8007)  DUMMY 
8007  FORMAT (Al) 

IF  (DUMMY  .EQ.  'Y'  .OR.  DUMMY  .EQ.  'N')  GOTO  8010 
GOTO  8002 
8010  RETURN 
END 


SUBROUTINE  CRASH ( DFLAG , CNTFL) 

VARIABLE  DECLARATION 

INTEGER  EFLAG ( 18 ) , IND , NOPTS , DFLAG , CNTFL 
INTEGER  S ETNUM , CHANUMB 
CHARACTER  *30  CRASHFI1 , CRASHFI2 
CHARACTER  *4  TSTCDNUM 
CHARACTER  *30  IFLNM2(18) 

CHARACTER  *2  C(18) 

CHARACTER  IFLNM1  *30,CHANUM  *2 , DUMMY  *1 
REAL  Y(- 1000  : 10000 ) , NPMS , YMN , YMX 
CHARACTER  *30  RESFIL(72) 

CHARACTER'  DAY*98HOUR*9 

INCLUDE  ' [ASGPROG j COMVAR .LIS' 

COMMON  /CBRES/RESFIL 

* /CBCH/CHANUMB 

* /CBCDNUM/TSTCDNUM 

COMMON/RSLTNTS/RESPF(2, 9, -1000:10000) ,NFPRES(2,9) ,NLPRES(2,9) 
COMMON/RHICS/RHIC (2*9) ,RT1(2,9) ,RT2(2S9) 

COMMON/CSTAT2/CHDEV ( 18 ) 

CRASH  WAVEFORM  PROCESSOR 

DFLAG  - 0 
S ETNUM  * 0 
KNUM  - CHANUMB 
IKJ  * 0 
KKJ-0 
N-0 

DO  30  1-1,18 
EFLAG ( I ) — 0 
30  CONTINUE 

REPORT  FILE 

CRASHFI1  -' CRASH ' //TSTCDNUM (1:4)//' .SUM' 

CRASHFI2  -’ CRASH' //TSTCDNUM( 1:4)//' .001’ 

OPEN ( 3, FILE-CRASHFI1, STATUS-' NEW' , ERR-  3000) 

OPEN (7 , FILE-CRASHFI2 , STATUS-' NEW , ERR—  3000) 

WRITE(3 , 210) 


C 

C 

c 

C 

c 


* SUBROUTINE  CRASH  * 


irkie&ic&ie'k'k'kirk'kirk'SeJeicjcirk'frJcie’kieirk'k-Jck-kirkick'k'kiekickickiefeTk’k’kjeik-ickic 
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210  FORMAT ( 43X, '**  CRASH  WAVEFORM  PROCESSOR  DATA  SUMMARY  **'// 

1 IX, 'THE  * SYMBOL  INDICATES  THAT  THE  ALLOWABLE  LIMIT  HAS 

2 , ' BEEN  EXCEEDED. '//IX, 'TIME  QUANTITIES  ARE  EXPRESSED  IN 

3 , 'MILLISECONDS' ) 

CALL  DATE (DAY) 

CALL  TIME (HOUR) 

WRITE(3 , 9995)  DAY, HOUR 
9995  FORMAT (// , 50X, A9 , 10X, A9) 

700  CONTINUE 

IF  (KNUM  .EQ.  0 .AND.  CNTFL  .EQ.  1)  GOTO  800 
SETNUM  - SETNUM  + 1 

WRITE (6 ,*)  'PROCESSING  SET  NUMBER  '.SETNUM 
35  XND-'  ' 

C 

C READ  CRASH  "X"  COMPONENT  DATA 
C 

701  IKJ  — IKJ  + 1 
KKJ-KKJ+1 

IFLNM1  - RESFIL(IKJ) 

705  DO  706  I —1000,10000 

Y(I)  - 0.0 

706  CONTINUE 
C 

C READ  X-AXIS  DATA  FROM  ARRAY 
C 

CALL  GETWAVE (KKJ ,0,5, NFP , NLP , DEL , ' 0 ' , Y) 

6 N-N+2 

KNUM  - KNUM  -1 
IFLNM2 (N- 1)-IFLNM1 
C 

C GETTING  THE  CHANNEL  NUMBER  FROM  INPUT  FILE 
C 

I-INDEX(IFLNM1 , ' . • ) 

910  CHANUM(1 : 1)  - IFLNMl(I+2  : 1+  2) 

CHANUM(2 : 2)  - IFLNMl(I+3  : 1+  3) 

C(N-l)  - CHANUM 
CALL  CLOSE (1) 

NPMS— 1 . 0/DEL/1000 . 

C 

C READ  CRASH  "Z"  COMPONENT  DATA 
C 

921  IKJ  - IKJ  + 1 
KKJ-KKJ+1 

IFLNM1  - RESFIL(IKJ) 

925  DO  926  1—1000,10000 

Y(I)  - 0.0 

926  CONTINUE 
C 

C READ  Z-AXIS  DATA  FROM  ARRAY 
C 

CALL  GETWAVE (KKJ, 0,5, NFP, NLP, DEL, 'O' ,Y) 

60  IFLNM2 (N)-IFLNMl 
KNUM  - KNUM  -1 
C 
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GETTING  THE  CHANNEL  NUMBER 

I— INDEX(IFLNM1 , ' . ' ) 

9600  CHANUM(1 : 1)  - IFLNMl(I+2  : 1+2) 

CHANUM(2 : 2)  - IFLNMl(I+3  : 1+3) 

C(N)  - CHANUM 
980  CONTINUE 

CALL  CLOSE (1) 

WRITE (3, 11)  IFLNM2(N-1) ,C(N-1) , IFLNM2(N) , C(N) 

11  FORMAT (//IX, 'RESULTS  FOR  "X"  COMPONENT  FILE  : • ,A30,5X, 

1 'CHANNEL  NO.  ',A2/9X,’AND  "Z"  COMPONENT  FILE  : ', 

2 A30.5X, 'CHANNEL  NO.  ' ,A2) 

WRITE (3, 13)  TSTPRF 

13  FORMAT (/IX, ’FACILITY  : • ,A25) 

950  XKJ  — IKJ  + 1 

IFLNMl-  RESFIL(IKJ) 

960  DO  961  I —1000,1000 

Y(I)  - RESPF(2 , SETNUM, I) 

961  CONTINUE 

HIC— RHIG ( 2 , SETNUM) 

T1-RT1 ( 2 , SETNUM)+375 
T2-RT2 (2 , SETNUM) +3 7 5 
T2T1-T2-T1 
HICEXP-930 

HICDEV— (HIG/HICEXP-1) *100 . 0 
IF(ABS (HICDEV) .LE.6.0)  GOTO  51 
EFLAG(N-l)-! 

EFLAG(N)=1 
IND-’  *' 

51-  WRITE(3 , 31)  HICDEV , IND 

31  FORMAT ( /IX,  '1.  HIC  DEVIATION  - \F10.2,lA2, 

1 ' (ALLOWABLE  LIMIT  » +/-  6%)’) 

CHDEV (SETNUM) -HICDEV 
WRITE (7, 1061) HICDEV 
1061  FORMAT (2X.F10. 2) 

WRITE(3 , 110)  HIC,T1,T2,T2T1 

110  FORMAT (/IX, 'THE  CALCULATED  HIC  NO.  - ',F6.1,’  FOR  T1  - ’,F7.3, 
1 ' AND  T2  - ' , F7 „ 3 , ' WITH  T2-T1  - ' ,F7.3) 

WRITE(3 , 20)  HICEXP 

20  FORMAT (/IX, 'THE  THEORETICAL  HIC  NO.  - ' ,F6.1) 

GOTO  700 

800  WRITE(3 , 1060) 

1060  FORMAT (//5X, 'FILE  NAME ',1 9X, ' CHANNEL  NO 20X, ' STATUS ' ) 

NA-Q 

NU-0 

DO  160  I— 1 ,N 
KND-0 

IF(EFLAG(I) .EQ.O)  GOTO  170 
KND-1 

170  IF(KND.EQ.O)  NA-NA+1 
IF(KND.EQ.l)  NU-NU+1 

IF(KND.EQ.O)  WRITE(3 , 1070)  IFLNM2(I) , C(I) 

IF(KND.EQ.l)  WRITE(3 , 1080)  IFLNM2 (I) , C(I) 

1070  FORMAT ( IX, A30 , 10X,A2 , 25X, 'PASS’ ) 
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1080  FORMAT ( IX, A30 , 10X, A2 , 25X, 'FAIL' ) 

160  CONTINUE 

WRITE(3 , 1050)  N , NA , NU 

1050  FORMAT (//IX, 'TOTAL  NO.  OF  CHANNELS  PROCESSED  - ',12 

1 /IX,' NUMBER  OF  ACCEPTABLE  CHANNELS  - ',12 

2 /IX, 'NUMBER  OF  UNACCEPTABLE  CHANNELS  - ' ,12) 

CALL  CLOSE(3) 

CALL  CLOSE (7) 

GOTO  9999 

3000  WRITE(6,*)  'ERROR!!!1 

WRITE (6 , *)  ' ***  "CRASH. SUM"  CANNOT  BE  CREATED  ***' 

9999  CONTINUE 
RETURN 
END 
C 

SUBROUTINE  FRSP3 (DFLAG , CNTFL) 

C 

C PROGRAM  FRSP3:  John  E.  Nickles 

C Transportation  Systems  Center 

C 23-October-1984 

C 

C Last  Revision  Date:  3-October-1985 

C 

C Abstract: 

C 

C This  program  calculates  the  frequency  response  of  a system 

C when  the  input  to  the  system  is  a sum- of- sines  waveform.  The 

C frequency  response  is  determined  by  calculating  the  ratio  of  two 

C DFT  outputs  at  the  signal  frequencies.  The  two  DFTs  are  the  DFT 

C of  the  recorded  system  output  and  the  DFT  of  the  reconstructed 

C system  input.  The  input  to  the  program  is  a raw,  time  series  data 

C set  of  the  recorded  system  output  signal.  The  system  input  data  can 

C be  recreated  if  the  waveform  specifications  and  the  sampling  fre- 

C quency,  with  which  the  system  output  data  set  was  made,  are  known. 

C The  data  is  windowed  before  the  DFTs  are  calculated.  The  program 

C uses  the  FAST  FFT  algorithm  to  compute  the  DFTs . 

C 

C Subroutines  Called: 

C 

C FAST  FFT 

C S GAUSS  Gaussian  window 

C 

C Revision  History: 

C 

C ll-September-1985  Added  scale  factor,  CSC 

C 3-October-1985  Integrated  FRSP3  with  waveform  generator  signal 

C processing  software  from  MGA  Research  Corp. 

C 

C 

c 

REAL  A( -1000 : 10000) , B(2050) , C(2050) , AR(20) , F(20) ,Y(20) ,Z(20) , 

* FDFT(20) , ARDB(20) ,U(0:1024) ,V(0:1024) ,W(2048) 

INTEGER  LN(20) 

CHARACTER* 3 0 OUTFIL , FRFIL , SUMFIL , STTFIL , CHAFIL( 108 ) 
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CHARACTER*9  DAY , HOUR , DUMMY* 1 
INTEGER  CNTFL, DFLAG , CHANUMB , SETNUM 
C 

COMMON  /CBFIL/CHAFIL 
* /CBCH/CHANUMB 

COMMON/SSSTAT2/SSF(18 , 14) , SSARDB(18 , 14) 

C 

INCLUDE  ’ [ASGPROG ] COMVAR .LIS' 

C 

XKJ  - -1 
KNUM  - CHANUMB 
SETNUM  - 0 
NN  - 2048 

G 

IF (mJM.EQ.O. AND. CNTFL. EQ.l) GOTO  999 
610  XKJ— XKJ -8-6 

SUMFXL- ' SUMSN * //CHAFXL (XKJ) (2:5)//' .SUM’ 

STTFIL* ’ SUMSN ' //CHAFXL (IKJ) (2 : 5)//* .001' 

IKJ-IKJ-6 

C 

650  OPEN ( 3 , STATUS— ’ NEW 8 , FXLE-SUMFXL) 

OPEN (7, STATUS- 'NEW' , FXLE-STTFIL) 

C 

CALL  DATE (DAY) 

CALL  TIME (HOUR) 

WRXTE( 3,1) DAY, HOUR 
1 F0RMAT(1X,A9 , 2X,A9 ,/) 

C 

700  CONTINUE 

IF  (KNUM  .EQ.  0 .AND.  CNTFL  .EQ.  1)  GOTO  800 

701  XKJ  - IKJ  + 6 
OUTFXL  - CHAFIL(XKJ) 

J— INDEX (OUTFIL, • . 9 ) 

712  FRFIL  -• FR1 ’//CHAFXL (IKJ) (2: 5)//’ . '//OUTFIL (J+l:J+3) 
750  OPEN (21, STATUS-’ NEW ' .FXLE-FRFIL) 

C 

C Read  output  signal  data  set. 

C 

XNDX-l+INT ( FLOAT ( XKJ ) /6 . 0 ) 

CALL  GETWAVE ( XNDX ,0,4, NFP , NLP , DEL , 'O' ,A) 

760  KNUM  - KNUM  -1 

N -(NLP  - NFP)  + 1 
DT  - DEL 

IF  (N  .GT.  NN)  N-NN 
L - 14 

FSG  - 273.4375 
ALPHA  - 9 
PI— 4 . *ATAN ( 1 . ) 

C 

C Index  translation 

C 

DO  5 I-O.NN+l 
5 B(I+1)-A(I) 

C 
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Recreate  input  waveform  data  file. 

OMEGA-2 . *PI*FSG 
DO  4 I— 1 , N 

T—DT* FLOAT ( I - 1 ) 

C(I)-0. 

DO  2 J-l.L 

X— S IN ( J *OMEGA*T ) 

2 C(I)-C(I)+X 

4 CONTINUE 

Zero  fill 

IF(N.LT.NN)GO  TO  8 
GO  TO  12 

8 DO  10  I— N+l , NN+2 
B ( I ) -0 . 

10  C(I)-0. 


Window  the  data 


12  CALL  SGAUSS(W,N, ALPHA) 
DO  14  I— 1 ,N 
B(I)-W(I)*B(I) 

14  C(I)-W(I)*C(I) 


Compute  DFTs 


CALL  FAST(B.NN) 

CALL  FAST(C ,NN) 

DO  20  1-0  sNN/2 

U ( I ) -SQRT ( C ( 2*1+1 ) *C ( 2*1+1 ) +C ( 2*1+2 ) *C ( 2*1+2 ) ) 
20  V ( I ) -SQRT ( B ( 2*1+1 ) *B ( 2*1+1 ) +B ( 2*1+2 ) *B ( 2*1+2 ) ) 


Compute  frequency  response 


30 


40 


DF— l./(NN*DT) 
DO  30  J— 1 , L 
FS— J*FSG 
RN-FS/DF 
K-INT(RN) 
R-RN-K 


! DF— frequency  increment  of  analyzing  DFT  output 
! J-th  harmonic  is  analyzed  during  J-th  pass. 

! FS  - frequency  of  J-th  harmonic 

! K - DFT  output  line  no.  corresponding  to  FS 
! R - remainder  fraction  after  truncation 


IF (R . GT . 0 . 5 ) K-K+l 

LN(J)-K 

Y(J)-U(K) 

Z(J)-V(K) 

F(J)-FS 

FDFT(J)-K*DF 

CONTINUE 

CSC— (0.993) *Y(1)/Z(1) 

DO  40  J— 1 , L 

AR(J)-CSC*Z(J)/Y(J) 

! AR  — ratio  of  output  DFT  to  input  DFT  at  J-th  harmonic. 
ARDB(J)-20.*ALOG10(AR(J))  ! frequency  response  in  dB  at  FS 
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WRITE (21, 150)  (F(J)  , ARDB(J)  ,J=*1,L) 

Write  summary  file. 

WRITE (3, 500)  OUTFIL, FRFIL 

500  FORMAT (////10X, 'RESULT  FOR  FILE  ' , 5X, A30 , /10X, 

* 'FREQUENCY  RESPONSE  FILE  1 , 5X, A30 ,///) 

WRITE(3 , 220) 

WRITE(3 , 221) 

WRITE(3 , 222) 

WRITE (3 , 230) (F(I) ,LN(I) , FDFT(I) ,Y(I) ,Z(I) ,AR(I) ,ARDB(I) , I— 1 ,L) 

Store  statistical  data  file. 

WRITE(7,509)(ARDB(I) ,1-1, L) 

509  FORMAT ( IX , <L>F9 .4) 

DO  510  1-1,  L 

SSF(INDX,I)-F(I) 

S S ARDB ( INDX , I ) -ARDB ( I ) 

510  CONTINUE 

WRITE(6 , 601)  OUTFIL 

601  FORMAT (IX, 'PROCESSING  COMPLETED  FOR  FILE  \5X,A30) 

CLOSE(21) 

GOTO  700 

End  of  loop 

800  CL0SE(3) 

CL0SE(7) 

999  RETURN 
120  FORMAT (A) 

150  FORMAT (2E13. 5) 

220  F0RMAT(/,3X, 'SIGNAL* ,4X, * DFT" , 6X, 1 DFT' ,7X, 'INPUT' ,6X, 'OUTPUT' , 

* 4X, ’SCALED8 ,2X, 'AMPLITUDE' ) 

221  FORMAT ( ' FREQUENCY' ,2X, 'LINE' ,4X, 'OUTPUT' ,7X, ’DFT' ,9X, ’DFT' , 

* 7X, 'DFT' ,5X, 'RATIO' ) 

222  FORMAT (4X, ' (HZ) ' ,5X, ’NO. ' , 2X , 8 FREQUENCY ’ ,5X, 'VALUE' ,7X, 'VALUE' , 

* 5X, 'RATIO' ,5X, ' (DB) • ,/) 

230  FORMAT ( 20 ( IX , F9 . 4 , 2X , 14 , 2X , F9 . 4 , 2X , 2 ( F10 . 4 , 2X) , F7 . 5 , 3X , F7 . 2 , :/)) 


END 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 


Subroutine  SGAUSS 


Subroutine  to  generate  a Gaussian  data  window 


John  E.  Nickles 
Transportation  Systems  Center 
14- August- 1984 


SUBROUTINE  SGAUSS (W,L, ALPHA) 
DIMENSION  W(L) 
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DO  5 I— 1 , L 

D— (2 . *(I-1)/L) -1 . 

5 W ( I ) —EXP ( - ( ALPHA*D*ALPHA*D ) /2 ) 

RETURN 

END 

C 

C 

C SUBROUTINE : FAST 

C REPLACES  THE  REAL  VECTOR  B(K) , FOR  K-1,2, . . . ,N, 

C WITH  ITS  FINITE  DISCRETE  FOURIER  TRANSFORM 

C 

c 

SUBROUTINE  FAST(B , N) 

C 

C THE  DC  TERM  IS  RETURNED  IN  LOCATION  B(l)  WITH  B(2)  SET  TO  0. 
C THEREAFTER  THE  JTH  HARMONIC  IS  RETURNED  AS  A COMPLEX 
C NUMBER  STORED  AS  B(2*J+1)  + I B(2*J+2) . 

C THE  N/2  HARMONIC  IS  RETURNED  IN  B(N+1)  WITH  B(N+2)  SET  TO  0. 
C HENCE,  B MUST  BE  DIMENSIONED  TO  SIZE  N+2 . 

C THE  SUBROUTINE  IS  CALLED  AS  FAST(B ,N)  WHERE  N-2**M  AND 
C B IS  THE  REAL  ARRAY  DESCRIBED  ABOVE. 

C 

DIMENSION  B(N+2) 

COMMON  /CONS/  PII,  P7 , P7TWO,  C22,  S22,  PI2 
IW  IS  A MACHINE  DEPENDENT  WRITE  DEVICE  NUMBER 
IW  - 6 

PII  - 4 . *ATAN ( 1 . ) 

PI8  - PII/8 . 

P7  - 1 . /SQRT ( 2 . ) 

P7TWO  - 2.*P7 
C22  - COS (PI8) 

S22  - SIN(PI8) 

PI2  - 2.*PII 
DO  10  1-1,15 
M - I 
NT  - 2**1 

IF  (N.EQ.NT)  GO  TO  20 
10  CONTINUE 

WRITE  (IW,  9999) 

99  FORMAT  (33H  N IS  NOT  A POWER  OF  TWO  FOR  FAST) 

STOP 

20  N4P0W  - M/2 

DO  A RADIX  2 ITERATION  FIRST  IF  ONE  IS  REQUIRED. 

IF  (M-N4POW*2)  40,  40,  30 
30  NN  - 2 

INT  - N/NN 

CALL  FR2TR(INT,  B(l) , B(INT+1)) 

GO  TO  50 
40  NN  - 1 
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PERFORM  RADIX  4 ITERATIONS. 

50  IF  (N4POW.EQ.O)  GO  TO  70 
DO  60  IT-1, N4P0W 
NN  * NN*4 
INT  - N/NN 

CALL  FR4TR(INT,  NN,  B(l) , B(INT+1) , B(2*INT+1) , B(3*INT+1) , 
* B(l) , B(INT+1),  B(2*INT+1),  B(3*INT+1)) 

60  CONTINUE 

PERFORM  IN-PLACE  REORDERING. 

70  CALL  F0RD1(M,  B) 

CALL  FORD2(M,  B) 

T - B(2) 

B(2)  - 0. 

B(N+1)  - T 
B(N+2)  - 0. 

DO  80  IT-4, N, 2 
B(XT)  - -B(IT) 

80  CONTINUE 
RETURN 
END 
C 

C-- -------------------------------------------------------  --------- 

C SUBROUTINE:  FR2TR 

C RADIX  2 ITERATION  SUBROUTINE 


C 

SUBROUTINE  FR2TR(INT,  BO,  Bl) 

DIMENSION  B0(INT+1),  B1(INT+I) 

DO  10  K— 1 , INT 

T - BQ(K)  + B1(K) 

B1(K)  - BO (K)  - B1(K) 

B0(K)  - T 
10  CONTINUE 
RETURN 
END 
C 

C- ---------------------------- 

C SUBROUTINE:  FR4TR 

C RADIX  4 ITERATION  SUBROUTINE 
C-------- - 

c 

SUBROUTINE  FR4TR(INT»  NN,  BO,  Bl,  B2 , B3,  B4,  B5,  B6,  B7) 
DIMENSION  L(15) , B0(NN*INT+2) , Bl(NN*INT+2) , B2(NN*INT+2) , 
*B3(NN*INT+2) , B4(NN*XNT+2) , B5 (NN*INT+2) , B6(NN*INT+2) , 

*B7 (NN*INT+2) 

COMMON  /CONS/  PII,  P7 , P7TWO,  C22,  S22,  PI2 

EQUIVALENCE  0-15, L(l)),  (L14,L(2)),  (L13,L(3)),  (L12,L(4)), 

* (Lll , L(5) ) , (L10 , L(6) ) , (L9 , L(7) ) , (L8,L(8)),  (L7,L(9)), 

* (L6 , L(10) ) , (L5 , L(ll) ) , (L4 , L(12) ) , (L3,L(13)),  (L2,L(14)), 

* (LI , L(15) ) 
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c 

C JTHET  IS  A REVERSED  BINARY  COUNTER,  JR  STEPS  TWO  AT  A TIME  TO 
C LOCATE  THE  REAL  PARTS  OF  INTERMEDIATE  RESULTS,  AND  JI  LOCATES 
C THE  IMAGINARY  PART  CORRESPONDING  TO  JR. 


10 

20 

30 

40 


50 


60 


70 


L(l)  - NN/4 
DO  40  K-2,15 

IF  (L(K-l) -2)  10,  20,  30 
L(K-l)  - 2 
L(K)  = 2 
GO  TO  40 
L(K)  - L(K-l)/2 
CONTINUE 

FIOVN  - PII/FLOAT (NN) 


JI 

JL 

JR 

DO 

DO 

DO 

DO 

DO 

DO 

DO 

DO 

DO 

DO 

DO 

DO 

DO 

DO 

DO 


- 3 

- 2 

- 2 

120 

120 

120 

120 

120 

120 

120 

120 

120 

120 

120 

120 

120 

120 

120 


J1-2,L1,2 
J2-J1,L2,L1 
J3-J2,L3,L2 
J4-J3,L4,L3 
J5-J4,L5,L4 
J6-J5,L6,L5 
J7-J6,L7,L6 
J8-J7,L8,L7 
J9-J8,L9,L8 
J10— J9 ,L10 ,L9 
J11-J10 , Lll , L10 
J12-J11,L12,L11 
J13-J12 , L13 , L12 
J14-J13 , L14 , L13 
JTHET— J 14 , L15 , L14 
TH2  - JTHET  - 2 
IF  (TH2)  50,  50,  90 
DO  60  K— 1 , INT 

TO  - B0(K)  + B2(K) 

T1  - B1(K)  + B3(K) 
B2(K)  - B0(K)  - B2(K) 
B3(K)  - B1(K)  - B3(K) 
B0(K)  - TO  + T1 
B1(K)  - TO  - T1 
CONTINUE 

IF  (NN-4)  120,  120,  70 
KO  - INT*4  + 1 
KL  - KO  + INT  - 1 
DO  80  K-KO.KL 

PR  - P7*(B1(K) -B3(K) ) 
PI  - P7*(B1(K)+B3(K) ) 
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B3(K)  - 
B1(K)  - 
B2(K)  - 
BO (K)  = 
CONTINUE 
GO  TO  120 


B2(K)  + PI 
PI  - B2(K) 
BO (K)  - PR 
BO (K)  + PR 
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90 


100 


110 


120 


ARG  - TH2*PI0VN 
Cl  - COS (ARG) 

51  - SIN(ARG) 

C2  - Cl**2  - Sl**2 

52  - C1*S1  + C1*S1 
C3  - C1*C2  - S1*S2 

53  - C2*S1  + S2*C1 
INT4  - INT*4 

JO  - JR*INT4  + 1 
KO  - JI*INT4  -l-  1 
JLAST  - JO  + INT  - 1 
DO  100  J»J0, JLAST 
K - KO  + J - JO 
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R1 

R5 

T2 

T6 

T3 

T7 

TO 

T4 


B1(J)*C1 
B1(J)*S1 
B2(J)*G2 
B2(J)*S2 
B3(J)*C3 
B3(J)*S3 
B0(J)  * 
B4(K)  + 


T2  * B0(J)  - 


T6 

T1 

T5 

T3 

T7 


- B4(K)  - 

- R1  + T3 


+ 

T2 

T6 

T2 

T6 


- R3 

- R1 

- R5 


+ T7 

- T3 

- T7 


BQ(J) 

- 

TO 

+ 

T1 

B7(K) 

- 

T4 

+ 

T5 

B6(K) 

* 

TO 

- 

T1 

B1(J) 

- 

T5 

- 

T4 

B2(J) 

- 

T2 

- 

T7 

B5(K) 

- 

T6 

+ 

T3 

B4(K) 

- 

T2 

+ 

T7 

B3(J) 

T3 

- 

T6 

CONTINUE 

JR  - JR 

2 

JI  - JI 

- 

2 

B5(K)*S1 

B5(K)*C1 

B6(K)*S2 

B6(K)*C2 

B7(K)*S3 

B7(K)*C3 


IF  (JI-JL)  110,  110,  120 
JI  - 2*JR  - 1 
JL  - JR 
CONTINUE 
RETURN 
END 


C--  -----------  ------------------- 

C SUBROUTINE:  FORD1 

C IN-PLACE  REORDERING  SUBROUTINE 
C-- ----- ------------- 


SUBROUTINE  FQRD1(M,  B) 
DIMENSION  B(2**M+2) 


K 
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KL  - 2 
N - 2**M 
DO  40  J— 4,N, 2 

IF  (K-J)  20,  20,  10 
10  T - B(J) 

B(J)  - B(K) 

B(K)  - T 
20  K - K - 2 

IF  (K-KL)  30,  30,  40 
30  K - 2*J 

KL  - J 
40  CONTINUE 
RETURN 
END 
C 

C - --------- 

C SUBROUTINE:  FORD2 

C IN-PLACE  REORDERING  SUBROUTINE 

C --------- - 

c 

SUBROUTINE  FORD2(M,  B) 

DIMENSION  L(15) , B(2**M+2) 

EQUIVALENCE  (L15,L(1)),  (L14,L(2)),  (L13,L(3)),  (L12,L(4)), 

* (Lll ,L(5) ) , (L10 , L(6) ) , (L9 ,L(7) ) , (L8,L(8)),  (L7,L(9)), 

* (L6 ,L(10) ) , (L5 ,L(11) ) , (L4,L(12) ) , (L3,L(13)),  (L2,L(14)), 

* (LI , L(15) ) 

N - 2**M 

L(l)  - N 
DO  10  K-2,M 

L(K)  - L(K-l)/2 
10  CONTINUE 

DO  20  K— M , 14 
L(K+1)  - 2 
20  CONTINUE 
IJ  - 2 

DO  40  J1-2,L1,2 
DO  40  J2-J1,L2,L1 
DO  40  J3-J2.L3 ,L2 
DO  40  J4— J3 ,L4,L3 
DO  40  J5-*J4,L5,L4 
DO  40  J 6— J 5 , L6 , L5 
DO  40  J7-J6 , L7 , L6 
DO  40  J8-J7,L8,L7 
DO  40  J 9-J  8 , L9 , L8 
DO  40  J10-J9,L10,L9 
DO  40  Jll— J10 , Lll , L10 
DO  40  J12-J11 , L12 , Lll 
DO  40  J13— J12 , L13 , L12 
DO  40  J14-J13 , L14 , L13 
DO  40  JI— J14 , L15 , L14 
IF  (IJ-UI)  30,  40,  40 
30  T - B(IJ-l) 

B(IJ-l)  - B(JI-l) 

B(JI-l)  - T 


Reprinted  with  permission  of  the  IEEE  Press 
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T - B(IJ) 

B(IJ)  - B(JX) 

B(JI)  - T 
40  XJ  - XJ  + 2 
RETURN 
END 

C Subroutine  GETXND  determines  the  index  of  the  data  array  which  is 
C closest  to  the  calculated  time  of  occurrence  of  the  maximum  data 
C value.  The  approach  taken  is  to  average  the  times  the  data  passes 
C through  a certain  level  which  is  expressed  as  a percentage  of  the 
C theoretical  maximum  (200  g's).  Interpolation  is  used  to  improve 
C the  time  calculation.  It  is  assumed  that  noise  riding  on  the  data 
C will  not  affect  this  calculation.  (This  should  be  true  as  long  as 
C the  percentage  chosen  places  the  level  in  an  area  of  the  curve 
C which  has  a steep  slope.) 

C 

SUBROUTINE  GETXND (NFP , NLP , PCT , Y , DEL , XOMAX , TOMAX , DIOMAX , VMAX , TDACL) 
REAL  Y( -1000  : 10000) 

C 

chklev-200 . 0*pct/100 . 0 
televl-0 . 0 
tclev2~0 . 0 
do  100  X-nfp,nlp-l 
til-(I+l)*del*1000. 
ti-I*del*1000. 

if  (abs(y(X)) .le.chklev.and.abs(y(I+l)) .ge.chklev)  then 
tclevl-ti+( ( (chklev-abs(y(X) ) )/(abs(y(I+l) ) -abs(y(X) ) ) ) 

& *(til-ti) ) 

end  if 

if  (abs(y(I+l) ) . le . chklev. and. abs(y(I) ) .ge.chklev)  then 
tclev2-ti+( ( (abs(y(I) ) -chklev)/(abs(y(I) ) -abs(y(I+l) ) ) ) 

& *(til-ti)) 

end  if 
100  continue 

tomax- ( tclevl+tclev2)/2 . 0 
tdacl-tclev2 - tclevl 
do  200  J-nfptnlp-X 
tjX-(J+X)*del*X000. 
tj -J*de 1*1000 . 

if  ( t j . le . tomax . and . t j 1 . ge . tomax)  then 
dtj -tomax- tj 
dt j 1-t  j 1 - tomax 

C If  the  calculated  time  of  maximum  falls  exactly  between  two  samples, 

C the  index  of  maximum  is  arbitrarily  assigned  to  the  earlier  sample, 
if  ( dt j . le . dt j 1 ) then 
iomax-J 
diomax-dtj 
end  if 

if  ( dt j 1 . 1 t . dt j ) then 
iomax-J+1 
diomax— -dtj 1 
end  if 

vmax— y(j  ) + ( (dtj/ (tj  1- tj  ) )*(y(j+l)  -y(j ) ) ) 
end  if 
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200  continue 
return 
end 


C 

C Subroutine  THCALC  obtains,  through  multiple  steps,  the  value  of  the 
C theoretical  HIC  if  the  curve  were  to  be  made  up  of  two  components 
C (x  and  z)  which  have  the  same  amplitude  and  duration  of  the  actual 
C data.  (The  parameters  necessary  for  determining  the  duration  and 
C amplitude  are  obtained  in  subroutine  GETIND.)  The  steps  taken  to 
C obtain  the  theoretical  HIC  value  are: 

C 1)  determine  the  equations  of  the  two  half- sine  curves  which  pass 

C through  the  points  determined  via  GETIND 

C 2)  generate  arrays  of  data  which  are  the  equivalent  of  the 
C equations  should  they  be  sampled  at  the  interval  of  DEL 

C 3)  obtain  the  resultant  of  the  two  arrays 

C 4)  calculate  HIC,  Tl,  and  T2  using  a modified  version  of  the  HIC3 

C subroutine 

C 

SUBROUTINE  THCALC (CCTD , THIC , TT1 , TT2) 

COMMON  /HSTMX/TOMX , YMXX , TOMZ , YMNZ , TCLEV1 , TCLEV2 , PCT , TDACLX , TDACLZ 
DIMENSION  YX( -1000 : 10000) ,YZ( -1000 : 10000) ,R( -1000 : 10000) 

DIMENSION  S(0:10000) ,JL(15000) ,KL(15000) 

DOUBLE  PRECISION  RC,RM,SC,SM 
INCLUDE  ' [ AS G PROG ] COMVAR .LIS’ 

DATA  NAR/15000/ 

C 

pi-3.1415927 
do  10  1—1000,10000 
yx(I)-0 . 0 
yz(I)-0.0 
r ( I ) -0 . 0 
10  continue 

ix01-INT(. 01/del) 

tx01-del*ix01 

tzOl-txOl+cctd 

C From  A - X * Sin(w  * t)  and  knowing  the  times  and  values  of  the 
C points  at  which  a certain  level  is  crossed, 

C w - (2  * arc=cosine(pct/100) )/dt 

wx-(2 . 0*ACOS (pct/100 . 0) )/( tdaclx/1000 . ) 
wz-(2 . 0*ACOS (pct/100 . 0) )/(tdaclz/1000 . ) 
icflagx-0 
icflagz-0 
do  20  l-0,nlp 
t-FLOAT ( I ) *de 1 
yx(I)-0.0 
yz(I)-0 . 0 

if  (t.ge.txOl.and.icflagx.eq.O)  then 
yx(I)-ymxx*sin(wx*(t- txOl) ) 
if  (yx(I) .gt.0.0)  then 
yx(I)-0.0 
icflagx-1 
end  if 
end  if 

if  ( t . ge . tzOl . and. icflagz . eq . 0)  then 
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yz(I)-ymnz*sin(wz*(t- tzOl) ) 
if  (yz(I) .It. O.O)  then 
yz(I)-0.0 
icflagz-1 
end  if 
end  if 

r(I)»SQRT((yx(I)**2)+(yz(I)**2)) 

20  continue 

C The  following  routine  is  taken  from  the  subroutine  HIC3 

C The  code  is  repeated  here  to  eliminate  the  need  to  write  UDS  files 

C 

C Compute  sum  (integral)  array  using  trapezoid  rule 
S(O)-0.O 
SM-0.D0 
RM—ABS (R(0) ) 

DO  102  1-1 , NLP 
RC-ABS(R(I>) 

SC-SM+.5*(RC+RM) 

S(I)-SC 

SM«SG 

RM-RC 

102  CONTINUE 

C Initialize  variables 
FMAX«0 „ 0 
JMAX=0 
KMAX=0 
JL(l)-0 
KL(1)— 1 
NmJ-NLP 
IBR-1 
IER-1 
INC— 1 
1BL-NAR 

C Begin  a round  of  region  evaluations 

103  IL-IBL-INC 

DO  104  IR-IBR.IER.INC 
JLI— JL(IR) 

IF  (JLI.GT.0)  THEN 
NJ-NWU 
ELSE 

JLI-ABS(JLI) 

NJ-MAX0(NWU-1,0) 

END  IF 
KLI-KL(IR) 

IF  (KLI.GT.Q)  THEN 
NK-NWJ 
ELSE 

KLI-ABS(KLI) 

NK-MAX0 (NWU- 1,0) 

END  IF 
KUI-KLI+NK 

IF  ((NJ.EQ.0) .AND. (NK.EQ.0))  THEN 

C This  is  a single  point  region  - Evaluate 

F-(S(KUI)-S(JLI))*(KUI-JLI)**(-.6) 
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IF  (F.GT.FMAX)  THEN 
FMAX-F 
JMAX-JLI 
KMAX-KUI 
END  IF 
ELSE 

C This  is  a finite  region  with  feasible  points  - Bound 
JUI— JLI+NJ 
FS— S (KUI ) - S ( JLI ) 

IF  (FS*MAX0(KLI-JUI,l)**(-.6) .GT.FMAX)  THEN 
C Region  is  possibly  optimal  - Evaluate  and  partition 

F-FS* (KUI -JLI )**(-. 6) 

IF  (F.GT.FMAX)  THEN 
FMAX-F 
JMAX-JLI 
KMAX-KUI 
END  IF 

IF  ((IL+4*INC)*INC.GT.IR*INC)  THEN 
WRITE(6 , 993) 

993  FORMAT ( /IX , ' ***  ERROR  - HIC3  array  size  exceeded’/) 

CALL  CL0SE(3) 

DFLAG  -1 
GOTO  670 
END  IF 

CALL  PART ( JLI, JUI, NJ .NWU-NJ ,KLI , KUI ,NK,NWU-NK, IL, INC, 

* JL,KL) 

END  IF 
END  IF 

104  CONTINUE 

IF  (ABS (IBL-IL-INC) .GT.O)  THEN 
C Some  regions  remain  - Set  up  boundary  index  arrays 
ITMP-IBL 
IBL-IER 
IER-ITMP 
IBR-IL 
INC—  1*INC 
NWU-NWU/2 
GO  TO  103 
ELSE 

C Done  - No  regions  remain  - Report  results 
THIC— DEL*FMAX**2 . 5 
TT1— 1000*JMAX*DEL 
TT2— 1000*KMAX*DEL 
END  IF 

670  CONTINUE 
RETURN 
END 
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modules : 


The  fourth  segment,  ”REC"„  contains  the  following 


RECTAN 
ZERO  2 
COMPRE 


The  source  code  listing  for  these  modules  is: 
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C * * 

C * SUBROUTINE  FOR  RECTANGLE  WAVEFORM  * 

C * * 

Q ***1ln^^fr*'!^***•5fc'5^r}lr**Vr**^>fc,*■****■st*,****■*Vr*****Vr***■*•******•*'*‘'*'*, 

c 

C Subroutine  rectan  performs  signal  processing  calculations 
C upon  the  rectangular  waveform  that  has  been  recorded.  This 

C program  performs  four  time  checks  and  two  amplitude  checks  on 
C the  rectangle  waveform  output.  All  rectangle  waveforms  are 

C processed  simultaneously. 

C 

C Rectan  then  outputs  both  upper  and  lower  overshoot,  attenuation, 

C and  timeshift.  Overshoots  and  steady- state  values  are  determined 

C by  subroutine  zer02. 

C 

SUBROUTINE  RECTAN ( I CHAN , DFLAG , CNTFL) 

C 

C define  data  types 

C 

CHARACTER  IFLNMl  *30,CHANUM  *2, DUMMY  *1 
INTEGER  Nl(13) ,N2(13) ,N3(13) ,N4(13) ,NOPTS, DFLAG 
INTEGER  NCHAN , ZIC , ZDC , ZERINC ( 100 ) , ZERDEC ( 100 ) , ZFLAG , 

* 0LFLAG,LB(14) ,TZIC,UBND,LBND, START, FINISH, 

* ING0V(13) ,IPS0V(13) ,TB,TR,CHANUMB, CNTFL 
REAL  XMIN , XMAX , YMIN , YMAX , TO 

REAL  Y(- 1000  : 10000) , DAN (10000) 

REAL  NPMS , POSOV ( 13 ) , NEGOV ( 13 ) , YMN , YMX , PAVGOV , NAVG0V , 

* SSPAV(13) ,SSNAV(13) ,AVG,TL(13) ,TU(13) 

CHARACTER  *30  CHAFIL(108) 


C 

C 

C 


INCLUDE  ' DISK$0VR: [ASGPR0G] COMVAR.LIS ' 


define  common  block  symbolic  names 


COMMON  /CBXY/  XMIN , XMAX , YMIN , YMAX 

* /CBTO/  TO 

* /CBNPTS/  NOPTS 

* /CBY/  Y 

COMMON  /OVER/ 

/RANGE/ 


POSOV , NEGOV , IPSOV , INGOV , SSPAV , SSNAV 

FINISH 

TL,TU 

N1,N2,N3,N4 
START , LB 

ZIC , ZDC , ZERINC , ZERDEC , ZFLAG 
EFLAG 


/SHFT/ 

/STPT/ 

/PTN/ 

/COUNT/ 

/FLAG/ 

/CBFIL/  CHAFIL 
/CBCH/  CHANUMB 
C0MM0N/RSTAT1/RP0S (18 , 10) ,RNOS(18,10) , RSSPAV(18 , 10) , 

& RSSNAV(18 , 10) , RTU(18 , 10) ,RTL(18,10) ,RFLNM1(18) , 

& RSATT(18) ,RCNUM(18) ,RTPERF(18) 

CHARACTER* 30  RFLNM1 
CHARACTER*25  RSATT , RTPERF 
CHARACTER* 2 RCNUM 
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DIMENSION  5SNMP(13) ,SSNMN(13) 

C 

C Rectangle  waveform  processor. 

C prepare  to  input  data  and  begin  processing. 

C 

IKJ  - -4 
KNUM  -CHANUMB 
DFIAG  - 0 
WRITE(6 , 1) 

1 FORMAT (/IX, ' **  RECTANGLE  WAVEFORM  PROCESSOR  **') 

I CHAN  - 0 

C 

C START  OF  THE  LOOP  FOR  PROCESSING 

C 

700  CONTINUE 

IF  (KNUM  .EQ.  0 .AND.  CNTFL  .EQ.  1)  GOTO  800 
IKJ  - IKJ  + 6 
YMAX  - 200.0 
YMIN  - -200.0 
DO  8700  I— 1000, 10000 
Y(I)  - 0.0 
8700  CONTINUE 

5000  IFLNMX  - CHAFIL(IKJ) 

C 

C READING  FROM  DATA  FILE  CONTAINING  WAVEFORM 

C 

5001  INDX-1+INT( FLOAT ( IKJ )/6.0) 

CALL  GETWAVE ( INDX ,0,1, NFP , NLP , DEL , !0’ ,Y) 

J-l 

DO  8705  I— NFP , NLP 
DAN(J)— Y(I) 

J-J+l 

8705  CONTINUE 

DO  8710  I™1 , NLP+1 
Y(X)-DAN(I) 

8710  CONTINUE 

KNUM  - KNUM  - 1 
I CHAN  - ICHAN  + 1 
NPMS-1 . 0 /DEL/1000 . 

C 

C The  first  processing  performed  determines  the  zero  crossing 
C of  the  waveform.  When  the  waveform  crosses  zero  from  negative  to 
C positive  the  point  number  is  stored  in  the  zerinc(I)  array. 

C When  the  waveform  crosses  zero  from  positive  to  negative  the  point 

C number  is  stored  in  the  zerdec(I)  array.  The  point  numbers  are  only 

C assigned  to  the  two  arrays  above  if  there  is  a definite  sign  change. 

C It  is  not  assigned  for  a positive  to  zero  to  positive  change 

C in  the  waveform.  All  crossings  are  stored  sequencially 

C in  their  respective  arrays. 

C 

ZIC-0 

ZDC-0 

NOPTS— ( NLP - NFP ) +1 
DO  55  1-1,100 
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ZERINC ( I )-0 
55  ZERDEC(I)-0 

DO  11  I-l.NOPTS 
IF(Y(I) ) 20,30,40 

20  ZFLAG--1 
GO  TO  50 

30  ZFLAG-0 
GO  TO  50 

40  ZFLAG— 1 

50  IF  (I.EQ.l)  GOTO  100 

60  IF ( ZFLAG. EQ.OLFLAG)  GOTO  100 
IF(OLFLAG.EQ.O)  GOTO  100 
IF(OLFLAG.EQ. -1)  GOTO  90 
ZDC— ZDC+1 

ZERDEC ( ZDC ) -I 
GOTO  100 
90  ZIC-ZIC+1 

ZERINC(ZIC)— I 
100  OLFLAG— Z FLAG 

11  CONTINUE 

IF  (Y(NOPTS)  .LT.  0.0  .AND.  Y(NOPTS)  .GT.  -250.0)  GOTO  8715 
GOTO  8720 

8715  ZIC  - ZIC  + 1 
ZERINC ( ZIC ) — NOPTS 
8720  CONTINUE 

IF  (ZIC  .GT.  10)  ZIC  - 10 
C 

C The  second  processing  performed  is  the  division  of  the  waveform 
C into  10  separate  sections.  The  point  numbers  are  assigned  to  the 
C LB(I)  array.  LB(1)  is  assumed  to  occur  EXACTLY  10  msec  before  LB(2) . 

C LB(2)  and  each  successive  LB(I)  corresponds  to  the  end  of  a complete 

C period  of  a rectangle  waveform. 

C 

DO  15  1-1,13 
LB ( I ) -0 
15  CONTINUE 
IMAX-NOPTS 
TZIC-ZIC 

IMIN— ZERINC (TZIC) 

YTST— . 5*YMAX 

31  DO  12  I— IMIN , IMAX 

IF(Y(I) .GT.YTST)  GOTO  21 

12  CONTINUE 
IMAX-ZERINC (TZIC) - 1 
TZIC— TZIC- 1 

IMIN— ZERINC (TZIC) 

GOTO  31 

21  UBND-TZIC 
TZIC— 1 
IMIN— 1 

IMAX-ZERINC (TZIC) 

YTST- . 5*YMIN 

61  DO  41  I— IMAX , IMIN , - 1 
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IF(Y(I) .LT.YTST)  GOTO  51 
41  CONTINUE 

IMIN-ZERINC ( TZIC ) +1 
TZIC-TZIC+1 
IHAX— ZERXNC(TZIC) 

GOTO  61 
51  LBND-TZXC 
J-LBND-1 

START-10- (UBND-LBND) 

FXNXSH-10 
START-2 
FXNXSH-ZIC 

DO  70  X— START , FINXSH+1 
J-J+l 

LB ( I ) — ZERXNC ( J ) 

70  CONTINUE 
LB(12)—NOPTS 
C 

C The  third  processing  performed  on  the  rectangle  waveform  determines 
G the  steady- state  averages  and  overshoots  for  each  of  the  20 

C sections  of  the  rectangle  waveform.  The  positive  steady- state  values 

C are  stored  in  the  SSPAV(X)  array.  The  negative  steady-state  values 

C are  stored  in  the  SSNAV(X)  array.  The  positive  overshoot  values  are 

C stored  in  the  POSOV(X)  array.  The  negative  overshoot  values  are 

C stored  in  the  NEGOV(X)  array.  The  average  positive  and  negative 

C values  for  overshoot  are  calculated  and  stored  as  PAVGOV  and 
C NAVGOV  respectively. 

C 

C Initialize  outputs  as  default  values 

C 

DO  105  1-1,13 
SSPAV ( I ) — 1 . 0E+4 
SSNAV(I)  — 1.0E+4 
POSOV ( X ) — 1 . OE+4 
NEGOV(I)  — 1.0E+4 
105  CONTINUE 

LB(1)«LB(2)-(10.0*NPMS) 

DO  6 I-l.ZDC 

XF(ZERDEC(X) .GT.LB(START) ) GOTO  6 
KK-I 

6 CONTINUE 
K-KK 

DO  110  1-1, FINXSH+1 
LL-LB(I) 

LU— ZERDEC(K) 

XF(LU.EQ.O)  LU-NOPTS 
K-K+l 

NPT— LU-LL+1 

CALL  ZER02(Y(LL) ,NPT, SSPAV(I) ,P0S0V(I) ,0,NPMS) 

110  CONTINUE 
K-KK 

JK-FINISH 
DO  120  1=1 , JK 
LL-ZERDEC(K) 
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K-K+l 
LU— LB ( 1+1 ) 

NPT— LU-LL+1 

CALL  ZER02(Y(LL) ,NPT, SSNAV(I) ,NEGOV(I) ,1,NPMS) 

120  CONTINUE 
PAVGOV-O . 0 
NAVGOV-O . 0 

DO  130  I-START-1, FINISH 
PAVGOV-PAVGOV+POSOV ( I ) 

IF( I. EQ. FINISH)  GOTO  130 
NAVGOV-NAVGOV+NEGOV ( I ) 

130  CONTINUE 

PAVGOV-PAVGOV/ ( FINI SH - START+2 ) 

NAVGOV— NAVGOV / ( FINI SH - START+2 ) 

150  NEGOV ( 13 ) — 1 . OE+4 
C 

C The  fourth  processing  performed  determines  the  time  shifts  in  the 

C waveforms  leading  edges.  All  edges  are  compared  to  the  first 

C leading  edge  which  is  the  first  rectangle  waveform  of  the  10  rectangle 

C waveforms  expected  for  each  channel  recorded.  Time  shifts  are  not 

C compared  between  channel- to -channel  variations. 

C 

C The  timeshifts  of  the  positive  section  are  stored  in  the  TU(I)  array. 

C Timeshifts  for  the  negative  sections  are  stored  in  the  TL(I)  array. 

C 

DO  7 1-1,10 
TL(I)  — 1.0E+4 
TU(I)— 1.0E+4 
7 CONTINUE 

YTST1-0 . 9*YMIN 
YTST2-0 . 9*YMAX 
TL(1)-(LB(1)-1)/NPMS 
LS TART- START - 1 

TU ( LSTART ) - ( ZERDEC ( LSTART ) - 1 ) /NPMS - ( 6+ ( LSTART - 1 ) *10 ) 

DO  300  I-START, FINISH 

IF(POSOV(I) .EQ.1.0E+4)  GOTO  310 
TL(I)— ( (LB(I) -1)/NPMS) - (10*(I-1) ) 

310  IF(NEG0V(I) .EQ. -l.QE+4)  GOTO  300 
Kl-KK+I-1 

TU(I)— ( (ZERDEC(Kl) -1)/NPMS) - (6+(I-l)*10) 

300  CONTINUE 
C 

C The  fifth  processing  performed  determines  the  normalized  steady- state 

C and  overshoot  values  for  each  of  the  20  sections . These  values  are 

C stored  over  the  values  in  the  already  existing  NEGOV, POSOV, 

C SSNAV  and  SSPAV  arrays. 

C 

C SSNM  is  a "normalization  factor"  which  takes  into  account  the 
C fact  that  the  starting  point  for  most  changes  in  the  waveform 
C is  not  at  zero  but  rather  at  the  last  steady- state  level'. 

SSNMP(l)— 0 . 0 
SSNMN(l)— SSPAV(l) 

DO  409  1-2, FINISH 
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SSNMP(I)— SSNAV(I-l) 

SSNMN(I)-SSPAV(I) 

409  CONTINUE 

DO  420  1-1, FINISH 

IF  (POSOV(I) .EQ.1.0E+4)  GOTO  410 
IF  (SSNMP(I) .EQ. -1.0E+4)  SSNMP(I)-0.0 

POSOV(I)— (POSQV(I) =SSPAV(I) )/(SSPAV(I) -S SNMP (I) )*100 . 0 
SSPAV(I)— (SSPAV(I) -200 . 0)/2 . 0 

410  IF(NEGOV(I) .EQ. -l.QE+4)  GOTO  420 
SSNM— SSPAV(I) 

IF  (SSNMN(I) .EQ.1.0E+4)  SSNMN(I)»Q.O 

NEGOV(I)— (NEGQV(I) -SSNAV(I) )/ (SSNAV (I) -SSNMN(I) )*100 . 0 
SSNAV(I)-(SSNAV(I)+200.0)/-2.0 
420  CONTINUE 

C 

C The  data  is  now  stored  until  it  is  needed  for  comparisons 
C with  allowable  limits.  The  data  are  the  upper  and  lower 
C overshoots,  the  upper  and  lower  attenuations,  and  the  upper 
C and  lower  timeshif ts . 

C 

RFLNMl ( INDX) -IFLNM1 
RSATT ( INDX) — SENATT 
X— INDEX ( IFLNM1 , 5 . ') 

910  GHANUM— IFLNM1 ( 1+2 : 1+3 ) 

RCNUM ( INDX) -CHANUM 
RTPERF ( INDX) -TSTPRF 
DO  520  K-1,10 

RPOS ( INDX , K) -POSOV ( K) 

RNOS ( INDX , K) -NEGOV (K) 

RS S PAV ( INDX , K) -S SPAV ( K) 

RS  SNAV ( INDX , K) -SSNAV ( K) 

RTU(INDX,K)-TU(K) 

RTL(INDX,K)-TL(K) 

520  CONTINUE 

CALL  CLOSE(l) 

WRITE (6 ,*)  'PROCESSING  COMPLETED  FOR  FILE  ' , IFLNM1 
GOTO  700 
800  CONTINUE 
RETURN 
END 
C 

SUBROUTINE  ZERO 2 (X , NPTS , Z, AVG, IFL, XNPM) 

C 

DIMENSION  X(2Q00) ,Y(2000) 

K— NPTS -XNPM 
IF (IFL. EQ. 1)  GOTO  520 
AVG— l.E+4 
DO  500  I— 1 ,K 
500  AVG— AMAX1 (X ( I ) , AVG ) 

GOTO  510 
520  AVG-l.E+4 

DO  530  I-l.K 
530  AVG— AMIN1(X(I) ,AVG) 

510  GTST— ABS (X ( 1 ) ) 
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IGTST-NPTS 

Y(l)-X(l) 

DO  30  1-2 , IGTST 
DO  40  INNER— 1 ,1-1 
IF(X(I) . GT.Y( INNER) ) GOTO  50 
40  CONTINUE 

Y(I)-X(I) 

GOTO  70 
50  CONTINUE 

DO  60  K— I , INNER+1 , - 1 
Y(K)— Y(K-l) 

60  CONTINUE 

Y(INNER)-X(I) 

70  CONTINUE 
30  CONTINUE 

TEST-1. E10 
IW— IGTST/2 
NEND-IGTST-IW 
DO  10  I— 1 , NEND 
D— ABS(Y(I) -Y(I+IW) ) 

IF(D.GT.TEST)  GOTO  10 
TEST-D 
INDEX— I 
10  CONTINUE 

Z-Y ( INDEX+ ( IW+1 ) /2 ) 

K— NPTS-XNPM 

LEA-NPTS/2 

SUM-0. 

DO  707  I-LEA.K 
SUM— SUM+X(I) 

* 

707  CONTINUE 

Z— SUM/(K-LEA+1) 

RETURN 

END 

C 

C *'***'***•**■*•****'**'***■**■*************■***■*■*•'*•*■***•*■***'**■*■**■*•*■*■*■)!:'*•* 

C * * 

C * SUBROUTINE  FOR  COMPARING  RECTANGLE  WAVEFORM  RESULTS  * 

C * * 

c ******l^***************,^**'*'***************'******-5t***'***-*-)!:*-* 

C 

C Subroutine  Compre  receives  as  input  the  overshoots,  attenuations,  and 

C timeshifts  from  the  rectangle  waveform  that  was  processed.  These  values 

C were  output  from  the  rectan  subroutine.  These  values  are  compared 
C to  their  allowable  limits  and  the  results  determine  if  a particular 

C channel  is  acceptable  for  recording.  All  channels  are  processed 

C sequencially. 

C 

C The  first  comparison  performed  compares  the  deviation  of  actual  time 
C from  theoretical  time  and  allowable  limits. 

C The  result  is  stored  as  EFLAG(1,I)  where  I-channel  number.  There 

C are  a total  of  five  more  comparisons  for  each  channel.  A table  is 

C provided  for  clarity. 
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c 
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c 

c 

c 

c 

c 
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c 
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c 

c 
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c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 


c 

c 

c 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk  kkkkkkkkkkk  kkkkkkkkkkkk  kk  kkk kkkkk 
k k 

comparison  * result  of  status  * description  of  comparison 
number  * is  stored  as  * 

* k 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


k 

* 

* 

* 


* 

* 

k 

* 


EFLAG  (1,1) 


* * 

* deviation  of  actual  time  from  * 

* theoretical  time  and  allowable* 

* limit  * 


* 

* 

* 


* 

* 

* 


EFLAG  (2,1) 


* 

* 

* 


deviation  of  slope  from  ideal  * 
slope  of  zero  * 


* 

* 

* 


* 

k 

* 


EFLAG  (3,1) 


* 

* 

* 


deviation  of  intercept  from 
ideal  of  zero 


* 

k 

* 


* 

* 

* 


* 

* 

* 


EFLAG  (4,1)  * attenuation  deviation 

* 


k 

k 

k 


k 

k 

k 


k 

k 

k 


EFLAG  (5,1) 


* 

* 

k 


overshoot  deviation 


* 

* 

* 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
k k k k 

*6  * EFLAG  (6,1)  * channel- to -channel  time  lag  * 


SUBROUTINE  COMPRE(XCHAN,DFLAG) 


identify  data  types 

INTEGER  EFLAG (6, 18) ,IND(25) ,DFLAG 

REAL  0(18,26) ,A(18 , 26) ,T(18 , 26) , SLOPE, INTCP,DEV(26) 
CHARACTER  *30  RECFI1 , RECFI2 , RECFI3 , RECFI4 , RECFI5 
CHARACTER  *4  TSTCDNUM 
CHARACTER  *2  C(18) 

CHARACTER  *9  DAY, HOUR 

DEFINE  COMMON  BLOCKS 

COMMON  /CBCBNUM/TSTCDNUM 

C0MM0N/RSTAT1/RP0S(18 , 10) ,RNOS(18,10) ,RSSPAV(18 , 10) , 

& RSSNAV(18 , 10) ,RTU(18 , 10) ,RTL(18,10) ,RFLNM1(18) , 

& RSATT(18) ,RCNUM(18) ,RTPERF(18) 

C0MM0N/RSTAT2/RT ( 18 , 20) 

C0MM0N/RSTAT3/RSLPT ( 18 ) , RINTCT ( 18 ) 

C 0MM0N /RS  TAT4/RA (18 , 20) 
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C OMMON /RS  TAT  5 /RO S ( 1 8 , 20) 

CHARACTER*30  RFLNMl 
CHARACTER* 2 5 RSATT , RTPERF 
CHARACTER* 2 RCNUM 
C 
C 

INCLUDE  ' [ASGPROGJCOMVAR.LIS’ 

C 

C Begin  output  to  data  file  with  title  of  results 
C 

DFLAG  - 0 

RECFI 1— ' REC ' //TSTCDNUM ( 1 : 4) // ' .SUM' 

RECFI2-,REC,//TSTCDNUM(1:4)//' .001' 

RECFI 3— ' REC ' //TSTCDNUM (1:4)//' .002' 

RECFI4- 1 REC e //TSTCDNUM (1 : 4) // • . 003 ' 

RECFI 5- ' REC ' //TSTCDNUM (1 : 4) // • . 004 ' 

OPEN ( 3 , FILE-RECFI1 , STATUS- • NEW ' , ERR-3  200 ) 

OPEN ( 7 , FILE-RECFI2 , STATUS- ' NEW ' , ERR-3200 ) 

OPEN (8 , FILE-RECFI3 , STATUS-' NEW' , ERR-3200) 

OPEN ( 9 , FILE-RECFI4 , STATUS- ' NEW ’ , ERR-3200 ) 

OPEN (10 , FILE-RECFI5 , STATUS-' NEW , ERR-3200) 

WRITE (3 , 10) 

10  FORMAT ( 43X, '**  RECTANGLE  WAVEFORM  PROCESSOR  DATA  SUMMARY  **'// 

1 IX, 'A  VALUE  OF  +/-  10000  INDICATES  THAT  THE  QUANTITY  ' , 

2 'COULD  NOT  BE  DETERMINED ' /IX , ' THE  * SYMBOL  INDICATES  THAT' 

3 THE  ALLOWABLE  LIMIT  HAS  BEEN  EXCEEDED //IX, ' TIME  ’, 

4 'QUANTITIES  EXPRESSED  IN  MILLISECONDS '/IX, 'AMPLITUDE  ', 

5 'QUANTITIES  EXPRESSED  AS  PERCENT  OF  FULL  SCALE') 

CALL  DATE (DAY) 

CALL  TIME  (HOUR) 

WRITE (3, 12)  DAY, HOUR 
12  FORMAT (//,50X,A9,12X,A9) 

DO  20  I— 1 , ICHAN 
DO  200  J-1,26 
DEV(J)— 1.0E+4 
200  CONTINUE 

DO  300  J-1,6 
EFLAG(J , I)— 0 
300  CONTINUE 

WRITE (3, 11)  RFLNMl (I) ,RCNUM(I) 

11  FORMAT (//IX, 'RESULTS  FOR  FILE  : ', A3 0 , 5X, ' CHANNEL  NO.  ' ,A2) 

WRITE (3 , 13)  RTPERF ( I ) 

13  FORMAT (/IX, 'FACILITY  : ' ,A25) 

DO  30  K -1,10 
L-(K-1)*2+1 
M-(K-l)*2+2 
0 ( I , L)  — RPOS  ( I , K) 

0 ( I , M) — RNOS ( I , K) 

A(I ,L)— RSSPAV (I ,K) 

A(I,M)-RSSNAV(I,K) 

T ( I , M)  — RTU  ( I , K) 

T(I,L)-RTL(I,K) 

30  CONTINUE 
C 
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Perform  least  squares  fit  on  timeshift  data  for  each  channel  processed 


N-0 

SUMXSQ-0 
SUMXY-0 
SUMY-0 
SUiffi-0 

DO  40  J-1,20,2 

IF(ABS(T(I,J)) .EQ. (1.0E+4))  GOTO  40 
X-(J-1)/2*10 
55  SUMX-SUMX+X 

SUMY— SUMY+T ( I , J ) 

SUMXY-SUMXY+X+T ( I , J ) 

SUMXSQ— SUMXSQ+X'**2 
N-N+l 

40  CONTINUE 
DO  41  J-2,20,2 

IF(ABS(T(X , J) ) .EQ. (1.0E+4))  GOTO  41 
X— (J/2-l)*10+6 
SUMX-SUMX+X 
SUMY— SUMY+T ( I , J ) 

SUMXY— SUMXY+X+T  ( I , J ) 

SUMXSQ-SUMXSQ+X**2 
N-N+l 

41  CONTINUE 

calculate  slope  and  intercept  for  least  squares  fit 

SLOPE-  ( SUMX*SUMY=N*SUMXY)  / ( SUMX**2  -N+SUMXSQ) 

INTCP- ( SUMY - SLOPE*SUMX) /N 
WRITE (3 , 1000) 

1000  FORMAT ( /IX, '1.  TIME  DEVIATION  FROM  THEORETICAL  TIME 
* (ALLOWABLE  LIMIT  - +/-  1.0  MSEC)’) 

DO  49  J-1,20 

49  IND(J)— ' • 

DO  50  J-1,20,2 

IF(ABS(T(I , J) ) .EQ. (1.0E+4) ) GOTO  50 
X-(J-1)/2*10 

75  DEV ( J ) — T ( I , J ) - SLOPE*X- INTCP 

IF(ABS(DEV(J) ) .GT. (1.0))  EFLAG(1,I)-1 
IF(ABS(DEV(J)).GT.(1.0))  IND(J)-'  *' 

50  CONTINUE 

DO  51  J-2,20,2 

IF(ABS(T(I , J) ) .EQ. (1.0E+4) ) GOTO  51 

X-(J/2-l)*10+6 

DEV ( J ) -T ( I , J ) - SLOPE*X= INTCP 

IF(ABS (T(I , J) ) .GT. (1.0))  EFLAG(1,I)-1 

IF(ABS(T(I,J)).GT.(1.0))  IND(J)-'  *e 

51  CONTINUE 
WRITE (3, 1021) 

WRITE(3 , 1002) 

1002  FORMAT ( 7X , ' 1 ' , 9X, ' 2 ' , 9X, ' 3 ' , 9X, ' 4' , 9X, ’ 5 • , 9X, ' 6 ’ , 9X, ' 7 ' , 9X, 

1 '8' ,9X, '9' ,8X, '10') 

1003  FORMAT (6X, 'll' ,8X, '12' ,8X, '13' ,8X, '14' ,8X, '15' ,8X, '16' ,8X, '17' , 
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1 8X, '18' ,8X, '19' ,8X, '20' ) 

WRITE(3 , 1004) (T(I,J) , J-1,10) 

WRITE (7 , 1004) (T(I , J) , J-1,10) 

1004  FORMAT (2X, 13 (F10 . 3) ) 

WRITE (3 , 1005) (IND(J) , J-1,10) 

1005  FORMAT (2X, 13 (4X,1A2,4X)) 

WRITE (3, 1021) 

WRITE (3, 1003) 

WRITE (3 , 1004) (T(I , J) , J-11,20) 

WRITE (7 , 1004) (T(I , J) , J-11,20) 

DO  2004  J-1,20 
RT(I,J)-T(I,J) 

2004  CONTINUE 

WRITE (3, 1005 )(IND(J) , J-11,20) 

DO  52  J-1,20 
52  IND(J)— ' ' 

WRITE (3, 1010) 

1010  FORMAT (/IX, '2.  TIME  LINEARITY  (ALLOWABLE  LIMIT  - +/-  1 %) ' , 

1 2 OX, ' 3 . TIME  OFFSET  (ALLOWABLE  LIMIT  - +/“  l.°  MSEC)') 

time  deviation/slope  intercept  comparison  to  allowable  limits 

IF(ABS( SLOPE) .GT. (0.01))  EFLAG(2,I)-1 
IF(ABS(SLOPE) .GT. (0.01))  IND(l)-'  *’ 

IF(ABS(INTCP) .GT. (1.0))  EFLAG(3,I)-1 
IF(ABS(INTCP) .GT. (1.0) ) IND(2)-'  *' 

SLOPE— SLOPE*100 

WRITE (3, 1011)  SLOPE, IND(l) , INTCP , IND ( 2 ) 

WRITE (8, 1012)  SLOPE, INTCP 
RSLPT ( I ) —SLOPE 
RINTCT(I)— INTCP 
1012  FORMAT ( 2X, F6 . 2 , 5X, F6 . 2) 

1011  FORMAT (/10X, 'SLOPE  - ' , F6 . 2 , 1A2 ,48X, ' INTERCEPT  - ' ,F6.2,lA2) 

IND(l)— ' ’ 

IND(2)-'  ' 

WRITE (3, 1020) 

1020  FORMAT (/IX, '4.  STEADY-STATE  AMPLITUDE  DEVIATION 
1 FROM  THEORETICAL  AMPLITUDE 
1 (ALLOWABLE  LIMIT  - 2.5  %) ' ) 

DO  95  J-1,20 


attenuation  comparison  to  allowable  limits 

IF(ABS(A(I , J) ) .EQ. (1.0E+4))  GOTO  95 
IF(ABS (A(I , J) ) .GT. (2.5))  EFLAG(4,I)-1 
IF(ABS(A(I,J)).GT. (2.5))  IND(J)-'  *' 

95  CONTINUE 

WRITE (3, 1021) 

1021  FORMAT (/59X, 'INTERVAL  NUMBER'/) 

WRITE (3, 1002) 

WRITE (3 , 1004)  (A(I,J) , J-1,10) 

WRITE(9 , 1004)  (A(I,J) , J-1,10) 

WRITE(3 , 1005) (IND(J) , J-1,10) 

WRITE (3, 1021) 
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WRITE (3, 1003) 

WRITE(3 , 1004)  (A(I,J),J-11,20) 

WRITE(9 , 1004)  (A(I,J) ,J— 11,20) 

DO  3004  J-1,20 
RA ( I , J ) —A ( I , J ) 

3004  CONTINUE 

WRITE (3 , 1005)  (XNB(J) ,J— 11,20) 

DO  96  J-1,25 
96  IND(J)-'  1 
WRITE (3, 1030) 

1030  FORMAT (/IX, ’5.  AMPLITUDE  OVERSHOOT  RELATIVE  TO  CALCULATED 
1 STEADY- STATE  AMPLITUDE  (ALLOWABLE  LIMIT  +/-  I3  %) ' ) 


overshoot  comparison  to  allowable  limits 
DO  85  J-1,20 

IF(ABS(0(I,J)) .EQ.1.0E+4)  GOTO  85 
IF(ABS (0(1, J) ) .GT. 13.0)  EFLAG ( 5 , 1 ) — 1 
IF(ABS(0(I,J)).GT.13.0)  IND(J)-'  *' 

85  CONTINUE 

WRITE(3 , 1021) 

WRITE (3, 1002) 

WRITE (3 , 1004)  (0(I,J),J-1,10) 

WRITE (10, 1004)  (0(1, J) ,J-1,10) 

WRITE(3 , 1005)  (XND(J) ,J-1,10) 

WRITE(3 , 1021) 

WRITE (3, 1003) 

WRITE(3 , 1004)  (O(I,J),J-ll,20) 

WRITE (10, 1004)  (0(I,J),J-11,20) 

DO  5004  J-1,20 
ROS (I , J)— 0(1 , J) 

5004  CONTINUE 

WRITE (3 , 1005)  (IND(J) ,J— 11,20) 

20  CONTINUE 

WRITE (3, 1040) 

1040  FORMAT (/IX, ' 6 . CHANNEL -TO -CHANNEL  TIME  DIFFERENCE  - ' , 

1 'ALL  CHANNELS  ARE  COMPARED  TO  THE  FIRST  CHANNEL  PROCESSED' 
1 ' (ALLOWABLE  LIMIT  - +/-  0.25  MSEC)') 

Kl-1 


channel -to -channel  time  variation  to  allowable  limits 

DO  110  I— 1 , ICHAN 
DO  101  J-1,20 
DEV(J)— 1.0E+4 
101  IND(J)— ' ' 

WRITE (3, 11)  RFLNMl(I) ,RCNUM(I) 

DO  100  J-1,20 

IF (ABS (T ( I , J ) ) . EQ . 1 . OE+4)  GOTO  100 
DEV(J)— T(I , J) -T(K1 , J) 

IF (ABS (DEV (J)) .GT.0.25)  EFLAG(6,I)-1 
IF(ABS (DEV( J) ) . GT . 0 . 25)  IND(J)-'  *' 

100  CONTINUE 

WRITE (3, 1021) 
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WRITE (3, 1002) 

WRITE (3, 1004) 

WRITE (3 , 1005) 

WRITE (3, 1021) 

WRITE (3, 1003) 

WRITE (3, 1004) 

WRITE(3 , 1005) 

CONTINUE 
WRITE(3 , 1060) 

FORMAT (//5X, ' FILE  NAME' , 19X, ' CHANNEL  NO.  ' , 20X, ' STATUS ' ) 

NA-0 

NU-0 


(DEV(J) ,J-1,10) 
(IND(J) ,J-1,10) 


(DEV(J) , J-11,20) 
(IND(J) , J-11,20) 


DO  60  I— 1 , I CHAN 
KND— 0 

DO  70  J-1,6 

IF(EFLAG(J,I) .EQ.O)  GOTO  70 
KND-1 

70  CONTINUE 

IF (KND. EQ.O)  NA-NA+1 
IF(KND.EQ.l)  NU-NU+1 

IF(KND.EQ.O)  WRITE(3 , 1070)  RFLNMl(I) ,RCNUM(I) 
IF(KND.EQ.l)  WRITE(3 , 1080)  RFLNM1 ( I ) , RCNUM ( I ) 

1070  FORMAT ( IX, A30 , 10X, A2 , 25X, 'PASS’ ) 

1080  FORMAT ( IX, A30 , 10X, A2 , 25X, 'FAIL') 

60  CONTINUE 

WRITE(3 , 1050)  ICHAN,NA,NU 

1050  FORMAT (//IX, 'TOTAL  NO.  OF  CHANNELS  PROCESSED  - ',12 

1 /IX, 'NUMBER  OF  ACCEPTABLE  CHANNELS  - ' ,12 

2 /IX, 'NUMBER  OF  UNACCEPTABLE  CHANNELS  - ’,12) 
CLOSE(3) 

CLOSE(7) 

CLOSE(8) 

CLOSE(9) 

CLOSE(IO) 

GOTO  3400 

3200  WRITE  (6  ,*)  'ERROR!!!  *"RECTAN. SUM"  CANNOT  BE  ACCESSED*' 
DFLAG  - 1 
3400  CONTINUE 
RETURN 
END 
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The  fifth  segment,  "STR",  contains  the  following 


modules : 

STAIR 

ZEROS 

COMSTR 

G0MPR2 

The  source  code  listing  for  these  modules  is: 
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c 

c 

c 


c 


cccccccccccccccccccccccccc 


cccccccccccccccccccccccccc 


c 


c 


SUBROUTINE  STAIR  c 


c 


c 


C 

C Subroutine  stair  receives  staircase  waveform  data  as  input  and  determines 
C the  leading  edges  for  each  level.  The  point  number  of  each  leading  edge 
C is  stored  in  the  zerinc(I)  array.  There  are  12  values  of  zerinc,  one 
C for  each  new  level  of  amplitude  encountered.  Zerinc (1)  is  the  edge 
C with  a theoretical  time  of  occurrence  of  10  msec.  Zerinc (12)  is  the 
C is  the  edge  with  a theoretical  time  of  occurrence  of  120  millisec  after 
C the  start  of  the  staircase  waveform. 

C 

C Subroutine  stair  calls  subroutine  zer03  to  determine  the  steady -state 
C and  maximum  overshoot  values  for  each  level  of  amplitude.  Subroutine 
C stair  sends  an  address  in  the  Y array  (Y(LL)),  the  number  of  points 
C between  two  consecutive  leading  edges  (NPT) , and  the  number  of  points 
C per  millisec  (NPMS)  as  input  to  subroutine  zer03.  Subroutine  zer03 
C sends  as  output  the  steady-state  value  for  each  level  (SSPAV(I))  and 
C the  maximum  overshoot  for  each  level  (POSOV(I)). 

C 

C Subroutine  stair  upon  receiving  the  steady  state  (SSPAV)  and  overshoot 
C (POSOV)  computes  a timeshift  (TL(I))  which  equals  the  actual  time 
C recorded  minus  the  theoretical  time  for  each  respective  amplitude 
C leading  edge  within  the  channel. 

C 

C Subroutine  stair  then  computes  the  normalized  steady- state  and  normalized 
C overshoot  values  and  stores  these  values  in  place  of  their  respective 
C actual  values  in  the  SSPAV  and  POSOV  arrays.  The  normalized  steady - 
C state  values  are  the  STEADY  STATE  ERROR  divided  by  the  full-scale 
C theoretical  amplitude  multiplied  by  100  to  give  a percentage  of  full 

C scale.  The  normalized  overshoot  values  are  the  OVERSHOOT  beyond  the 

C calculated  steady- state  value  divided  by  the  change  in  calculated 
C amplitude  multiplied  by  100  to  give  a percentage  of  actual  change  in 
C amplitude  recorded. 

C 

C Subroutine  stair  then  writes  the  results  to  a file.  Outputting  normalized 

C overshoot  (POSOV(I)),  normalized  steady  state  (SSPAV(I))  and  timeshift 

C (TL(I))  for  each  level. 

C 

C Principal  variables: 

C TL(I)  - timeshift  or  time  error  for  leading  edges  within  a channel 

C POSOV(I)  — overshoot  for  each  level  of  amplitude  encountered 

C SSPAV(I)  — steady-state  amplitude 

C NPMS  - number  of  points  per  msec 
C G — full-scale  output 

C NOPTS  - number  of  points  per  data  curve 

C ZERINC (I)  - point  number  corresponding  to  a leading  edge 

C LEVEL  ==  amplitude  of  a particular  portion  of  waveform 

C NPT  - number  of  points  between  leading  edges 

C Y(LL)  - the  address  in  the  Y array  that  is  passed  to  subroutine  zer03 

C that  defines  the  first  element  in  the  x array  in  subroutine 
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C zer03.  The  first  element  is  the  first  point  on  any  individual 

C level. 

C 

C 

SUBROUTINE  STAIR ( I CHAN , DFLAG , CNTFL) 


Define  data  variables  and  arrays 

INTEGER  NOPTS , LL , LU , NPT , TB , TR , DFLAG , CNTFL 
INTEGER  Nl(13) ,N2(13) ,N3(13) ,N4(13) 
INTEGER  ZIC , ZERINC(22) , 

* LB (14) ,TZIC,UBND,LBND, START, FINISH, 

* INGOV(13),IPSOV(13) 

REAL  XMIN , XMAX , YHIN , YMAX , TO , G , LEVEL , N9 
REAL  Y( -1000  : 10000) 

REAL  NPMS  , P0S0V ( 23 ) , YHN , YMX , PAVGOV , NAVGOV , 

* SSPAV(23) ,TL(23) ,TU(13) 

INCLUDE  ' DISK$0VR: [ASGPROG] GOMVAR. LIS ' 
CHARACTER  IFLNM1*30 , CHANUM  *2, DUMMY  *1 
CHARACTER  *30  CHAFIL(108) 

INTEGER  CHANUMB 


Define  common  block  symbolic  names 


COMMON  /FILE/ 

IFLNM1 

* 

/CBNPTS/  NOPTS 

* 

/OVER/ 

POSOV , NEGOV 

* 

/RANGE/ 

FINISH 

* 

/SHFT/ 

TL.TU 

* 

/STPT/ 

N1,N2,N3,N4 

* 

/PTN/ 

START, LB 

* 

/COUNT/ 

ZIC , ZERINC 

* 

/FLAG/ 

EFLAG 

* 

/CBFIL/  ' 

CHAFIL 

* 

/CBCH  / CHANUMB 

C0MM0N/SSTAT1/SP0S(18 , 13) , SSSPAV(18 , 13) , STL(18 , 13) , SFLNM1(18) , 
& SSATT(18) ,SCNUM(18) ,STPERF(18) 

CHARACTER*30  SFLNM1 
CHARACTER*2 5 SSATT, STPERF 
CHARACTER*2  SCNUM 
DIMENSION  AVE(23) 


stair  waveform  processor 

DFALG  - 0 
WRITE (6 , 1) 

1 FORMAT ( /IX, ***  STAIR  WAVEFORM  PROCESSOR  **' ) 

I CHAN  - 0 
KNUM  - CHANUMB 
IKJ  - -2 
26  CONTINUE 

DO  8700  1—1000,10000 
Y(I)  - 0.0 
8700  CONTINUE 
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IF  (KNUM  .EQ.  0 .AND. 
IKJ  - IKJ  +6 
YMAX  - 200. 

YMIN  - -200. 


CNTFL  .EQ.  1)  GOTO  300 


Initialize  variables 

DO  22  1-1,13 
TL(I)— 1 . 0E+4 
POSOV(I)— 1 . 0E+4 
22  SSPAV(I)— 1 . 0E+4 


Input  waveform  data 

5000  IFLNM1  - CHAFIL(IKJ) 

READING  FROM  DATA  FILE  ARRAY 

8 INDX-1+INT( FLOAT ( IKJ )/6.0) 

CALL  GETWAVE ( INDX ,0,3, NFP , NLP , DEL , ' O ' ,Y) 

6 I CHAN  - I CHAN  + 1 
KNUM  - KNUM  - 1 
NPMS-1 . 0 /DEL/1000 . 

C 

C store  edge  locations  in  zerinc  array 
C 

C when  waveform  reaches  1/2  the  expected  amplitude  for  each  respective 
C level,  it  is  that  point  that  determines  the  leading  edge. 

C 

ZIC-0 
G— 200 

LEVEL-G/10 

NOPTS  - (NLP  - NFP)  + 1 
DO  10  J-NFP, NOPTS 

IF(ZIC.GE.5)  GOTO  20 
IF (Y( J) .LT. LEVEL)  GOTO  10 
ZIC-ZIC+1 
ZERINC (ZXC)-J 
LEVEL-LEVEL+G/ 5 
GOTO  10 

20  IF(ZIC.GT.5)  GOTO  30 
LEVEL-G/2 

IF(Y(J) .GT. LEVEL)  GOTO  10 
ZIC-ZIC+1 
ZERINC (ZIC)-J 
GOTO  10 

30  IF(ZIC.GT. 6)  GOTO  40 

LEVEL—  G/2 

IF(Y(J) .GT. LEVEL)  GOTO  10 
ZIC-ZIC+1 
ZERINC (ZIC)-J 
LEVEL— .90*G 
GOTO  10 
40  K-J - 1 
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IF(Y(J) .LT. LEVEL. OR. Y(J) .LT.Y(K))  GOTO  10 
ZIC-ZIC+1 
ZERINC(ZIC)— J 
LEVEL-LEVEL+G/5 
10  CONTINUE 
C 

C call  subroutine  to  calculate  steady- state  and 
C overshoot  values  for  each  level 
C 

C initial  zero  level 
C 

LL-1 

LU— ZERING(l) 

NPT-LU-LL+l 

CALL  ZER03(Y(LL) ,NPT, SSPAV(l) , POSOV(l) , 0 ,NPMS) 


positive  levels 


DO  90  1-2 , 6 
J-I-l 

LL-ZERING(J) 

LU-ZERINC(I) 

NPT— LU-LL+1 

CALL  ZER03 (Y(LL) ,NPT,SSPAV(I) ,P0S0V(I) ,0,NPMS) 
90  CONTINUE 


intermediate  levels 


DO  11  1-7,8 
J-I-l 

LL-ZERINC (J ) 

LU— ZERINC ( I ) 

NPT-LU~LL+1 

GALL  ZER03 (Y(LL) ,NPT, SSPAV(I) , POSOV(I) , 1 ,NPMS) 
11  CONTINUE 


negative  levels 

DO  15  1-9,12 
J-I-l 

LL-ZERINC (J) 

LU-ZERINC(I) 

NPT— LU-LL+1 

CALL  ZER03 (Y(LL) ,NPT , SSPAV(I) , POSOV(I) , 0 ,NPMS) 
15  CONTINUE 

final  zero  level 

1-13 

LL-ZERINC (12) 

LU-NOPTS 
NPT— LU - LL+1 

CALL  ZERO 3 (Y(LL) ,NPT , SSPAV(I) , POSOV(I) , 0 ,NPMS) 
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compute  timeshift  using  points  on  the  leading  edges 

DO  110  1-1,12 
DAN— ZERINC ( I ) - 1 
TL(I)— ((DAN) /NPMS- 10*1) 

110  CONTINUE 

normalize  overshoot  and  steady- state  error  values 

AVE  is-  the  value  against  which  the  overshoot  is  normalized. 
Since  there  is  no  amplitude  change  before  the  initial  zero 
level  (theoretically)  the  value  used  for  normalization 
has  been  set  to  full-scale,  steady-state  value  (SSPAV(6)). 
DO  199  1-1,13 
IF  (I.EQ.l)  THEN 
AVE(I)— SSPAV (6) 

ELSE 

AVE(I)— ABS(SSPAV(I) -SSPAV(I-l) ) 

END  IF 
199  CONTINUE 

DO  200  1-1,13 


sspav  is  redefined  as  deviation  in  percent  of  full  scale 

POSOV ( I )— ABS  (POSOV ( I ) .-  SSPAV ( I ) ) /AVE ( I ) *100 . 0 
205  LEVEL- ( I - 1 ) *40 

IF(I.EQ.7)  LEVEL-LEVEL- 240 
IF(I.GE.8)  LEVEL-LEVEL-480 
SSPAV(I)-( (SSPAV(I) -LEVEL) /200.0)*100„0 
IF(I.GT.7.AND.I.LT.13)  SSPAV(I)— SSPAV(I) 

200  CONTINUE 


store  results 

I-INDEX(IFLNM1, ' . ' ) 

910  CHANUM(1 : 1)  - IFLNMl(I+2 : 1+2) 

CHANUM(2 : 2)  - IFLNMl(I+3 : 1+3) 

SFLNMl(INDX)— IFLNM1 
S S ATT ( INDX) -S  ENATT 
SCNUM ( INDX) -CHANUM 
STPERF ( INDX) -TSTPRF 
DO  210  K-1,13 

SPOS ( INDX , K) -POSOV (K) 

SSSPAV(INDX,K)— SSPAV(K) 

STL(INDX,K)-TL(K) 

210  CONTINUE 
CLOSE (1) 

WRITE (6,*) 'PROCESSING  COMPLETED  FOR  FILE  ' , IFLNM1 
GOTO  26 
300  CONTINUE 
RETURN 
END 

SUBROUTINE  ZERO 3 (X , NPTS , AVG , PK , IFL , XNPM) 
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DIMENSION  X(1000) ,Y(1000) 

K-NPTS -XNPM 
IF(IFL.EQ.l)  GOTO  520 
PK—  l.E+4 
DO  500  I-1,K 
500  PK-AMAX1(X(I) ,PK) 

GOTO  510 
520  PK-l.E+4 

DO  530  I-l.K 
530  PK-AMIN1(X(I) ,PK) 

510  K-NPTS -XNPM 
LEA— NPTS/2 
SUM-0 „ 

DO  747  LP— LEASK 
SUM— SUM+X ( LP ) 

747  CONTINUE 

AV  G— SUM/ ( K - LEA+1 ) 

RETURN 

END 

C 

G ccccccccccccccccccccccccccc 

G c c 

C c SUBROUTINE  COMSTAIR  c 

G c c 

G ccccccccccccccccccccccccccc 

G 

G Subroutine  Comstair  performs  a comparison  between  processed  stair 
C waveform  data  and  the  corresponding  theoretical  values  desired. 

C 

C Output  of  +/-  10000  indicates  quantity  could  not  be  determined 
C 

C The  symbol  * indicates  the  allowable  limit  has  been  exceeded 
C 

C Time  quantities  are  expressed  in  millisec 
C 

C Amplitude  quantities  are  expressed  as  percent  of  full  scale 

C 

G Subroutine  Comstair  receives  data  from  program  MGAPR0C 

G 

C NCHAN  - number  of  channels 
C C(X,J)  ” channel  number 
C T(I)  - timeshift 
C EFLAG  - channel  status 
C 

SUBROUTINE  COMSTR ( T , NCHAN , EFLAG , DFLAG ) 

C 

C define  data  types 
C 

CHARACTER  IFLNM1  *14 
CHARACTER  DAY*9,H0UR  *9 

CHARACTER  *30  STARFI1 , STARFI2 , STARFI3 , STARFI4 , STARFI5 , STARFI6 

CHARACTER  *4  TSTCDNUM 

INTEGER  EFLAG (9 , 18) , IND(23) , DFALG 
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REAL  0(18,23) , A(18 , 23) , T(18 , 23) , SLOPE , INTCP , DEV(23) , SLOPEA, INTCPA 
INCLUDE  ' DISK$OVR: [ASGPROG] COMVAR. LIS ' 

C 

C COMMON  BLOCKS 
C 

COMMON  /CBCDNUM/TSTCDNUM 

COMMON/S STAT1/SP0S (18 , 13) , SSSPAV(18 , 13) , STL (18 , 13) , SFLNM1(18) , 

& SSATT(18) ,SCNUM(18) ,STPERF(18) 

COMMON/SSTAT2/ST(18 , 13) 

COMMON/SSTAT3/SSLPT(18) ,SINTCT(18) 

COMMON/S  STAT4/S  SLPA ( 1 8 ) , S INTCA (18) 

COMMON/SSTAT5/SA(18 , 13) 

C0MM0N/SSTAT6/S0S (18 , 13) 

CHARACTER* 30  SFLNMl 
CHARACTER*25  SSATT , STPERF 
CHARACTER* 2 SCNUM 
C 

C compare  results  of  stair  waveform  processing  with  allowable  limits 
C 

DFLAG  - 0 

STARFI1-' STAIR* //TSTCDNUM// ' .SUM' 

STARFI2-' STAIR' //TSTCDNUM//' .001' 

STARFI3-' STAIR' //TSTCDNUM//’ .002' 

STARFI4— • STAIR ' //TSTCDNUM// ' . 003 1 
STARFI5-’ STAIR’ //TSTCDNUM//' .004' 

STARFI6-’ STAIR’ //TSTCDNUM//' .005' 

OPEN (3 , FILE-STARFI1 , STATUS-' NEW' , SHARED , ERR-3200) 

OPEN (7 , FILE-STARFI2 , STATUS-' NEW' , SHARED , ERR-3200) 

OPEN (8 , FILE-STARFI3 , STATUS-' NEW' , SHARED , ERR-3200) 

OPEN ( 9 , FILE— STARFI4 , STATUS- ’ NEW ’ , SHARED , ERR-3200) 

OPEN (10 , FILE-STARFI5 , STATUS- 'NEW' , SHARED , ERR-3200) 

OPEN (11 , FILE— STARFI6 , STATUS-' NEW’ , SHARED , ERR-3200) 

WRITE (3 , 10) 

10  FORMAT ( 43X, '**  STAIR  WAVEFORM  PROCESSOR  DATA  SUMMARY  **'// 

1 IX, 'A  VALUE  OF  +/"  10000  INDICATES  THAT  THE  QUANTITY  ’ , 

2 'COULD  NOT  BE  DETERMINED ' /IX , ’ THE  * SYMBOL  INDICATES  THAT' 

3 , ' THE  ALLOWABLE  LIMIT  HAS  BEEN  EXCEEDED . ' //IX, ’ TIME  ' , 

4 'QUANTITIES  EXPRESSED  IN  MILLISECONDS '/IX, 'AMPLITUDE  ', 

5 'QUANTITIES  EXPRESSED  AS  PERCENT  OF  FULL  SCALE') 

CALL  DATE (DAY) 

CALL  TIME (HOUR) 

WRITE (3, 12)  DAY, HOUR 
12  FORMAT (//,50X,A9,10X,A9) 

C 

C input  data 
C 

C 0(1 ,J)  - overshoot 
C A(I , J)  - steady-state  error 
C T(I,J)  - timeshift 
C 

IFLNM1  - ' STAIR . DAT ' 

DO  20  I-l.NCHAN 
DO  200  J-1,13 
DEV ( J)— 1 . OE+4 
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200  CONTINUE 

DO  300  J-1,9 
EFLAG(J,I)=0 
300  CONTINUE 

WRITE(3 , 11)  SFLNMl ( I ) ,SCNUM(I) 

11  FORMAT (//IX  / RESULTS  FOR  FILE  : \ A3 0 , 5X, ' CHANNEL  NO.  ' ,A2) 

WRITE (3, 13)  STPERF(I) 

13  FORMAT (/IX, ’FACILITY  : ' ,A25) 

DO  30  K -1,13 
0(1 ,K)— SPOS (I , K) 

A(I,K)=SSSPAV(I,K) 

T(I,K)=STL(I,K) 

30  CONTINUE 
C 

G perform  least  squares  using  timeshift  data  on  y-axis  and  theoretical 
C time  on  x-axis 
C 

C if  data  correlation  is  perfect  slope  * 0.0  and  y intercept  =0.0 

C 

N-0 

SUMXSQ-Q 

SUMXY-0 

SUMY-0 

SUMX-0 

DO  40  J-1,12  . 

IF(ABS (T(I , J)  ).o EQ.  (1 . QE+4)  ) GOTO  40 
X-J*10 

55  SUMX-SUMX+X 

SUMY~SUMY+T(I,J) 

SUMXY=SUMXY+X*T (I , J) 

SUMXSQ-SUMXSQ+X**2 

N-N+l 

40  CONTINUE 

SLOPE- ( SUMX*SUMY-N*SUMXY) / ( SUMX**2 - N*SUMXSQ) 

INTCP-C  SUMY-  SLOPE*SUMX)  /N 
WRITE (3, 1000) 

1000  FORMAT (/IX, '1.  TIME  DEVIATION  FROM  THEORETICAL  TIME 
* (ALLOWABLE  LIMIT  +/"  1.0  MSEC)') 

DO  49  J-1,13 
49  IND(J)-’  ! 

DO  50  J-1,12 

IF(ABS(T(I , J) ) .EQ. (1.0E+4))  GOTO  50 
X~J*10 
C 

C DEV(I)  - deviation  from  least  squares  fit 
C IND(I)  - indicated  channel  status 

C 

75  DEV ( J ) «T ( I , J ) “ SLOPE*X- INTCP 
C 

C EFLAG(L,I)  = status  indicator  for  a particular  comparison 
C L = test  being  compared 
C I = channel  number 

C there  are  eight  test  comparisons  for  each  channel 
C 


c 

c 

c 

c 

c 

c 

c 

c 

c 


EFLAG  ( 1 , J ) 
" (2 , J) 

" (3 , J) 

" (4,J) 

" (5 , J) 

" (6,J) 

" (7,J) 

" (8,J) 


time  deviation  from  least  squares  fit 
time  slope  deviation  from  least  squares  fit 
time  intercept  deviation  from  least  squares  fit 
amplitude  deviation  from  least  squares  fit 
amplitude  slope  deviation  from  least  squares  fit 
amplitude  intercept  deviation  from  least  squares  fit 
amplitude  deviation  from  theoretical 
amplitude  overshoot 


IF(ABS(T(I, J)) .GT. (1.0))  EFLAG(1,I)-1 
IF(ABS(T(I,J)).GT.1.0)  IND(J)-’  *' 

50  CONTINUE 

WRITE(3 , 1021) 

WRITE(3 , 1002) 

1002  FORMAT ( 7X , ' 1 ' , 9X, ' 2 * , 9X, ’ 3 • , 9XS ’ 4! , 9X, ’ 5 ' , 9X, ’ 6 ' , 9X , ' 7 ' , 9X, 

1 '8' ,9X, '9' ,8X, '10' ,8X, 'll' ,8X, '12' ) 

WRITE (3 , 1004) (T(I , J) , J-1,12) 

WRITE (7 , 1004) (T(I , J) , J-1,12) 

1004  FORMAT ( 2X , 12 ( F10 . 3 ) ) 

DO  2004  J-1,12 

ST(I , J)— T(I , J) 

2004  CONTINUE 

WRITE (3 , 1005) (IND(J) , J-1,12) 

1005  FORMAT (2X, 13 (4X,1A2,4X)) 

DO  52  J-1,13 

DEV(J)— 1 . OE+4 
52  IND(J)— ' ' 

WRITE (3, 1010) 

1010  FORMAT (/IX, '2.  TIME  LINEARITY  (ALLOWABLE  LIMIT  - +/•  1 %)', 

1 20X, '3.  TIME  OFFSET  (ALLOWABLE  LIMIT  - +/-  1-0  MSEC)') 

IF(ABS (SLOPE) .GT. (0.01))  EFLAG(2,I)-1 
IF(ABS (SLOPE) .GT. (0.01))  IND(l)-'  *' 

IF(ABS(INTCP) .GT.(l.O))  EFLAG(3,I)-1 
IF(ABS(INTCP) .GT. (1.0))  IND(2)-’  *’ 

SLOPE-SLOPE*100 

WRITE(3 , 1011)  SLOPE, IND(l) , INTCP , IND ( 2 ) 

WRITE (8, 1009)  SLOPE, INTCP 
SSLPT(I)— SLOPE 
SINTCT(I)— INTCP 
1009  FORMAT (1X,F6.2,2X,F6.2) 

1011  FORMAT (/10X, ’SLOPE  - ', F6 . 2 , 1A2 ,48X, ' INTERCEPT  - ’.F6.2.1A2) 

IND(l)— ’ ' 

IND (2)—'  ’ 

C 

C perform  amplitude  least  squares  fit 
C 

C adjust  theoretical  values  to  be  all  positive 
C perform  least  squares  fit  on  adjusted  data 
C correct  y - intercept  value 
C 

N— 0 

SUMXSQ— 0 

SUMXY-0 

SUMY-0 
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SUMX-0 

DO  41  J-2,12 

IF(ABS(A(I,J)) .EQ.1.0E+4)  GOTO  41 

X-(J-l) *20+100 

IF(J . EQ. 7)  X-X-120 

IF(J.GE.8)  X-X-240 

Al— A ( I , J ) 

IF(J.GE.S)  Al— A(I.J) 

SUMX-SUMX+X 
SUMY— SUMY+Al 
SUMXY-SUMXY+X*A1 
SUMXSQ— SUMXSQ+X**2 
N-N+l 

41  CONTINUE 

SLOPEA- ( SUMX*SUMY- N*SUMXY) / ( SUMX**2 - N*SUMXSQ) 
INTCPA- ( SUMY- SLOPEA*SUMX) /N+SLOPEA*10Q 
DO  51  J-2,12 

IF(ABS (A(I , J) ) . EQ. 1 . OE+4)  GOTO  51 
Al— A(X , J) 

IF(J.GE.8)  Al— A(I,J) 

X— (J-l)*20 
IF(J.EQ.7)  X— X-120 
IF(J.GE.8)  X-X-240 


amplitude  deviation  from  least  squares 

DEV ( J ) -Al - SLOPEA*X~ XNTCPA 
IF ( ABS (DEV ( J ) ) „ GT . 2 * 5 ) EFLAG(8,I)-1 
IF(ABS (DEV(J) ) .GT.2.5)  IND(J)-' 

51  CONTINUE 

C 

WRITE (6, 1252) 

1252  FORMAT ( ' STEADY-STATE  AMPLITUDE  DEVIATIONS  FROM  BEST  FIT 
* STRAIGHT  LINE  (%  FULL  SCALE)8) 

WRITE (6 , 1253) (J ,BEV(J) , J-2,12) 

1253  FORMAT (4 ( ' DEV( ' , 12 , 8 )-' , F8 . 3) ) 

C 

1021  FORMAT (/59X, 8 INTERVAL  NUMBER'/) 

C WRITE(3 , 1002) 

C 

C output  amplitude  deviations  and  indicators 
C 

DO  96  J-1,13 
DEV ( J ) — 1 . OE+4 
96  IND(J)-'  5 
WRITE(3 , 1110) 

1110  FORMAT (/IX, *4.  AMPLITUDE  LINEARITY  (ALLOWABLE  LIMIT  - +/-  2.5  %) 8 
1 , 15X, ' 5 . AMPLITUDE  OFFSET  (ALLOWABLE  LIMIT  - +/-  2.5  %) ' ) 

IF (ABS (SLOPEA) .GT. 0.025)  EFLAG(5,X)-1 
IF ( ABS ( SLOPEA) .GT. 0.025)  IND(l)-’  *• 

IF (ABS ( INTCPA) .GT.2.5)  EFLAG(6,I)-1 
IF ( ABS (INTCPA) .GT.2.5)  IND(2)-'  *' 

SLOPEA— SLOPEA*100 

WRITE(3 , 1011)  SLOPEA, IND(l) .INTCPA, IND( 2) 
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WRITE (9, 1009)  SLOPEA, INTCPA 
SSLPA(I)— SLOPEA 
SINTCA( I) -INTCPA 
IND(l)— ' ' 

IND(2)-'  ' 

WRITE (3, 1120) 

1120  FORMAT (/IX, '6.  STEADY-STATE  AMPLITUDE  DEVIATION  FROM  ' 
1 ' THEORETICAL  AMPLITUDE' 

1 ' (ALLOWABLE  LIMIT  - +/-  2.5  %)') 

DO  95  J-1,13 

IF(ABS(A(I , J) ) .EQ.1.0E+4)  GOTO  95 
IF(ABS(A(I,J)) .GT.2.5)  EFLAG(7,I)-1 
IF (ABS (A(I,J)) .GT.2.5)  IND(J)-'  *' 

95  CONTINUE 

WRITE (3, 1021) 


output  steady- state  results 
WRITE (3, 1002) 

WRITE (3, 1004) (A( I, J) , J-l, 12) 

WRITE(10 , 1004) (A(I , J) , J-l , 12) 

DO  3004  J-l ,12 
SA(I,J)-A(I,J) 

3004  CONTINUE 

WRITE (3, 1005 )(IND(J) , J-l, 12) 

DO  97  J-1,13 
DEV(J)— 1.0E+4 
97  IND(J)— ' ' 

WRITE (3, 1030) 

1030  FORMAT (/IX, '7.  AMPLITUDE  OVERSHOOT  RELATIVE  TO  CALCULATED 
1 STEADY-STATE  AMPLITUDE  (ALLOWABLE  LIMIT  +/-  13  %)') 

DO  85  J-2,13 

IF (ABS (0(1 , J) ) . EQ. 1 . OE+4)  GOTO  85 
IF (ABS (0(1, J)) .GT.13.0)  EFLAG(8,I)-1 
IF(ABS(0(I,J)) .GT.13.0)  IND(J)-'  *' 

85  CONTINUE 

WRITE(3 , 1021) 

WRITE (3, 1002) 


output  overshoot  results 

WRITE (3 , 1004)  (0(1, J) , J-2,13) 

WRITE (11, 1004)  (0(1, J) , J-2,13) 

DO  4004  J-2,13 
S0S(I,J)-0(I,J) 

4004  CONTINUE 

WRITE (3 , 1005)  (IND(J) , J-2,13) 

20  CONTINUE 
GOTO  320 
3200  CLOSE(3) 

WRITE (6 , *)  'ERROR  ! ! ! ' 

WRITE (6,*)  ' ***  "STAIR. SUM"  CANNOT  BE  OPENED  ***' 
DFLAG  - 1 
320  CLOSE(7) 
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CL0SE(8) 

CL0SE(9) 

CLOSE (10) 

CLOSE (11) 

RETURN 

END 

C 

G cccccccccccccccccccccccccccccccc 

C c c 

C c SUBROUTINE  C0MPR2  c 

C c c 

C cccccccccccccccccccccccccccccccc 

c 

C Subroutine  Compr2  checks  the  channel- to -channel  time  error. 

C All  channels  are  compared  to  the  first  channel  processed. 

C 

C Input  to  subroutine  Compr2  will  be  the  number  of  channels 
C and  the  timeshift  of  each  level  of  each  channel  processed. 

C Output  will  be  the  channel  status  of  all  channels  processed. 

C 

SUBROUTINE  C0MPR2 ( T , NCHAN , EFLAG ) 
define  data  types 

G0MM0N/SSTAT1/SP0S (18 , 13) , SSSPAV(18  , 13) , STL(18 , 13)  , SFLNM1(18) , 
& SSATT(18) , SGNUM(18) , STPERF(18) 

CHARACTER*  30  SFLNM1 
CHARAGTER*25  SSATT, STPERF 
CHARACTER* 2 SCNUM, IND(23) 

REAL  T(18 , 23) ,DEV(23) 

INTEGER  EFLAG (9, 18) 

WRITE (3 , 1040) 

1040  FORMAT (/IX s ' 8 . CHANNEL - TO - CHANNEL  TIME  DIFFERENCE,' 

1 8 ALL  CHANNELS  ARE  COMPARED  TO  THE  FIRST' 

1 • CHANNEL  PROCESSED8 

1 8 (ALLOWABLE  LIMIT  - +/*  0.25  MSEC)1) 

Kl-1 

DO  110  1-1, NCHAN 
DO  101  J-1,13 
DEV(J)— 1.0E+4 
101  IND(J)-’  ' 

C 

WRITE (3 , 11)  SFLNM1 (I) , SCNUM(I) 

DO  100  J-1,12 

IF(ABS(T(I , J) ) . EQ.1.0E+4)  GOTO  100 

C 

C calculate  channel- to -channel  time  deviation 
C comparing  to  the  first  channel  processed 
C 

DEV(J)— T(I , J) -T(K1 , J) 

IF(ABS(DEV(J)) .GT.0.25)  EFLAG (9 , I)-l 
IF (ABS (DEV( J) ) . GT . 0 . 25)  IND(J)-'  *' 

100  CONTINUE 
C 
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output  deviations 


WRITE(3 , 1021) 

WRITE (3, 1002) 

WRITE (3 , 1004)  (DEV(J) , J— 1 , 12) 

WRITE (3 , 1005)  (IND(J) ,J-1,12) 

110  CONTINUE 

WRITE(3 , 1060) 

1060  FORMAT (//5X, 'FILE  NAME' , 19X, ' CHANNEL  NO 20X, ' STATUS ' ) 
NA-0 
NU-0 

DO  60  I-l.NCHAN 
KND-0 

DO  70  J-1,9 

IF(EFLAG(J,I) .EQ.0)  GOTO  70 
KND-1 

70  CONTINUE 


determine  channel  status 

IF(KND.EQ.O)  NA-NA+1 
IF(KND.EQ.l)  NU-NU+1 

IF(KND.EQ.O)  WRITE(3 , 1070)  SFLNMl(I) , SCNUM(I) 

IF(KND.EQ.l)  WRITE(3 , 1080)  SFLNMl(I) , SCNUM(I) 

70  FORMAT ( IX, A30,10X,A2,25X, 'PASS') 

80  FORMAT ( IX, A30,10X,A2,25X, 'FAIL' ) 

60  CONTINUE 

output  channel  status 

WRITE(3 , 1050)  NCHAN , NA , NU 
1050  FORMAT (//IX, 'TOTAL  NO.  OF  CHANNELS  PROCESSED  - ’,12 

1 /IX, 'NUMBER  OF  ACCEPTABLE  CHANNELS  - ',12 

2 /IX, 'NUMBER  OF  UNACCEPTABLE  CHANNELS  - ' ,12) 

CLOSE(3) 

11  FORMAT (//IX, 'RESULTS  FOR  FILE  : ', A3 0 , 5X, ' CHANNEL  NO.  ' ,A2) 
1021  FORMAT (/59X, 'INTERVAL  NUMBER'/) 

1002  FORMAT ( 7X , * 1 ' , 9X, ' 2 ' , 9X, 9 3 • , 9X, ' 4' , 9X, ' 5 ' , 9X, ' 6 ' , 9X, ' 7 ' , 9X, 
1 '8' ,9X, '9' ,8X, '10' ,8X, 'll' ,8X, '12' ) 

1004  FORMAT (2X, 13 (F10 . 3) ) 

1005  FORMAT (2X, 13 (4X.1A2.4X)) 

RETURN 

END 
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6.  EPROM  CHECK-OUT  SOFTWARE 


This  section  presents  a narrative  description,  flow  diagram,  and  program 
listing  of  the  software  which  has  been  developed  at  MGA  Research  Corporation 
(MGA)  to  compare  the  content  of  an  EPROM  with  the  ideal  test  waveforms. 

The  ideal  test  waveform  data  was  generated  and  stored  in  files  on  the 
VAX  11/780  at  NHTSA.  The  data  files  were  then  transmitted  to  an  IBM  PC  at  MGA 
and  converted  to  a binary  format  suitable  to  be  burned  into  EPROMs.  Due  to 
the  size  of  the  waveform  data  files,  the  waveform  was  divided  into  segments 
such  that  each  segment  would  be  only  as  large  as  the  capacity  of  one  of  the 
EPROMs  in  the  waveform  generator  instrument.  In  total,  there  are  eight  EPROMs 
contained  within  the  waveform  generator  instrument  which  contain  two  sequences 
of  test  waveforms.  The  EPROMs  were  burned  using  an  IBM  PC  and  a PROLOG  EPROM 
programmer.  Once  they  were  burned,  the  contents  of  each  EPROM  was  read  back 
and  stored  in  a data  file.  The  data  file  was  then  compared  against  the  ideal 
test  waveform  data  using  the  EPROM  check-out  program. 

The  EPROM  check-out  program  was  written  in  the  BASIC  language  on  an  IBM 
PC.  The  program  first  establishes  the  required  arrays  for  the  data  and  then 
prompts  the  user  for  the  name  of  the  data  file  which  contains  the  EPROM  data 
(the  data  file  which  was  read  back  from  the  EPROMs).  Once  the  user  has 
responded,  the  program  checks  the  data  read  back  from  the  EPROMs  against  the 
data  files  which  contain  the  ideal  waveform  data.  If  there  are  any  discrep- 
ancies between  the  two  data  sources,  the  user  is  given  a "FAIL"  message 
otherwise  a "PASS"  message  is  given.  Shown  in  Figures  6-1  and  6-2  are  the 
flow  diagram  and  program  listing  respectively,  for  this  program. 
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READ  EPROM  DATA  FILE  NAME 


I 


OPEN  EPROM  DATA  FILE 


I 


K=1 


GIVE 

"FAIL" 

MESSAGE 


GIVE 

"PASS" 

MESSAGE 


OPEN  FILE  CONTAINING  PROPER 
1/8  th  OF  IDEAL  WAVEFORM 


COMPARE  THE  TWO  FILES 


K=K+1 


Figure  6-1  Flow  Diagram  of  EPROM  Check-out  Program 


10  CIS 

15  WIDTH  80 

20  DIM  C$(8)SX3$(16 ) 

30  C$  (1  )=".001" 

35  C$(2  )-"o002" 

40  C$(3)^"oQ03" 

50  C$(4  )s".004" 

55  C$(5)=".,Q05" 

60  C$(6  )~"o006" 

65  C$  (7  )”"c007*s 
70  C$(8 )s,,c008” 

100  INPUT  "PLEASE  ENTER  THE  EPROM  LEVEL  (LOW/HIGH)  tt",C9$ 
110  C10$s"8s"+C9$+".DAT" 

120  OPEN  "l"s#19C10$ 

121  FOR  KS1  TO  8 

130  C11$^C9$+C$(K) 

140  OPEN  "l"s#2,C11$ 

150  FOR  N-1  TO  63 

160  INPUT  #1,X1$ 

167  REM  PRINT  XI  $ 

170  X2$* INPUTS (50, #2) 

175  NN-1 

180  FOR  J=1  TO  16 

190  X3$(J )=MID$(X2$(NN,2) 

200  NN=NN+3 

210  NEXT  J 

220  FOR  J-1  TO  16 

230  X4$sX4$+X3$ ( J ) 

240  NEXT  J 

245  REM  PRINT  XI  $ 

250  IF  (X4$=X1 S )s0  THEN  300 

252  X4$-"" 

255  NEXT  N 

256  GOTO  260 

260  PRINT  »PASS",K 

261  CLOSE  # 2 

262  NEXT  K 

263  CLOSE  m 

264  END 

300  PRINT  "FAIL" 

310  END 


Figure  6-2  EPROM  Check-Out  Program  Listing 


140  copies 
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